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Chapter 1 

Maxwell's Equations and Gauge 
Fields 



1.1 Maxwell's Equations 

describe all phenomena related to (classical) electrical and magnetical fields. In 
free space they are given by: 

AivE = — divB = 

£p (1.1) 

curl E = —B curl B = (j, (j + e E) 

where E denotes the electric and B the magnetic field; g(x, t) is the charge den- 
sity and j{x) current density. The divergence and the rotation of a function R(x) 
is defined by 

3 d 

divR(x) = g^T R i( X ) = 9iRi(x) 

4 = 1 % 

3 d 

(curl Jl(su)) ( = ^2 e ijkQ—Rk = eijkdjRkix) 
j,k=i i 

The notation on the very right-hand side employs Einstein's summation conven- 
tion: summation (here from 1 to 3) is implied for each index that appears twice 
€ijk is the antisymmetric Levi-Civita-Symbol, defined by 

fl (ijA;) = (123),(231),(312) 
tijk = I -1 (ijk) = (132), (321), (213) 
y otherwise 

It follows that e ijk = e jki and e ijk = -e jik . 
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Maxwell's Equations and Gauge Fields 



For consistency the continuity equation must be valid: 

£> + divj = (1.2) 

Proof: 

q — Sq divE 

= e div ( (curl B - /i j) J 

1 

= — div curl B — divj 
div curl B = die ijk djB k = did^B k = (1.3) 

antisym. sym. 

The physical meaning of the continuity equation is charge conservation: 



Q(t) = / g(x,t) d 6 x 
Jv 

Q = / q d 3 x = - / divj d 3 x 



v 



Using Gauss' theorem, 



AwRd'x = fb Rds (1.4) 

one can see that the change of the total charge inside a volume V is equal to the 
total current which flows out of it: 



Q = - fP jds (1.5) 

JJdV 

Here dV denotes the surface (= boundary) of V. 

We now will derive differential equations of second order from Maxwell's 
equations: 

(curl curl E) { = e iifc a,-(curl E) k 

^ijkdj€ k l m OlE m 

Exploiting e ijk e klm = 5uS jm - 6 im Sji, one finds 

(curl curl E) { = 8u8j m djdiE m - 5 im 5 j id j diE m 



1.1 Maxwell's Equations 
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(curl curl E) { = did m E m - d^Ei 
= VidivE - AEi 
curl curl E = VdivE - AE 



Using 



it follows that 



Vdiv£ — AE = -curl B 



1 



-dtHo (j + £o E 



—Vg - AE = -fi d t j - eoHoE 



c 2 



1^1 



A E 



d . 



-v<? 



The d'Alembert operator \j$jfis — A^j is often represented by □. 
Similarly one can show that 



B 



-curl E = —curl (curl B — fj, j) 



1 



—B = - ( V divB -AB) + /x curl j 



(1.6) 



1 a 2 



(1.7) 



These equations are wave equations*. One can distinguish between several special 
cases: 

1. j, Q = => free propagation of waves, the differential equations are linear 
and homogeneous. 

2. j, q are given functions =>• inhomogeneous differential equations describing 
the electromagnetic fields caused by currents and charges. 

3. j, q depend (approximately) linearly on E, B . 

linear dielectric media (polarisable media), described by linear differen- 
tial equations. 

4. j, q depend nonlinearly on E, B. 

=> nonlinear differential equations, generally hard to solve (Nonlinear Op- 
tics). 
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Maxwell's Equations and Gauge Fields 



1.2 Scalar Potential and Vector Potential 

Maxwell's equation (11.11) states that divB = (= diBi). According to Eq. (11.31) 
we generally have div curl A = 0. We therefore can write the magnetic field in 
the formQ 

A = Vector potential 



B = curl A 



B does not possess any sources (in some speculative theories such sources are 
considered and represent magnetic Monopoles). 
Inserting this into (11.11 ) we find 

=^curl E = —B = —curl A 
=^curl (E + A) = 



we have generally curl grad 
^E + A = 



ei jk djd k (f) = 



E 



-V0 - A 



in a simply connected space 
6 = scalar Potential 



The potentials A and are not measurable, in contrast to E and B. The equations 
of motion for the potentials are given by 

Q 



-div grad — div A 
curl B = curl curl A = /i (j + so[— V0 — A] 



divE = — 




(1.8) 



If divA = would hold, <p would obey Poisson's equation and the evolution of 
A would be governed by a wave equation. In this case (!) <p represents the usual 
Coulomb potential. We will see below how to achieve this. 

1.3 Gauge transformations 

B = curl A does not fix A completely. Because curl grad x = holds for any x, 
the vector potential 

A' = A + V% 



'We disregard subtle topological effects which appear in multiply connected spaces. 



1.3 Gauge transformations 
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does lead to the same magnetic field. Likewise, <p and 

<p' = <p - d t x 

create the same electric field since 

E' = -V0' -A' = - V0 + Vx - A - Vx 
= -V0- A = E 

These transformations are called Gauge transformations. The physical fields E, B 
are not changed by a gauge transformation. Gauge transformations can be used 
to impose certain (convenient) constraints, called gauges, on the potentials. The 
most common gauges are 



divA = Coulomb gauge 



divA 



Lorentz gauge 



(1.9) 



Both are independent from each other and generally cannot be fulfilled simultane- 
ously. Coulomb gauge leads to the Coulomb potential but is not covariant under 
Lorentz transformations. Lorentz gauge is covariant and therefore very important 
for relativi stic situations. 

Remarks on the derivation of Coulomb gauge: 
Let A', 4>' be some potentials. We need to find a function x so that divA = after 
a gauge transformation: 

A = A' - Vx => divA = divA - Ax = 
This is a differential equation for x, 

Ax = divA' Poisson equation for x , (1-10) 

with the solution 



X(x,t) = - 



1 f divAV) 



47T 



d 3 x' 



\x — X' 



Proof: 



A x x = --7- f divAVJA*— L- ; 

47T / \X — X' 



d 3 x' 



We now use the very important theorem 



\x — x' 



-Att5{x — x) 
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Maxwell's Equations and Gauge Fields 




divA'(x')(-4ir)6(x - x ') 



divA'(a:) 



qed. 



Inserting the Coulomb gauge Eq. (11.91) into Eq. (11.81) leads to the field equa- 
tions for the potentials in Coulomb gauge, 



Interpretation: = Coulomb potential 

A = electromagnetic waves 
Attention: the physical interpretation of the scalar and vector potential depends on 
the gauge condition! 

1.3.1 Aharonov-Bohm effect 

The Aharonov-Bohm effect M7\ is one of the most intriguing effects in quantum 
mechanics and demonstrates how subtle nature can be. Consider an electron inter- 
ference experiment where the electrons are moving around a solenoid (see figure) 



Inside the solenoid the magnetic field is nonzero, but the electrons are confined 
to move in the region where the magnetic field vanishes. We therefore would ex- 
pect that the solenoid does not generate a phase shift in the electron's interference 
pattern. However, this is not the case. A vector potential that corresponds to a 
vanishing magnetic field away from the origin is given by 




(1.11) 



Q 




A 



AB — 




V x A 



■AB — 



for r ^ 



(1.12) 



1.4 Multipole Expansion 
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with r = y^x 2 + y 2 and = (— sin 0, cos 0, 0) where (j) is the angle of cylindri- 
cal coordinates in the x-y plane. 

The phase shift induced by a vector potential can be calculated using the WKB 
approximation (see, e.g., chapter 9 of Ref. |43 1) and is given by A(p = q/h <f> A ■ 
dx. If we integrate along a circular trajectory we have dx = re^ d<p and hence 

q f 2n 

Aip = - rd<p A AB • 

= (1.13) 

n 

Now it seems that the solenoid induces a phase shift even though the electrons 
always move in a field-free region. To make matters worse we can apparently (!) 
change this result by applying a gauge transformation: the new vector potential 

A' = A + V\ with x = (1.14) 

seems to be zero everywhere, so the phase shift calculated with A' should be zero, 
too. 

The resolution of this paradox is subtle. First of all, the phase shift (11.131) has 
been experimentally verified shortly after its prediction. But why is the phase shift 
predicted with the aid of A 1 wrong? The answer is that A' cannot be set to zero 
everywhere because the gauge function \ is iU defined at the point = = 2ix. 
In a proper treatment, which is mathematically advanced, a change of gauge that 
is defined everywhere will not change the AB phase shift. More on this can be 
found in Ref s. ||44ll45l . 

To conclude this section I also want to mention that the AB effect is an ex- 
ample of a topological effect. The phase shift does obviously not depend on the 
radius of the circular trajectory chosen. In fact, it can be shown that Aip does not 
depend on the shape of the trajectory at all as long as it is enclosing the solenoid. 
Effects that do not depend on smooth changes of the shape of an object but de- 
pend on its overall structure are called topological. Typically topological effects 
appear when the space in which particles can move are multiply connected. In the 
AB effect an electron can take two principal ways around the solenoid (left and 
right) so that the space in which it can move is two-fold connected. This is the 
mathematical reason behind the AB effect. 



1.4 Multipole Expansion 

Multipoles are a very important tool in Electrodynamics. They are, for instance, 
used for the derivation of selection rules of quantum mechanical atoms and molecules. 

Basically the multipole expansion is an expansion in terms of angular momen- 
tum eigenfunctions very similarly to quantum mechanics. 



8 Maxwell's Equations and Gauge Fields 

Scalar Multipoles 

Example: scalar potential in Lorentz gauge. 
Fourier transformation in time leads to 



M x ) = J e~ lM (j){x, t) dt 

=► (A + k")Ux) = , fc = ^ 

e c 

This is the inhomogeneous Helmholtz equation. Green's function for Helmholtz's 
equation fullfills 

(A + k 2 )G(x, x 1 ) = -5{x-x) 
For the boundary condition G — > for \x — x'\ — > oo it is given by 

i p ik\x—x'\ 

G(x - x') 



Air \x — x' 

oo 

' =0 m=-2 basis function 

where 

(1 \ ^ sin 37 
— <9 X spherical Bessel functions (1.15) 

x J x 

( \ f 1 ^\ ^ COS 37 

/i, =ji — i—xf ( — £? x ) spherical Hankel functions (1-16) 

\x J x 

The solution is given by 



x) = G(x, x 



B 



1.4 Multipole Expansion 
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Particularly important for practical applications is the field 4> u (x) at a position x 
outside of a charge distribution. Setting the origin of the coordinate system equal 
to the center-of-mass of the charges one has r = \x\ > \x'\ = r' and therefore 



r , r> = r. 



ik 



oo I 



(*) = -££ Y lm {^)hf\h 

fcn 

u i=0 m=-i 



(1.17) 



In many applications the extension of the charge distribution is small against 
the wave length A = ^ = — (e.g. an atom or molecule: r' — 1 A, A = 1 ^m). 
The argument kr' of ji is therefore small. For x < 1 we have 



x 



3i{x) 



(21+ 1)!! 



k l 



(2Z + 1)! 

The spherical multipole moment of order /m is defined by 



Q, m := / r n Y; m g u dtt'dr' 



Special cases: 



Q 



00 



Air 

Q 



g u (acf) d 3 x' 
Q = total charge of ion/molecule at u 



Qn 
Qio 




8tt 



Qu{x')(x' -iy') d 3 x' 




71 



{Px ~ iPy) 




(1.18) 



(1.19) 



where p = (p x ,P y ,Pz) is the (Cartesian) electric dipole moment: 



10 



Maxwell's Equations and Gauge Fields 



A point dipole creates a potential of the form 

px 



X 



47T5n|a;| 3 



Physical meaning of a point dipole: to charges of equal magnitude but opposite 
sign at a very small distance L. 

Dipole moment = q(L) ■ Ln = p 

Let L go to and q(L) — > oo, i such a way that q(L) ■ L remains finite. 
The electric field of a point dipole at the origin is given by 

E Dip (x) = - , with x = — 

The general multipole expansion of electrostatical fields is 

. 1 



x 



where 



^x^xlvn XiXj i 
r r 3 2 j r$ 



Qij = J g(x)(3xiXj — 5ijX 2 ) d 3 x 



is the (traceless) tensor of the quadrupole moment. 

Eq. dl.171) can also be derived by an expansion in spherical harmonics. They 
are eigenstates of the angular momentum operator and fullfill the orthogonality 
relation 

y; m ,(tf, <p)Y lm ($, <p) sin?? dd dip = 8 w 8 mm , (1.20) 



'0 JO 

and the completeness relation 

oo I 



V V Y? m (tf,tft)Y lm {#,tp) = <% - <p') 5(cosi) - costf') (1.21) 
z — ' z — ' ^ v " 

Because of Eq. (11.211 ) we have for any scalar function /($, </?): 

m<p) = f w,<p')5(n-a) dsn 

= E / Y *mW)Y lm (n)f(n') dsn 

l,m 

= Y ^ V) " Fl ™ with coefficients^ := J Y* m (ti)f{ti) dSl' 



1.4 Multipole Expansion 



11 



In our case 

<Pu(x) = y~V; m (/t9, (p)(j>u,,lm( r ) 
l,m 

with<^ m (r) = [ Y^Q'^r'^'^^dQ' 



This can be inserted into the Helmholtz equation 

(A + k 2 )^(x) = —g u (x) 

^0 

Laplace operator spherical coordinates: 

A = -jrr 2 - -L ; L = x x -zV 

L 2 Y lm = 1(1 + l)Y lm 

i d 2 d 1(1 + 1) l2 \ ~ . . i _ . . 



r 2. y r [j r r z j £q 

Solving this differential equation with suitable boundary conditions (outgoing ra- 
dial wave for r — > oo, i.e. ~ cos(A;r)i sin(Ar)-) leads to expression (11.171) for 



Multipole expansion for vector fields 

It is also possible to make a multipole expansion of divergence-free vector fields 
(e.g., divi? = 0). One essentially replaces the spherical harmonics Yi m by 

Xi m := =LYi m (d, (f) '"Vector spherical harmonics'" 

V Kl + 1) 

X lm fullfills 

J X* m Xii m >dTl = 5w5 mm > 
J X* m ■ (x x X im )dfi = 

A complete set of vector functions for divergence-free fields is then given by 

Fi(kr)X lm and V x (g t (kr)X lm 
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Maxwell's Equations and Gauge Fields 



with Ft(kr) = F^h^ikr) + F^hf\kr) 
9l (kr) = g^h^ikr) + gf ] hf \kr) 



h^\kr) = spherical Hankel functions h\ z> (kr) = h\ L> *(kr) 
The electromagnetic field in free space can then be expanded according to 

B = [ a£ ^' m ) F i( kr ) x im + a m {l, m)V x g t (kr)Xi 



(2), 



,(1)* 



■ lm I 



Lm 



Lm 



-a E (l,m)V x Fi(kr)Xi m + a m (l,m)gi(kr)X lm 



Fields proportional to CL rri tire called spherical TM fields ('"transverse magnetic'"), 
and fields proportional to a E are called spherical TE fields ('"transverse electri- 
cal'"). 

These fields are associated to the corresponding multipole moments if kr 1 
is valid inside the charge distribution One then finds 



Mi 



a E (l,m) ~ Qi m 
1 



a m (l,m) ~ M« m 



(l+l) 



2C + Jl^J 



d J x 



M Zm is called magnetic multipole moment. More on the vector multipole expan- 
sion can be found in Ref. HI, in particular in chapter 16.2 and 16.6. 



1.5 Electrodynamics in Dielectric Media 

Up to now we have only considered electromagnetic fields in free space in the 
presence of given charge or current distributions. However, one often wants to 
describe electromagnetic fields inside of some medium (e.g., glass or a crystal). 

Almost all media consist of atoms and molecules. In such systems the charges 
are in a bound state, i.e., they cannot be substantially displaced by an applied 
external field. However, they are polarisable or even have a permanent dipole 
moment. Atoms, for instance, are polarisable by an external electric field: 




1.5 Electrodynamics in Dielectric Media 
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Molecules (water, for instance) can also have a permanent dipole moment. This is 
often the case for ionic bounds. If one is only interested in describing macroscopic 
fields which are averaged over a volume containing many molecules, the evolution 
of the macroscopic fields E, D, B and H is governed by 

divD = q div.B = 

(1.22) 

curl E = -B curl H = j + D 

macroscopic Maxwell equations 

The influence of the specific medium is contained in the relations 

D = e E + P ; B = fi (H + M) 

Here D is the dielectric displacement and H the magnetic field, P is the po- 
larisation, and M magnetisation of the medium. There are several special cases: 

• In a vacuum we have P = and M = 0, leading to the ordinary Maxwell 
equations. 

• If molecules have a permanent (electric or magnetic) dipole moment but are 
randomly aligned inside the medium, one still finds M = P = because 
of the averaging. If one aligns the molecules M or P ^ is possible. The 
classical example for a model of this behaviour is the Ising model. 




M=0 MtO 



• If the atoms/molecules are polarisable, one finds 

P = xE ; D = eE = e e r E = (e + x)E 

e is the dielectric constant of the medium and % i ts susceptibility. This ef- 
fect can (particularly in crystals) depend on the direction of the fields (bire- 
frigence): 

£,•»•: dielectric tensor. 
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• According to the model of Fedorov optically active media are described b>H 

M = ; D = e(E + (3 curl E) 

The meaning of optical activity is that the polarization of a light beams 
rotates during the propagation through the medium. 

• In nonlinear media one has, for instance, 

D = e E + k\E\ 2 E 

1.6 Derivation of macroscopic Maxwell equations 

The fields E, D, B, and H describe macroscopic fields which result from an 
averaging over '"microscopic"' fields E m { c and -B m ; c . The latter fullfiil the usual 
Maxwell equations 

div£ mic = — divS mic = 

curl E mic = -B mic curl B mic = /j (j mic + e E mic ) 

Charge and current distributions are divided into two parts, 

£?mic £?free £?bound 
J mic J free 3 bound 

The free quantities describe particles which are not bound to atoms/molecules. In 
classical mechanics their charge distribution can be written as 

£freeO) = ^ qi5(x - Xi) 

i 

jfy ee (x) = qiXiS(x — x^ Sum over free particles. 

i 

The bound quantities correspond to the atoms and molecules of the medium. 

Abound (x) = ^ Q n (x) Sum over molecules 

n 

Q n (x) = qiS(x — Xi) Sum over the particles (e _ , nuclei) of the molecule 

i 

JboundO) = ^j n {x) ; j n (x) = ^2 qi Xi5{x - Xi ) 



2 There are several models for optical activity. We use Fedorov's model just as an example. 



1.6 Macroscopic Maxwell equations 
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In most situations experimental measurements cannot resolve electromagnetic 
fields of individual molecules (size ~ 1 A). In optical experiments the spatial res- 
olution is in the order of a wavelength (~ 6000A). One therefore can perform a 
spatial averaging over E mic and B mic without getting a worse description of the 
experiments. 

E(x) = (E mic (x)) = J f(x')E mic (x - x') dV 
B(x) = (B mic (x)) 

In these expressions f(x) is a function which varies slowly on molecular scales, 
but the support of which is small against the wavelength. Because we consider an 
averaging procedure we of course demand J f(x) d 3 x = 1. 

if(x) 



700 nm 



The equations for E and B then become 

d 



J f{x')B^{x-x') dV 

j f( x ')^-B^ c/l {x-x') dV 

/ f(x')divB mic (x-x') dV 
J . ' 



=0 

= (divB mic (x)> = 
and analogously di\E(x) = ^ mic ^ 

cur\E(x) = -B(x) 
curl B(x) = /io(j mic (a;)) + ~^E(x) 
The mean values can be calculated as follows. 

(^mic) = (£>free) + ( Abound) 

= : g(x) + (f?bound) q(x) : macroscopic charge density 

Abound) = ^ i^{x)) (1-23) 
n(molecules) 
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The size of each molecule is much smaller than the averaging area. We therefore 
find for molecule n 

(Q n (x)) = J f(x')g n (x - x) dV 

= J f(x') £ qi 5(x - x 1 - Xi ) dV 

i(e~ , nuclei) 

i(e~ , nuclei) 

If x n is the center of the molecule we have \x{ — x n \ <C support of /. One 
therefore can perform a Taylor expansion of /: 

f{x £Cj) f {x X n (tCj •£«,)) 

« f(x - x n ) - (x { - x n )Vf(xi — x n ) H 

{Qn(x)} « 9i{/(a5 - a?n) - («i - X n )Vf(Xi -x n )-\ } 

i(e~ , nuclei) 

(g n {x)} « /(as - x n ) % -V/(aj - a: n ) ^(cc, -#„) + • 

i(e — , nuclei) 




g n isthe total charge of the nth molecule. If the molecules are not ionized, q n = 0. 
p n is the molecule's dipole moment. 

(Abound) = - P « ' V ^( a; ~ Xn ) 

n (molecules) 

Define the macroscopic polarization P{x) = f(x — x n )p n ; q n = 



(£>boundO)) = -divP(ai) 



Inserting this and Eq. (11.231) into the microscopic Maxwell equations one gets 

divP- = — (g mic ) 

= — q divP 

D = e q E + P =>■ divl? = e 
A similar treatment of 

curlB = /ioOmic) + \e 



1.7 Plane Waves 
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leads to 



(imic) =3+ curl M + P 
M = (^m n S(x - x n )) 



m r 



E—^—(xi x miXi) = \^ —^—Li 
2m; V ' ^ 2m; 

i(e~, nuclei) i(e~ .nuclei) 



Here m n denotes the magnetic dipole moment of the nth molecule, (details can 
be found in H), section 6.7). This concludes the derivation of the macroscopic 
Maxwell equations Eq. (I1.22I ). 



1.7 Plane Waves 

In vacuum (g = j = 0) the following relations hold, 

UE = UB = ; divS = divB = 

Ansatz for solution: 

E = £k£k exp(ikx — iujt) 

fid 2 \ 
UE = £ k e k I — — - A I exp(ikx - iuot) 

where E k G C is the amplitude and 

e fc the polarisation vector. 

ui 2 

=>- DE = Sk^ki — t + & 2 ) exp(ikx — iut) = 
cr 

This is a solution of DE = if u = c\k\, where u k : = c\k\ is the frequency of a 
photonH with momentum hk. Its energy is given by Planck's equation E = fkj k - 
The polarisation vector e k is determined by the condition divl? = diEi = 0: 

divE = £ k (e k )iiki exp(ikx — iukt) 

= i£ k (e k k) exp(ikx — iaj^t) 
i 

= 



3 The notion of a photon as the elementary quantum of the radiation fields is introduced in 
Quantum Electrodynamics (QED). However, the properties of a photon are closely related to those 
of a classical electrodynamical field. We therefore use this notion already here. 
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=>- e k must be perpendicular to the wave vector k\ 

In M 3 there are two linearly independent vectors perpendicular to k =>- there 
are two independent polarization vetors. The specific direction of e£ is arbitrary. 
One possibility is 

'cos if sin $\ / cos ip cos $\ ^ /— sin 9 

k = k J sin^sini!? J ; ejj^ = I sin <p cos $ J ; = -xe^ 1 ' = I cos 92 

cos $ J \ — sin $ J V 

The real vectors and e^, correpsond to linearly polarized light. Generally the 
polarization is a (complex) superposition of these two vectors, 

e k = ae k 1] + (3e {2) ; with \a\ 2 + \(3\ 2 = 1 , a,/3eC 
An example is circularly polarized light: 

e (±) - —( e {1) ±ie {2) ) 
~ \/2 

Different conventions are used in the literature for the names of circular light 
beams. We use that of Ref. [TT1ITT1 : 

left circular or positive helicity, a + light 
right circular or negative helicity, a~ light 
a + - and a~ light plays an important role in selection rules for atomic transitions: 
cr ± : Am = ±1, where m is the magnetic quantum number of the atomic state. 
General superposition vectors e k = ae^ 1 + j3e k 2 ^ are called elliptically polarized. 

In Quantum Electrodynamics the polarization vector is connected to the spin 
of a photon. Photons have total spin 1, but their spin in z-direction can only 
take the values s z = ±1. This corresponds to a + - and a~ light. The selection 
rule Am = ±1 is then simply a consequence of angular momentum conservation 
when the atom absorbs a photon. The value s z = is not realized in nature 
because of the vanishing rest mass of photons (see, e.g., Ref. |0, section 14.5). 

The magnetic field of a plane wave can be derived from the electric field with 
the aid of the Maxwell equation curl E = —B: 

B = -curl {S k e k e {ikx ~ iuJkt) + c.c.} 
= -e l3k d J (e k ) l £ k e^ kx -^+c.c. 
= -e^e^EJk^- 1 ^ + c.c. 
= -i£ k {kxe k ) l e^ kx ~^ t) +c.c. 



)? t I 5?) , 



1,8 Paraxial approximation, focused light beams 
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To solve this we set Bi ~ exp(— iuj k t). We then find 

B = —(kx e k )e {ikx ~ iuJkt) + c.c. 
= — (k x E) 

The general solution in free space is a superposition of plane waves 
E(x, t)= d 3 k 4 CT)e fe ) exp(ikx - nut) + c.c. 

J (7=1,2 

/£(?) 
d 3 k Y -^( k x e k ] ) exp(ikx - iut) + c.c. 



The following special cases are worth to be mentioned: 

1 . The polarization vector is composed of a real part e and a imaginary part 



e", 



E = £ k e k exp(ikx — iu k t) + c.c. 

= £ k (e' + ie") exp(ikx — iuj k t) + c.c. 

where we assume e' 2 + e" 2 = 1, £ 6 K. 

==>- -E = £fce'2 cos(fc:r — icufct) — £ k e"2 sin(fca; — io^i) 

This is a running wave. 

2. Consider a superposition of two counter-propagating waves with £ k = £_ k 
and e k = e k : 

E = E k + E_ k 

— £ k e k exp(ikx — iuj k t) + £_ k e_ k exp(— ikx — iu k t) + c.c. 

==>- E ~ cos (kx) cos (uj k t) 

This is a standing wave. 

1.8 Paraxial approximation, focused light beams 

The paraxial approximation allows to describe the shape of a laser beam or the 
propagation of a laser pulse in a much simpler way. It reduces the second-order 
wave equation to a first-order differential equation along the axis of propagation. 
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As we will see in Chapter it can also be used to calculate the optical properties 
of a medium. 

The idea behind the paraxial approximation is to consider a light field which is 
a running wave with a spatially (and/or temporally) varying envelope £(x, t). The 
variation of the envelope should be small on the scales of one wavelength (and/or 
one optical period of time). We write the field as 

E = eS(x,t)e- iuJt e ikz ; k = - (beaming-direction) (1.24) 

c 

and consider (here) solutions of the wave equation in free space. It follows that 
= DE 



-d 2 - A ) 5(x,t)e lkz - lut e 



j-^ {d 2 t £ - 2iud t £ - u 2 £) - (Aj_£ + d\£ + 2ikd z £ - k 2 £) j 



^ikz—iuot ^ 



\d?£ - 2ik- 

c 2 c 



-d 2 £ - 2ik-d t £ - A ± £ - d 2 z £ - 2ikd z £ } e ikz - tu)t e 



-2ik Qft + d z J - A ± J £ + ^d 2 £ - d 2 £j e ikz ~ iuJt e (1.25) 

with Aj_ = d 2 + d 2 . We assume that the envelope £(x, t) varies slowly with z 
compared to e lkz , and with t compared to e~ lut . This implies that 

kd z £ > d 2 £ , -d t £ > d 2 £ . (1.26) 

We then can neglect the second derivatives with respect to z and t and are led to 
the paraxial wave equation 



2ik I -d t + d z ) £ = -A ± £ . 



(1.27) 



When the light field is monochromatic the envelope is constant in time and we 
arrive at the Paraxial Helmholtz equation 



2ikdX — —A i £ 



This equation has the same structure as the free Schrodinger equation in two di- 
mensions (z=t). 

We know how to solve this equations from quantum mechanics. If the solution 
has a transverse Gaussian profile of width w at z = 0, the envelope of the laser 
beam is given by 

£(z, q) = E msix ^-e k ^ with q(z) = z - iz R and q = a/x 2 + y 2 



1 . 8 Paraxial approximation, focused light beams 2 1 

The quantity z R = kw^ is called Rayleigh length. The intensity of the beam is 
proportional to \E\ 2 : 







1 


zo z 




1 ^ 

,Wn\/2 




w ^"^^-^ 



The Rayleigh length is connected to the degree of focusing (the minimal 
width) w of the beam. Its physical interpretation is that it corresponds to the 
length over which the area of the light beam increases by a factor of two. The 
smaller zr is, the smaller is also the minimal width w , but also the range over 
which the beam is strongly focused. 

Other important solutions ( Laser modes) are the Hermite-Gaussian beams: 

£ ltm = ^(./^^(y/Voje-^e-^ 1 ^) ( i.28) 

with a(z) = u>q(1 + 2iz/ zr) and Hi the Hermite polynomials. 

|£ u l |£2, 3 I 




-3-2-10123 -3-2-10123 

x/w x/w 



Of particular interest are complex superpositions of these modes. For instance, 
the superposition £ L = £ 1>0 + z£ 01 takes the form 



(1.29) 
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The presence of the factor x + iy indicates that light in this mode possesses an 
orbital angular momentum . In the figure below this is indicated by the color 
change which represents the variation of the phase of the field (red corresponds to 
zero, the colors then circle through the rainbow colors). 




Formally it can be shown (see, e.g., Ref. [9|, Sec. 14.3) that the angular momen- 
tum operator for electromagnetic field is given by 



(1.30) 



It is easy to check that £ L has an orbital angular momentum of M 12 = h along 
its propagation axis. Such a field can be generated using specially shaped phase 
plates. An interesting application of these fields is that in the interaction between 
atoms and light not only energy and momentum but also angular momentum is 
conserved. One therefore can transfer orbital angular momentum from light to 
atoms. This effect has been used to generate vortices in Bose-Einstein condensates 

mm. 



Chapter 2 

Quantum Field Theory of Light - 
QED 

2.1 Many-Particle Theory of Quantum Mechanics 

We assume that the reader is familiar with the Schrodinger equation for a single 
particle. If one considers two non-interacting free particles, the Hamiltonian is 
given by the sum of the individual Hamiltonians, 

H= Pi_ + A (2.1) 
2m 2m 

The corresponding two-particle wavefunction has the form 

ip(x 1 ,x 2 ) = (f) n (x 1 ) ■ (f) m (x 2 ) (2.2) 

In this wavefunction the particles are distinguishable (through their quantum num- 
bers n, m). However, experiments demonstrate that identical particles are NOT 
distinguishable. In the theoretical description this is guaranteed if i'(x 1 , x 2 ) leads 
to the same measurement results as ip(x 2 , £Ci). To examine this more closely we 
introduce the exchange operator E, 

Ei/)(xi, x 2 ) = ip(x 2 , aji) (2.3) 

Because of E 2 = II the eigenvalues of the exchange operator are ±1. The corre- 
sponding eigenfunctions fullfin 

Eip B (xi,x 2 ) = +^b(x 1 ,x 2 ) = ifj B (x 2 ,x 1 ) (2.4) 
E^ F (xx,x 2 ) = -iPf(x 1: x 2 ) = il) F {x 2 ,Xx) (2.5) 



All experimental results indicate that all particles in nature are eigenstates of E: 
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• ipB Bosons 

• ipp ^ Fermions 

The well-known Pauli principle for Fermions is a direct consequence of the anti- 
symmetry of lj)F- 

Spin- Statistics Theorem 

• Bosons <^ interger spin 

• Fermions half-integer spin 

The proof of this theorem essentially relies on causality and can be found in 
Ref. Q. 

Many-particle Wavefunction ij){xx, ■ ■ ■ , x n ) should be symmetric or antisym- 
metric under exchange of any pair of particles: 

1p(^Xi, . . . , Xi, . . . , Xj, . . . , X n ) izlft^Xi, . . . , Xj, . . . , Xi, . . . , Xn) (2.6) 

This can be achieved by symmetrization: sum over all permutations and normalize 

l/j B (Xx, . . . , X n ) = — = ^ *P ( X P(1)> • ■ • » X p(n)) ( 2 - 7 ) 

' P 

i/j F (x u ...,x n ) = — L=^(-l) p ^ (x p (i),...,Xp( n) ) (2.8) 

• p 

whereby p characterizes the permutation: (— l) p = 1 for an even number of ex- 
changes, and —1 for an odd number of exchanges. 

As an example we use the product wavefunction of free particles, 

lj){Xx, X 2 , X 3 ) = 1pl{xx)lp m {x 2 )lp n (x 3 ) , (2.9) 

i.e., particle 1 in state 1, particle 2 in state m, and particle 3 in state n. This means 
that the particle are distinguishable. 

For Fermions the corresponding indistinguishable wavefunction is 

l/;f mn (xx,X2,X 3 ) = -^{xl)i{xx)^ m {x 2 )^ n {x 3 ) - ^l(x 2 )^m(xx)^ n (x3) 

+tpl(x 2 )ljj m (x 3 )ljj n (x 1 ) ~ ^(^^m^^n^i) 
+1pl(x 3 )lp m (x 1 )iJ n (x 2 ) - 1pl(Xx)lp m (x 3 )'4) n (x 2 )} 

For a large number of particles this procedure becomes rather tedious. 



2.1 Many-Particle Theory of Quantum Mechanics 
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2.1.1 Fock Space and Number Representation 

In single-particle quantum mechanics we usually do not need to know the full 
wavefunction ij)(x) in position space. To completely fix the state it is sufficient to 
know a complete set of quantum numbers. For the hydrogen atom, for instance, 
these are the quantum numbers n,l and m (ignoring spin). To know the state it is 
sufficient to specify \n, I, m) instead of ip n im(x) = (x\n, I, m). 

Likewise, in many-particle quantum mechanics it is sufficient to know, how 
many particles are in a state: the ket |l f , l. m , l n ) can unambiguously be associated 
with the state "one particle in state 1, one in m, one in n", i.e., with the wavefunc- 
tion ^Ln( X ^ X 2,X 3 ). 

Example for Bosons: 



+^(xi)^(x 3 )^ m (x 2 ) 
+ipi(x 2 )ipi(x 3 )ip m (x 1 )} 
|2/, l m ) 

Superposition of such states are possible and can also be prepared in an exper- 
iment. Example 

|^) =a|2, J l m )+ / 0|3,,O m )+ 7 |l,,2 Tn ) 
However, also superpositions of states with a different number of particles, 

like 

\i/>) = a|l,) + /3|2,)+ 7 |3i,4 m )... , 

are possible. Formally superpositions of states with different particle numbers are 
an element of Fock space, which is defined as follows: 

• Let H {1) be the Hilbert space for 1 -particle wavefunctions 

• = Hilbert space for 2-particle wavefunctions . . . 
Fock space is the direct sum of all N -particle Hilbert spaces: 

f := n {0) © n {1) © n {2) © ■ • • © n (N) (2.10) 

Tv® describes the vacuum and correspondingly includes only a single state, 
denoted by |0) ((0|0) = 1). 
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Fock space is needed to describe processes in which the number of particles 
is not coserved. For instance, pair production of e + and e~ from a 7 photon is 
described byQ 

|1 7 ) ® |0 e+ ) |0 e _) — > |0 7 ) ® |l e+ ) |l e _) 

A second example is the excitation of one atom from the ground state |1 9 ) to the 
excited state |l e ) by absorbing a photon: 

I Iff) ® \h) — ► |le> ® |0 7 ) 



2.1.2 Creation and Anihilation Operators 

These two types of operators allow a simple representation of many-particle wave- 
functions. The creation operator a \ for a particle in state I (also called mode I) and 
the corresponding anihilation operator ai are defined by 

a\\nf\n^..) = V^+l\{n (l) + 1),, n<™> . . .) (2.11) 

a, K (0 ,nL m) . . .) = V^O |(n« - . . .) (2.12) 

These are the definitions for Bosons. 

These operators fullfill the following commutation relations (For simplicity 
we omit the indices). 

\a\ a] \n) = a) a \n) — aaJ \n) 



= a) \/n \n — 1) — ay/n + 1 \n + 1) 
= n |n) — (n + 1) |n) 
= (-l)|n> 

This implies [a, a^] =1. The action of these operators is equivalent to the corre- 
sponding operators for the harmonic oscillator. The general commutation relations 
are 

[a h a} n ] = Si m (2.13) 



u l t u m 



(2.14) 



The physical reason for this close analogy to the harmonic oscillator is that the 
energy levels of free particles in one mode have the same structure: the creation 



'Actually an additional particle is needed to make energy and momentum conservation possi- 
ble. However, here we are only interested in the underlying principle. 



2.1 Many-Particle Theory of Quantum Mechanics 



27 



of n particles each with energy E corresponds to a total energy nE . This corre- 
sponds to an n-fold excitation of a harmonic oscillator of energy E = ftw 
Obviously one finds for Bosons 

aiW W _^ W (iv-i) 
a, |0) = = Number! 

Here the first two equations define a mapping from Fock space into Fock space. 
The last equation defines the vacuum state and essentially states that one cannot 
anihilate a particle if there is none. 

Because of Pauli's principle it is not possible to have two Fermions in the same 
state. Therefore, the Fermion creation operators need to fullfill (a|) 2 = 0. This 
can be achieved through the anticommutation relations 

{a\,a) m } = {a u a m } = {a u a^} = 5i m (2.15) 

where the anticommutator is defined by {A, B} := AB + BA. 



2.1.3 The Field Operator 

How can one calculate expectation values using creation operators? Let us first 
consider an ordinary 1 -particle operator 0(x, p). In the theory of a single particle 
one has 

(O) = [ ^*{x)6i){x)£x 



If we have N particles in this state the expectation value of O is of course given 
by 

N 

(0) = Y,(0) i 



1 -particle 

8=1 



This suggests to define the expectation value of a 1 -particle operator in the frame- 
work of many-particle theory as follows, 



iV 



(0 N ) := ^(dfapi) ) 

i=l 

r f 

= / d ^ Xl ■ ■ - d3x nij*N(Xl,- ■ -XN)dN4>(Xl,- ■ -Xn) (2.16) 

i=l ^ 

Consider for instance the energy and the corresponding equations for the eigen- 
values, 

H(x ,p)ij) a {x) = EMx) (2.17) 
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For a n -particle state \n a , . . . , n a >, ■ ■ ■ , n a n) the energy is given by 

E = E a n a + ... + E a n a , + ... (2.18) 



The particle numbers n a ,n a >, . . . can be calculated with the aid of creation/anihilation 
operators. To do so we introduce the operator 



N a = a\a a 



(2.19) 



This is the Number operator for particles in state a. Its physical meaning is re- 
vealed by the following calculation. 

(. . . , n a , ... \N a \ n a , ...) = (..., n a , . . . \a\a a \ ...,n a ,...) (2.20) 
= (..., n a , . . . |4v^l •■•>("«- l)o, • • •) 

= (..., n a , ... | y/(n a - 1) + lv^l • • • , n fl , . . ■> = n a {. . . \ . . .) = n a 
And therefore we have for E: 

E= d 3 x^l(x)H(x 1 p)'ilj a (x)(...,n a ,... J n a/ \ala a \...,n a ,...,n a/ } (2.21) 



+ . . . 

In this expression the operators a a , do only act on the many parti- 
cle states, not on the wavefunction ip a (x). 

Contrary wise, H(x,p) acts only on ip a {z) Dut not on |. . . , n , • • •). 
The number representation of many-particle states can be considered 
as a kind of bookkeeping procedure to simplify the combinatorics as- 
sociated with the symmetrization/antisymmetrization of many-particle 
states. 

We therefore can rewrite the expression for E as 

E = (..., n a , n a r, ... | J d 3 x tp* a (x)H(x, y)i) a (x)a\a a \ n a , n a >, ...) + .. . 

(2.22) 
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The field operator is defined by ty(x) := J2 a a a^a{x) and allows us to rewrite the 
mean energy as 



E= (...,n a ,n a >,...\ J d 3 x¥(x)H(x,y)V(x)\... 



n a ,n a >, . . .) 



(2.23) 



Proof: 



;. . . I / a *h*\ ...> = <... I / a E«(x)^E*(«W ■ ■ ■> C"4) 

a o 

= / ^ (• • • \aW\ ■ • •) = Ea (• ■ ■ l a « aa l ' ' ■) = E < 



a n a = E 



One generally can express the expectation value of any 1 -particle operator 0(x, p) 
in a Fock state \ip) by 



(0(x,p)) = (ip\ j d 3 x &(x)0(x,p)V(x)\iJ>) (2.25) 
This holds for any state in Fock space, e.g., 

|V) = | la) + |2a') ^ (Vl^*^) = (lal^lla) + (2 a , |^|2 a ,) (2.26) 

Remarks on the field operator: 

• In a sense, ^f(x) is both an operator and a 1 -particle state: 

^!{x) = Y, a ^a{x) (2.27) 

a a : Operator in Fock space 
i}j a (x): 1 -particle state 

One could (and sometimes does) write ip a {x) as (x\ip a ) \. pait i c i e 

- = (x| ^a a |^) Eint = (sl^Etat (2-28) 



Interpretation of ^(x): What is the meaninf of the state ^(x) |0)? In the 
following we will see that it is a very special 1 -particle state. To see this we 
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first consider a perfectly localized state in Ti.^ (actually, perfectly localized 
states are not normalizable, but this doesn't matter for us): 

{0„}= basis of^ (1) 



<f) xo (x) = S(x - x ) = ^ c « ( 

n 

d 3 x (j)* n (x)(f) xo (x) = J d 3 x (f>* n (x)8(x - x ) = <fr* n (x ) 
^ 'Pxoix) = J ^2 ( t ) n( x o)4 ) n(x) (completeness relation) (2.29) 

n 

The Fock-space representation of this state is given by 

^C(xo) |1„) = ^C(x )4 |0) = &(x ) |0 n > (2.30) 

n n 

Thus, ^(xq) creates a particle at the position xq. 
ty(x ) anihilates a particle at the position x . 

We therefore can (formally) introduce a state |l x ) = ^(x)\0) which corre- 
sponds to one particle at position x. This allows us to find the mode function 
of a given state in much the same way as in ordinary quantum mechanics 
when the state is given by a Dirac ket \<f>). The corresponding wavefunction 
is then given by (x\4>). The same state can be represented in Fock space as 
"one particle in mode 0", i.e., 1 1 ^) . To calculate the corresponding wave- 
function one simply forms (1^ j 1^,) = (0|\I/(a;)|l^). In general, |l x ) projects 
onto the one-particle wavefunction of a given state in Fock space. A 
two-particle wavefunction can be calculated like (l x l y \ijj) and so on. 

Commutation relations for Bosons: 



[^(x),^ f ( 



X 



^2i>i{x)4j* n {y) [ai 



l.n 



^2?p n (x)ijj*(y) = S(x - y) (for Bosons) (2.31) 





[*{x),&{y)] = 


= S(x - y) 






[&(x),tf{y)] 


= [*(x),*(j/)]=0 


For Fermions we have analogously 






{a u a\ n } = 8i m z 


= aia m + a, 








= S(x - y) 





{^(x),^(y)} = {^(x),^(y)} = 



(2.32) 
(2.33) 



(2.34) 
(2.35) 
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2.1.4 Second Quantization 

The construction of field operators described above is mathematically equivalent 
to a second quantization procedure applied on the Schrodinger wavefunction. 
Canonical quantization of a classical particle: 

Using the Lagrange function L(x, x) one can define the canonical momentum by 

P = ^ (2-36) 

ox 

Quantization is achieved by re-interpreting x and p as operators and demanding 
the commutation relations 

[x,p]=ih (2.37) 

The same procedure can be done for a field: 

For a Schrodinger particle one can introduce the Lagrange density 

h 2 

C(ip, if)*, ip, ip*) = ihip*ip V^* Vip - ip*Vip (2.38) 

2m 

The Euler-Lagrange equation for this Lagrangean leads to the complex conjugate 
of the Schrodinger equation, 

n <9£ n DC dC 

dip d(diip) dip 

h 2 

-> dtihip* = d 2 — ip* - Vip* (2.40) 
2m 

One can define the canonically conjugate momentum density by 

dC 

7r( x ) = = ihip*(x) (2.41) 

dip(x) 

If one applies again the canonical quantization scheme one finds 

ip(x),ir(x) ^(x),il(x) 
[W(ar), n(j/)] = ihS(x - y) -> [^(x), ^(^J = 5 ( x _ y ) (2 .42) 

This means that, by applying the canonical quantization procedure a second time, 
we end up with the same commutation relations as that derived from the construc- 
tion of the field operator from creation/anihilation operators. 

The method of second quantization is most popular method used to construct 
a many-particle theory of fields. However, it has the disadvantage that the con- 
struction of Fock space and the meaning of the field operator is not as obvious as 
in the previous procedure. 
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2.1.5 Summary of the construction of a many-particle theory 

Recipe to construct a many-particle theory (quantum field theory) using the ex- 
ample of the Schrodinger field: 

(i) Choose stationary solutions (p n (x) of a 1 -particle theory with (f n \(p m ) = 5 nm 
Te scalar product should be time-independent since otherwise the meaning 
of the modes, and therefore of the creation and annihilation operators, would 
be time-dependent. The interpretation of operators in Schrodinger picture 
would then be more difficult. 



ihd t ip = Hip, (ip\ip') = J d 3 x ip*ip' 
d t (ifj\ip'}= I (V^' + ^V)^ = / 1^(Hi/j*)i/j' + iP*^-Hifj')\d 3 x 



{ip* (HiP') - tp* (Hip')} = (2.43) 



where we have used the hermitean conjugate Schrodinger equation and the 
hermiticity of H with respect to the scalar product. 

(ii) Replace the general solution of the 1 -particle equations of motion, 

Hx) = J2c n ip n (x) , (2.44) 

n 

by the field operator 

^(x) = ^2a n ip n (x) with [a n , aJJ = 5 nm (2.45) 

({a n , a] n } = 5 nm for Fermions, respectively). 
^(x) creates a particle at position x. 

2.1.6 Quantization in Coulomb Gauge 

The Lagrangean density C of the electromagnetic field is given by 

£ oc E 2 - c 2 B 2 = (-V0 - A) 2 - c 2 (curl A) 2 (2.46) 

C does not depend on derivatives of the fields E and B. It therefore does only 
lead to Maxwell's equations if it is considered as a function of the potentials. The 
field equations in Coulomb gauge (divA = 0) are 

□A = ^„i-iv0, A0 = -^ (2.47) 

c 2 e 
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The equation for <p is not dynamical! 
In vacuum the corresponding equations read 



A0 = 



(2.48) 



Boundary conditions: 
Solution: 



DA = 



A(x,t) = A (+) (x,t) 



A { ~\x,t) (2.49) 

mode function 

/ d 3 k ^ £akN ak exp (-iuj k t) exp (ikx) + v cx > (2.50) 
. - A(-\x,t) 



A i+) (x,t) 



Here we have introduced 



the expansion coefficients A ka 

the polarization vectors £ a k 

the normalization coefficients N ak 

AH (x, t) denotes the so-called „positive frequency part" 
which only contains anihilation operators. 



'oc exp (-iu k t)) 



We know that e tkx is a complete basis (because of Fourier transformation, and 
leaving aside issues of normalization). Since all possible factors of e lhx already 
appear in A^ + \x,t), it is sufficient to consider only A^ + '(x,t). Since the clas- 
sical field A(x,t) is real, the corresponding quantum field operator A must be 
hermitean. Because of [a, a)} = 1 it is clear that a cannot be hermitean. Thus, 
the full operator A cannot be composed out of anihilation operators alone. This is 
another (and more important) reason to use A^ + '(x, t) instead of A(x, t) for the 
quantization procedure. 

Let now A^ + \x,t) and A'^ + \x,t) be solutions of DA = 0. A conserved 
scalar product is then given by 



A (+) (x,t) A' {+ \x,t) 



A^ 



A^> A' i+) } 
AH 
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The prefactor appears out of dimensional considerations: 

=► [<A|A')] = 1 



Cm 

Proof of time-independence of the scalar product: 



d t (A\A) = -i e jr J d?x {X* A' + A* A' - A* A' - A* A'} = 



use -4 A - AA = 



d t (A\A) = -ij c 2 J d 3 x {(AA*) A' - A* (AA')} 

= -i j c 2 J <?x {A* (AA') - (AA*) A'} 
= 

where we used J (V.V.A*) A' 3 = - J V,A* VA'j = J A)\X,A'j 
Normalization of plane waves N ka £ko exp (ikx — iujf.t) = (pka(x, t) 
(<Pka\<Pk>a>) = 5 aa <5(k - k') 

= -ijr J & z xN h<J N Val e* ka e hla , 

with J d 3 x e i[k - k ' )x = (27r) 3 5(fc - fe') it follows that 

(fkal^k'a') = ^(2ir) 3 £j^£ko^5(k - k') N ka N ka , 2u k 



^(2tt) 3 2^ N 2 J ay 5(k - k') 
n 

5 a y 5(k - k') 



2(2vr)3e a; fc 
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for t = it follows that H for A( 



xr. 



A(x) = [ d 3 k y^A ka J— e ka exp(ikx) + c.c. (2.51) 

Quantization: replace the amplitude by the anihilation operator with 

a ka ,al, a , =5 a yd(k-k') 

a ha creates a photon ( = quantum of light ) with momentum hk and polariza- 
tion vector e ka . 

In contrast to the classical theory of electromagnetism, intensity and energy of a 
light beam appear in quantized portions. 

Photons are Bosons and have spin 1, but they have only two possible values 
for the magnetic quantum number, m — +1, — 1. These values correspond to left- 
and right-circularly polarized plane waves, m = ( =>- longitudinally polarized 
photons e ka || k ) is not possible because of the vanishing mass of photons 
(see, e.g., Sec. 14.5 of Ref. @). 

The operators for the electric and magnetic field can easily be derived from 
that for the vector potential: 



E(x) 



= ~A(x) 

t=o 



= i 

t=0 



B{x) = votA(x) =i [d 3 ky a ka J— -4 (k x e ka )e ikx + H.c(.2.53) 

J a V ^ > oUJk 

The Hamiltonian can be found using the classical expression for the energy 
density: 



H = J d 6 x [E (x) + c z B (x) 

fak ( j . t 



J <i'kY^ hut jo^Ob, + i j (2.54) 



2 We set t = to derive the time-independent vector potential operator in Schrodinger picture. 
One can also consider t ^ to derive the corresponding operator in interaction picture. 
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(<£4> sum over harmonic oscillators) 

The Heisenberg equations of motion for operators (in Heisenberg picture) are 
generally given by 



iW t O{t) = 0(t),H 



(2.55) 



so that we can deduce 



iHEi(x, t) 



Ei(x, t), j d?k ^ fi^k |aL a *ff + \ 

J ^ V 2(27T)3e 



a>)i ~ 4v e ^K'a'h J ^ k ^ [ a la a ^ + \ 



2(27T)3e 

e ik ' x (e k , (J :) i 5{k - k')5 aa ,a ka - e- ik ' x {el a ,)ia[ lal {-8 a(T ,5{k ~ k ')) 



E = J d 3 kJ2 



n 



Let us now calculate (curl B 



„IKX i T ,,-lKX * 

a ka e £ka + a ka e e ka 



(curlB)i = i j ^k^Siji^-ak^k x e ka )ie ikx < 



h 







2(2vr) 3 £ ^ 



h 



(2.56) 



curl B = - J d 3 k ^ 



(k x (fc x e fc(T )a fcCT e ifca: + x (fc x e* Ja^e"^ 3 

(2.57) 



2(27r)% o; fc 



With fe • £fo = we also find 



kx (kx e ka ) = k(k ■ e ka ) - e ka k 2 = -e ka k 2 



(2.58) 
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so that 



curlB = J d?k^ 



_h 

2(27r) 3 £ w fc 



k 2 {a ka e lkx e ka + a\ a e- ik * ! e* klT } (2.59) 



1 • 



curl B = —E = fi e E 



(2.60) 



This is one of the inhomogeneous Maxwell equations for j — (otherwise we 
have curl B = \E = fi {j + EqE}). Analogously one can show that 



ihB = [B, H] = . . 

B = -curl E 
Because k ■ e ka = k ■ (k x e k(T ) = 

divB = 

and 



div£ = 



Q_ 



(2.61) 



(2.62) 



(2.63) 



The field operators fullfill Maxwell's equations in a vacuum 



The equations of motions are easily solved if one considers the individual modes, 

iha ku = [a ka , H] 



/dV£/iM„,*(k-k'Hv 



The solution is 



Q>ko 



(2.64) 



aka(0) 



where a k(T (0) is the Schrodinger operator for the corresponding mode . The solu- 
tion has the same structure as the respective classical wave. 

Remark regarding the factor of\ in the Hamiltonian: 

This term is constant and therefore does not have an effect on the field dynamics. 
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It represents the energy of the vacuum state and diverges mathematically. How- 
ever, it depends on the boundary conditions imposed on the mode functions. 
Example: To mirrors at distance L (see Fig. I2.ll ). The (ideal) mirrors can be 
thought of as representing boundary conditions for the electromagnetic field. As 
a consequence, only discrete wavenumbers with k n = ^f-n, n £ N instead of 
\k\ £ E are possible. Consequently, the sum over the modes in the Hamiltonian 
only runs over discrete wavenumbers and the vacuum energy is changed (though 
still divergent). 




Figure 2.1: Sketch to the Casimir effect 



If one subtracts the vacuum energy in free space from that between the two mir- 
rors one finds a finite energy difference which depends on the distance L. This 
leads to an attractive force, the {Casimir force), between the mirrors. The Casimir 
effect is extensively discussed in Ref. 0. 

2.1.7 Remarks on the Quantization Procedure for other Gauges 

The fact that we have only quantized A but not <p is connected to our choice of 

Coulomb gauge. 

Reason: DA = A0 = 

==>• (f) is not dynamical, the corresponding conjugate momentum vanishes. This 
can be understood with the aid of the Lagrangean: 




\ J [(-V0 - A) 2 - c 2 {rotA) 2 ]d 3 x 

\ j [(V0) 2 + 20divA + A 2 - c 2 (rotA 2 )]rf 3 a; 



2.1 Many-Particle Theory of Quantum Mechanics 



39 



Coulomb gauge: 



divA = — r = 



Lorentz gauge: 

divA + — = O^L = - / [(V0) 2 - 200— + A 2 - c 2 (rotA 2 )]d 3 x 
cr 2 J <r 

If one now considers the action S = J Ldt one can integrate by part the term oc 
and arrives at a new Lagrangean density £ oc 2 . One therefore gets dC'/d<fi ^ 
so that the scalar canonical momentum is nonzero. 

=>- in Lorentz gauge we do not only have to circularly polarized photons but 
also a scalar photon 0) and a longitudinally polarized one ( < S=>- divA ^ 

0). The two new photons are unphysical and have to be removed by imposing 

appropriate constraints: (divA + 0/c 2 ) \tp) = 0) for physical states \ip). 

A unified prescription to quantize constrained systems has been developed 
(among others) by Dirac, but it is rather sophisticated (see, for instance, chapter 1 
ofRef. ED). 
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Chapter 3 



Mode theory and linear optical 
elements 



The invention of the laser has provided us with a source of highly coherent light 
and laid the foundation for many outstanding experiments in quantum optics which 
otherwise would have been impossible to perform. From an idealized theoretical 
point of view, a laser prepares a many-photon state in which all photons are in 
the same spatial (and polarization) mode. Hence, the state can approximately be 
written as 



where is the creation operator for a photon in the spatial mode 0(x). The 
notion of a mode has been introduced in Chapter |2] and will be studied in detail 
in this chapter. For many applications it is sufficient to think of a laser mode as a 
plane wave; the mode function then reduces to that given in Eq. (12.511 ). i.e., 



However, this is not always sufficient since a real laser beam has a finite width 
which increases over time. A better description can be given using the paraxial 
approximation CSec. II .81) which allows to describe TEM laser modes. The proper 
description of photon modes will be the topic of the first part of this chapter. At 
the end of the chapter we will also consider the possibility to work with non- 
orthogonal modes. This includes the Hong-Ou-Mandel dip and Kerr nonlinear 
media. 

Ignoring the shape of a photon mode for the moment, Eq. (13.11) still represents 
a tremendously simplified form of the real state of a laser. In general, a laser emits 
photons in different modes and does not produce a pure state but a mixed state 
(see Chapter |71). 



M = /(ai)|0> 



(3.1) 




(3.2) 
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However, Eq. (13.11) nevertheless describes many experiments remarkably well 
and is the basis of many theoretical investigations in quantum optics. As a con- 
sequence, one often needs only to consider one or two modes which reduces the 
study of a quantized light field to one or two harmonic oscillators. This simple 
model allows for a convenient description of many experiments which do only 
employ so-called linear optical elements: mirrors, beam splitters, phase shifter, 
polarization rotators, quarter-wave and half- wave plates, polarizing beam splitters, 
and, to some extent, so-called nonlinear crystals. The second part of the chapter 
will provide an introduction to this important field. 

3.1 Mode expansion of the quantized radiation field 

The general construction of a quantum field theory presented in Chapter |2] in- 
cluded the introduction of a field operator which is defined as a sum over sta- 
tionary mode functions times their corresponding annihilation operator. For the 
electromagnetic field this field operator corresponded to the positive-frequency 
part A( + )(x) of the operator for the vector potential. 

However, this is not the most general expansion. A better way is to admit 
a general orthonormal set of mode functions A„(x). One example for such a 
mode function might be a laser pulse at some fixed time. Clearly, such a state 
will propagate in time so that the meaning of what the corresponding creation 
operator actually creates will change over time. While this may appear odd in 
Schrodinger picture it is actually quite natural if we formulate the theory in the 
Heisenberg picture (or in the interaction picture if the electromagnetic field is 
coupled to atoms). Following the derivations in Sec. 12. 1.61 the positive-frequency 
part of the vector potential operator in Heisenberg picture can be written as 



As has been shown in Sec. 12. 1.61 this operator is a solution to Maxwell's equations 
since the mode functions are solutions and the annihilation operators are constant. 
As in any quantum field theory, the field operator can also be expanded in another 
set of mode functions, 




(3.3) 



a 




(3.4) 
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The classical fields A n (x, t) are simply a set of solutions of Maxwell's equations 
in vacuum; plane waves are the most prominent example. Furthermore, the re- 
quirement that they are positive-frequency modes restricts these solutions to be of 
the general form 

A n (x,t) = I d 3 A;^A n (k,a)£ kiCT e ik - x e-^ (3.5) 

J cr 

in terms of an expansion in plane waves. Note that this also guarantees that the 
field modes are transverse, i.e., V • A n = 0. 

There is a useful formal way of dealing with positive-frequency fields. Taking 
the time derivative of Eq. (13.51) one finds 



ihd t A n (x,t) = / d 3 khu k ^2 Ai(k,<r) e^e* 

= tub J d 3 k A n (k, a) e k;(T e' k ' x e 



kx e ~iu) k t 



= fiwA n (x : t) (3.6) 

This is a first-order differential equation in time that be used instead of the (second- 
order) wave equation to describe the dynamics of the Maxwell field. In its deriva- 
tion we have introduced the frequency operator lj which simply multiplies any 
Fourier component k with the corresponding frequency uj k = c|k| . In other words, 
it is diagonal in momentum space. Another operator that is diagonal in momen- 
tum space is the operator — zV. It is then easy to show that the frequency operator 
can formally be written as uj = cy/—A. Such a non-polynomial function of a dif- 
ferential operator is called a pseudo-differential operator. The frequency operator 
is Hermitean with respect to the usual Schrodinger scalar product. With its aid we 
can write any positive-frequency solution to the wave equation or to Eq. (13.61 ) in 
the elegant form 

A n (x,t) =e- <firt A n (x,0) . (3.7) 

The commutation relations of the annihilation operators a n are determined by 
the following 

Theorem: For two normalized transverse positive-frequency solutions of the Maxwell 
equations, represented by vector potentials A n and A m , the commutator of the re- 
spective operators is given by 



a ni a m _ — (A n |A m ) , (3.8) 



where the conserved scalar product between two classical solutions of Maxwell's 
equations is given by 

(A n |A m ) = ^ J d 3 x (a; • A m - A; • A m ) . (3.9) 



44 



Mode theory and linear optical elements 



It is worthwhile to note that this commutator is time independent even if the modes 
are not stationary because of the time independence of the scalar product. We 
remark that the scalar product between two positive-frequency modes can be re- 
duced toQ 

(A n |A m ) = j d?x{A* n - {-i)uA m -iCjA* n - A m ) 
^ / d 3 x (At ■ {-i)uA m - iA* ■ uA m) 



h 
H 



d 3 xAl ■ uA r 



Proof of theorem: The field operators for the electromagnetic field can gener- 
ally be expanded in normalized mode functions as H 

A = ^A n (x,t)a n + A£(x,t)aJ, (3.10) 

n 

E = ^E n (x,t)a n + E;(x,t)a+ (3.11) 
B = ^B n (x,t)a n + B;(x,t)4 (3.12) 

n 

with E„ = — A n and B n = V x A n . The operator a n can then be represented in 
the form 

f a n =(A n |A) . (3.13) 



This also shows that the annihilation operator is time independent since both the 
mode function and the field operator are solutions to Maxwell's equations, so that 
their scalar product is time independent. Eq. (13.131) can be shown as follows. 

(A n |A) = ^ ^ a m(A n |A m ) -j- q| n (A n |A m ) 
fit 

= ^ a m 5 mjn + a t m xO (3.14) 



For later use we also remark that it can also be written in the form 
(AJA m ) = / d^E'.^X 



2 In this chapter mode theory is presented for the free electromagnetic field using the electric 
field operator E. However, after the Power-Zienau-Woolley transformation (see Chapter it is 
usually more convenient to deal with the dielectric displacement D instead. This poses no problem 
since for free fields these operators are related by D = EqE so that D = —A/eq. 
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The scalar product (A n | A.* r ) vanishes because the modes are positive-frequency 
modes, i.e., they can be represented as A n (x, t) = exp(— iut) A n (x, 0) It there- 
fore follows that 



h 



d 3 x (A* n ■ iuA* m - iuA* n ■ A* mJ 



(3.15) 



With Eq. d3T3l we find 

[a n ,al] = [(A n |A), (A|A m )] (3.16) 
^ J d 3 xdV[(A;(x)A(x) - A;(x)A(x)), (A(x')A m (x') - A(x')A m (x'))] 

J d 3 xdV[-(A;(x)E(x) - A;(x)A(x)), (A(x')A m (x') + E(x') A m (x'))] 



. 2 

ie 



h 

We can now exploit that the equal-time commutators [A, A] and [E, E] vanish, 
and that the remaining commutator is given by 

[A,(x),^(x')] = ^^(x-xO (3.17) 

with the transverse delta function 

^(x - x ) = J (<% - ^) exp(zk ■ (x - x')) (3.18) 



This results in 



[a n ,al] = ^y"d 3 xdV^(x-x') x (3.19) 
-A- ni j(x)A m j(x ) — A j(x)A mj -(x ) 



All mode functions appearing in this expression are transverse (V • A n = etc.). 
Since for transverse fields the identity 

J dV ^(x - x')A mj (x') = A m ,(x) (3.20) 

holds, this expression is identical to the scalar product between the modes and 
therefor concludes the proof of Eq. (13.81) □ 



46 



Mode theory and linear optical elements 



The result (13.81) allows for a convenient expansion of all quantum field oper- 
ators in terms of mode operators. For instance, the Hamiltonian of the quantized 
radiation field reads 

H = f y"d 3 x: (e\x) + c 2 B\x)^j : 

= | / d 3 x : {i- A n^n - A>t) . (-A m a m - A* m al) 

n,m 

+c 2 (V x A n a n + V x A* n a ] n ) ■ (V x A m a m + V x A* m a m )^j : 
= — ^ ^ j^ndm / (yA n ■ A m + c V x A n ■ V x A T 

n.m 

+ala m J (A* n • A m + c 2 V x A; • V x A m 

+«L.an J (A n ■a; + c 2 vxa„.vxa; 
+<4<i J (k ■a: + c 2 vxa;.vx a,;) } (3.21) 



Now, 



c 2 / V x A ■ V x A' = c 2 / V x V x A ■ A' 



c 2 



J (-AA + V(V- A)) ■ A' 
(u 2 A) • A' 

^A ■ CjA' (3.22) 

where the second-to-last line follows from the definition of the frequency operator 
and the transversality of the field modes, V ■ A = 0. The last line is a consequence 
of the Hermiticity of the frequency operator with respect to the Schrodinger scalar 
product. These expressions allow us to further reduce the Hamiltonian. The coef- 
ficient for a n a m (and for a^a^) vanishes: 

J [k n ■ A m + c 2 V x A n ■ V x A m ) = J {{-iuA n ) ■ (-iuA m ) + c 2 V x A n ■ V x A m ) 

= J (-u)A n ■ uA m + c 2 V x A„ ■ V x Am) 
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In a similar manner, we find for the other coefficients 



a; • A m + c 2 v x a; • v x A m = / a; • A m + wa; ■ uA m 



J (A; • A m + iuA* n ■ A m ) 
J {-ik* n ■ uA m + iuA* n ■ A 
-ik* n ■ uA m + iA* n ■ wA m ) 



-(A n pA m ) (3.24) 



so that 



H = - ^2 (4 a m(A n |^A m ) + a| n a n (A m |^A n )) 

n,m 

= ala m ( A n \ioA m ) (3.25) 



Often monochromatic modes are chosen so that CjA n = co n A n . For orthonor- 
mal modes it then follows that (A m \uA n ) = uJ n S n ,m so that the Hamiltonian is 
reduced to 

H = a\a n hw n . (3.26) 

n 

3.1.1 General construction of single-photon states 

From the above it is clear that the annihilation operator for a specific mode func- 
tion A n can be written as 

Q>n = (A n |A) 

= y I d 3 x ( a: • a - a: • a 

h 
h 

= 2 J± I d 3 xA * . £A« . (3.27) 



J d 3 x (A* • {-iuk {+) + iuA^) - iuA* n ■ (A (+) + A ( -); 
J d 3 x (a; ■ (-iwA« + iwAH) _ zA; • (wAW + wA^- 
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Now, E(+) = -<%A(+) = iwAW, so that 



2eo 



d 3 xA! • E (+) 



(3.28) 



This means that by specifying the properly normalized field mode A n (x, t) we 
can construct the corresponding annihilation operator using the positive-frequency 
part of the electric field operator. 

A single photon in this mode is described by the state 

m = 4|0) • (3.29) 

If this photon interacts with an atom in the ground state \g), the latter will be 
excited and the photon will be absorbed at the position x of the atom. The cor- 
responding interaction Hamiltonian can be written as — d • E^ + ^(x )|e) (g\ + H.c, 
so that the Rabi frequency of a single photon becomes 



Ml = (0|® (e|d-E( + )(x )|e)(s|at|0>® \g) 
= d-(0|EM(x )4|0> 



= d-[E«(x ),a, 
The commutator can be evaluated by 



(3.30) 



[E 



(+), 



x o),4" 



[E! +) (x ),^||d 3 xA n (x). E (-)(x)] 



2e 
—ih 
Q 2e 

—ih 



d 3 xA n ,,(x)[ES +) (x ),E^(x)] 
d 3 xA nj (x)^^(x - x ) 



(-)/ 



i yd 3 x(wA nj (x))^(x-xo) 



(») 



ia>A nji (x ) 
= -9 t A„ )i (x ) 
= E nj j(x ) 

In step (i) we have used the fundamental commutation relation 



[Ej +) (x ),E^(x)] = ^l^(x-x 



hw 



(3.31) 



Step (ii) follows from the fact that the mode function A„ is transverse. For any 
transverse vector field the transverse delta distribution acts like the ordinary delta 
distribution: 

' d 3 xA nij (x)4(x - x ) = A n>l (x ) . 
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We thus have found the result that the electric field of a single photon in mode 
A n is simply given by E n = — A n . While this may appear self-evident it is an 
important check that the mode formalism is working correctly. 



3.1.2 Gaussian modes in paraxial approximation 

In this Section we will derive explicit expressions for normalized mode functions 
in the paraxial approximation that has been introduced in Sec. 11.81 However, 
since we are dealing with the frequency operator we will first rederive the paraxial 
approximation. 

The paraxial approximation relies on the assumption that the field mode is 
almost monochromatic. In other words, its Fourier transform is narrowly peaked 
around a certain wavevector fcx, which we will take to be fc^ = k^e z This is 
essentially the same statement as assuming that the field E n can be written as 

E n = £ n (x, t)e lkLZ ~ iuJLt , lu l = ck L , (3.33) 

where £ n (x,t) is slowly varying both in space and in time so that \d z £ n \ <C 
fci]£ n | and \d t £ n \ "C 0Ji\£ n \- An immediate consequence of this is that 

V-E n = 

= (V • £ n + ik L ■ £ n ) e^-^ 

« ik L ■ £ n e ikLZ ~ iuJLt , (3.34) 

so that we can write £ n (x,t) ~ e£ n (x,t), where e is a constant polarization 
vector orthogonal to k L . 

Being a positive-frequency mode of the free electromagnetic field, E n fulfills 
the exact equation id t E n = CjE n . In Fourier space the right-hand side becomesH 

d 3 k 

l3; 



uE n (x) = u / ——e lkx E n {k) 



J -^c\k\e lk - x E n (k) , (3.35) 



with 



E n (k) = J d 3 x'e- ik x 'E n {x') 

'^e^eUx')^-^ 



ee 



-iuiLt 



ee~ iuJLt £ n (Ak) (3.36) 



J d 3 x'e- iAk - x '£ n (x') 



3 We suppress the time dependence to keep the notation simple. 
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with Ak = k — k,L- Since by assumption E n (k) is narrowly peaked around k^ 
we can expand the frequency that appears in Eq. (13.351) as 



c\k\ = cv k 



2 



c^(k L + Akf 



= cJk 2 L + 2k L Ak z + Ak 2 z + Ak 2 x + Ak 2 . (3.37) 

Since | Afe| < fc^we can expand the square root in lowest order in the components 
of Ak to obtain 



c\k\ ^ck L + cAk z + — {Ak 2 x + Akl) . (3.38) 

2k l 



Thus, 

cbE n (x) = ee-^t J J^ c \ k \e^ x £ n (Ak) 

/d 3 Ak 
— — -c\k L + Ak\e iAk - x £ n (Ak) 

~ ee ik ^-^ J (ck L + cAk z + ^- {Ak 2 x + Afc 2 )) e iAfc -£ n (Afc) 

« ee ifc ^-^* fcA^ - zed, - Jl- (gg + <9 2 ) J (3.39) 

This corresponds to the paraxial approximation. The corresponding paraxial wave 
equation follows by inserting this in the dynamic equation, 

(id t -u)E n = ee^-^^LSn + id^-uEn 

« ee 1 ^-^ L L + id t - ck L + icd z + ^- {d 2 x + d 2 y )jM0) 

so that 



(id t - ^ee^-^Sn « ee* 1 '"^ + c<9,) + ^- + d 2 )^) £ n 



(3.41) 

In Sec. ll.8l we have studied monochromatic solutions to these equations which 
provide an excellent description of a laser beam. However, monochromatic fields 
are not normalizable so that they are not useful to describe the state of a single 
photon. To get normalizable states we need to admit pulses. To find such states 
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we can exploit the close similarity of the paraxial wave equation and the free 
Schrodinger equation. It is then easy to see that a Gaussian mode is given by 



£ 



E 



■nit) 



2 v (t) 



(-c(t-t )+z-z y 

2c 2 T 2 



v(t) 



2 ,c(t 



k, 



. (3.42) 



In this expression r denotes the pulse duration which corresponds to the time in 
ehich the intensity / = 2£ c|£| 2 (or rather the exponential factor that appears in 
/) drops to 1 /e of its maximum value. Likewise, wq denotes the transverse 1 /e 
width of I for t = t and z = z . E max is the maximum electric field in this mode 
and can be calculated by normalizing the mode, 



2c 



it 1 

2e 



2e 

fVjjT, 



d^|E„| 



d 3 x \£„ 



hkr, 



so that 




(3.43) 



(3.44) 



where V = c2tixw\ is the volume of the mode in vacuum. This is remarkably 
similar to the normalization factor of plane-wave modes if the electromagnetic 
field is quantized in a box. 



3.2 Coherent States 

Coherent states are probably the most important states in quantum optics and rep- 
resent a quantum mechanical description of classical fields. 

Consider a single mode a, [a, at] = 1 of the electromagnetic field. There are 
several definitions of coherent states \a): 

• Eigenstate of a: a \a) = a \a), a E C 

• state of minimal uncertainty (excercises) 

• Displacement of the vacuum 
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What is a \a)l Expand |a) in number states: 

oo oo 

l«) = Y] Oi n \n) = V -j= (a ] ) n |0) 

n=0 n=0 

oo oo 

) = a„a |n) = a n v^ |w — 1) 



n=0 n=l 

With n' = n — 1 one finds 



oo 



a | a) = a n i +1 y/n' + 1 |n') = a |a) = a a„/ |n') 

n'=0 n'=0 

It follows that a n+ i = a , whereby a has to be fixed by the normalization 
condition. From 

a 2 

"2 = -^=«0 



one can deduce 



a 3 

«3 = -^«n etc. 



a n = — f=«o (3.45) 



This results in 



oo „ oo 

a , , a 

la 



n=0 V^l" n=0 

Using (n|m) = 5 nm and Eq. (13.451) one is led to 

|2n 



(a|a) 

v m=u n=U 

This leads to 

12 



) = 1 = Kl 2 V -^==- (n| V -^== |m) = |a | 2 V -^-p- = K| 2 exp |a 



2 

We thus can represent the coherent state |a) as 



«o = ex P ( — ■ L ^ — ) (3.46) 



a) = exp f — y'-) £ -j= \n) = exp f~ ^jr) exp(aa f ) |0) 

^ ' n=0 ^" ^ ' 



(3.47) 
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| a) can also be constructed using the displacement operator 



D(a) = exp (acr — a*a 



To do so we can use the Baker-Campbell-Hausdorf equation: 

e (A+B) = e A e B e -±[A,B} e -\[A\A,B}\ 

Inserting A = aa) and B = —a*a we find 

-\{A,B] = - l -{-\a\ 2 )[a\a] = -^e^ 

=>[A, [A,B]] = 
Together with Eq. (13.491) the then can be brought into the form 



(3.48) 



(3.49) 



D(a) 



e aa e~ a a e 



(3.50) 



This allows us to prove the following properties of the displacement operator: 

• D(a) is unitary: 

D\a) = e a * a ~ aat = e - {aa ^ a * a) = D-\a) = D{-a) 

• We have the relation 

\a) = D(a)\0) (3.51) 

It becomes clear from Eq. (13.511) that the vacuum can be considered as the 
special coherent state a = 0. 



It was shown in the excercises that 



E , 

B ~ i[a 



a + a) 
t _ 



Thus, 



E = (a\ E \a) = (a\ a + a' \a) = a + a* = 2Rea 

B = (a\ B \a) = (a\ z(a' — a) \a) = i(a* — a) = 2Ima 



(3.52) 
(3.53) 



(3.54) 
(3.55) 



4 Strictly speaking this holds only for a certain choice of polarization vector and phase. For 
a single mode we have A oc acxp(ikx)e + exp(— ikx)e* and therefore E = A oc 
iu}aexp(ikx)e + (— iuj)a' exp(— ikx)e* as well as B — curl A oc iacxp(ikx)k x e + 
(— i)at cxp(— ikx)k x e*. Choosing, for instance, circularly polarized light (e oc e x + ie y ) 
and x = one finds that the y components of the fields fullfill the above relations. 
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(AE) 2 



(E 2 ) - E 2 

(a\ a 2 + aa^ + a) a + (a)) 2 \a) — (a — a*) 2 

a 2 + (a*) 2 + a*a + (a\ aa* \a) — (a — a*) 2 

a 2 + (a*) 2 + a* a + (a\ a) a + [a, a'] |ot) — (a — a*) 2 

a 2 + (a*) 2 + a*a + (a\ a'a + 1 |a) — (a — a*) 2 

a 2 + (a*) 2 + 2aa* + 1 - a 2 - (a*) 2 - 2aa* 

1 



(3.56) 



=^> AE is independent from a. Analoguous calculation for AB: 



{ABf 



1 



(3.57) 



=^ AEAB = 1 



(3.58) 



The general form of Heisenberg's uncertainty relation for two (arbitrary) her- 
mitean operators is given by AEAB > \([E, B])\/2. From this we can see that 
coherent states have minimum uncertainty. 

One can characterize coherent states by a complex amplitude a E C and by an 
uncertainty 1 . This translates into the following graphical picture 



Figure 3.1: Visualization of a coherent state by a complex amplitude with uncer- 
tainty 1 (small circle). The dashed circle describes the (free) time evolution of the 
amplitude. 

Free time evolution of coherent states: 



ik Im a 




a 



a(t)) = U(t) |a(0)) oc £/e Qat |0) = Ue°^WU\0) 
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U\0) = |0) \a(t)) = Ue aaf U*\0) 

Ue A U ] = U(l + A+^- + ...)tf 

= UU ] + UAlP + hjA{U ] U)AU ] + ... 

= 1 + UAtf + -(UAtf) 2 + ... 

= exp(UAU + ) (3.59) 

Ue aat U j = e aUa ' u ' = exp^e^V) (3.60) 

=> \a(t)) = \ae~ lwt ) (3.61) 

Hence, coherent states remain coherent states and their amplitude moves on a 
circle around the origin of the complex plane. 

• Coherent states are not orthogonal to each other: 

-lal 2 -|/3| 2 *o 

(a\f3) = e - e 2 e Q,/3 (3.62) 

The proof can be done using the expansion in number states or the displace- 
ment operator and Baker-Campbell-Hausdorf equation (Eq. (13.491 ) on page 

|53. 

• Coherent states are overcomplete: 

J d 2 a\a) (a\ = tt! (3.63) 

Proof: Expansion in number states and integration over the phase angle of 
a. 

3.3 Classical and Quantum Mechanical Interference 

The difference between classical and quantum mechanical interference is at first 
sight rather subtle, but can be well understood using coherent states. Set \ip) = 
\a) = D(a) |0);D t (a) = D~\a). 

D(a) can be used to switch to a unitarily equivalent „picture" in which the 
state |^>) = D~ 1 {a) \ip) coincides with the vacuum. The operators then transform 
as 



6 = D\a)OD{a) 



(3.64) 
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\0\j>) = <V|0|V» 



E(x) = D\a)E(x)D[ 



a 



(a k(7 ...e lkx + U.c. 



ftw, 



dk ^ a ka 

a 

To calculate DUa)aD(a) we use the theorem 



2(27r)3 eo 

+ H.c. + iJ77.eJ hx D^aD + H.c. 



\- e a e ikx D\a)aD(a) + H.c. 




k n = B: 



[A,k n ] (3.65) 



A = a*a — aa , B = a 

k — a; ki — [a*a — aa\ a] 

k 2 = h = . . . = o 



a 



D'(ot)aD(a) = a + a 



E(x) 



E aii other modes + (iy/^~-^ a e lkx a — %^f7T.e* a e lkx cv) 
+ (iaJT~.e a e ikx - iJT~.e* a a*e~ ikx ) 



E{x) = E{x) + E* assica \x) 



(3.66) 



The state \a) in absence of external classical fields is thus unitarily 
equivalent to the vacuum in the presence of the field £^ lass - (x) (and 
Bf-Cs)). 

We can exploit this equivalence to understand the difference between the interfer- 
ence of photons and classical Maxwell fields: 

• Classical interference: 

E = E 1 + E 2 , e.g. Ei oc e ikiX 

I oc EE* = l^l 2 + \E 2 \ 2 + 2Re(E 1 E* 2 ) 
2Re( J Bi J B;) oc cos(a:(Jfei - fe 2 )) 
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• Quantum description of classical interference: 

= 1 0:1,0:2) 

= D 1 (a 1 )D 2 (a 2 ) |0> 

D 1 (a 1 ) = e aia i- a *i ai 
D 2 (a 2 ) analogously 

The product of operators D\D 2 creates the sum of classical fields, 

D\d[ED 1 D 2 = D\(E + E 1 )D 2 = E + E l + E 2 

Transformation using D\d\ from to the vacuum: 

/oc {ip\E^E^\ip) = (o\(e h + e* 1 + e* 2 )(e {+) + e 1 + e 2 )\o) 

= \Ei + E 2 \ 2 

classical interference fringes 

In quantum mechanics one can also consider superpositions of coherent states. 
This does not describe classical fields anymore: 

\if>) = |oi) + \a 2 ) = (£>i(ai)+D 2 (a 2 )) |0) 
Interference fringes for this state: 

I oc (i)\E { ~ ] E {+) \i)) 

= (0\(Dl + Dl)E h) E {+ \D 1 + D 2 )\0) 

= (0|(D^ (_) E m Dx + DlE { - ) E {+) D 2 + D ] 2 E { - ) E {+) D l + D ] 2 E^ ] E {+) D 

= {Ex] 2 + \E 2 \ 2 + (0\(D\E { - ) E i+) D 2 + D ] 2 E { - ] E {+) D^ti) 

= {Ex] 2 + \E 2 \ 2 + ( ai \E^E( +) \a 2 ) + (a^E^ E {+) \ ai ) 



=^J oc \E X \ Z + \E 2 \ 2 + E* 1 E 2 {a 1 \a 2 ) + E* 2 E 1 {a 2 \a l ) 

!«i! 2 I"2l 2 

( ai |a 2 ) = e-^e—^e a ^ 
^|( ai |a 2 )| 2 = ... = e -l Ql - Q2 l 2 

==>- For \ct\ —a 2 \ ^> 1 quantum interference is strongly suppressed. For a classical 
field with a large number of photons | ct | 2 ^> 1 it is therefore unimportant. The 
state |oi) + |o 2 ) is also called Schrodinger cat state since it corresponds to a 
superposition of macroscopic states (for |«j| 2 3> 1). 
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3.4 Linear optical elements 

The mode expansion given above clearly indicates that a complete quantum theory 
of light must contain all possible modes of the electromagnetic field. However, 
in certain situations that are quite common in quantum optics experiments it is 
possible to ignore all but a few or even a single mode. The reason behind this 
is that the free wave equation is linear so that superpositions of mode functions 
again correspond to a mode function. Hence, if we do not observe the shape of the 
mode function but only explore the quantum statistical properties of light (e.g., cy 
counting the clicks of a detector) we may assume that the photons always occupy 
one specific mode whose annihilation operator a corresponds to a single harmonic 
oscillator. 

There are of course a number of conditions when this approximation can be 
applied. Obviously the mode functions of all photons involved must almost per- 
fectly overlap. We therefore need a photon source that always generates photons 
in nearly the same mode. The Hong-Ou-Mandel dip (see below) demonstrates 
what can happen if this condition is not fulfilled. 

Furthermore, the single mode approximation can break down if light interacts 
with a medium in such a way that photons are distorted in a different way. This 
is in particular the case in nonlinear media where the presence of other photons 
distorts the mode of a photon. In that case the assumption that all photons occupy 
a given mode clearly breaks down. It may even happen that the notion of a single- 
photon mode becomes meaningless because genuine two-photon states (e.g., more 
complicated entangled states) are needed for a proper description. 

Despite these restrictions the single- (or few-)mode approximation is extremely 
useful if only linear optical elements are used to manipulate photons. Ignoring the 
dynamics of the mode then allows to describe the quantum aspects of photons with 
a very simple harmonic oscillator model that is tremendously helpful to sharpen 
the intuition about quantum mechanics. 

In this section we will discuss a simple description of lossless linear optical 
elements in the context of the few-mode approximation. Linear optical elements 
include phase shifters, polarizing and non-polarizing beam splitters, and polariza- 
tion rotators (e.g., quarter-wave plates). 

3.4.1 The basic principles of lossless linear optical elements 

Lossless passive linear optical elements are characterized by three features: 

• Absence of losses: this is an idealization that works quite well for linear 
optical elements. It states that no photons are lost when light passes through 
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the element. In quantum mechanics this translates into the unitarity of the 
evolution. 

• Passivity: there is no source of energy in these elements that could change 
the energy of the light field. In other words, the energy of the incoming light 
beam is conserved. Furthermore, the usual case in quantum optics is that 
the frequency of the light is not changed either, energy conservation then 
implies that the total number of photons in passive linear optical elements 
is conserved. 

• Linearity: there are no nonlinear effects associated with an optical element. 
This means the evolution of the quantum state does not depend on the pho- 
ton number or the intensity of the light beam. 

Taken together these conditions imply that lossless passive linear optical elements 
may be described by a unitary transformation U that conserves photon number and 
energy. One can phenomenologically construct U by making a suitable ansatz for 
an effective Hamiltonian inside the medium and calculate the corresponding time 
evolution. In the following we will do this for a number of common elements. 

3.4.2 Phase shifter 

A phase shifter can simply be a piece of dielectric material like glass that can 
change the phase of a light beam relative to the phase it aquires during free evo- 
lution. In practise, however, thinks are more complicated |j: one often wants to 
change the phase continuously which cannot be accomplished easily with a piece 
of glass. Instead, a mirror is fixed on a piezo crystal and put somewhere into the 
path of the light beam. Applying a voltage to the crystal then allows to change the 
position of the mirror by very small amounts. The resulting change in the optical 
path's length then generates a phase shift. 

For simplicity, and to motivate the theoretical description for the phase shifter 
given below, we stick to a model in which the phase shift is generated by a dielec- 
tric medium that slightly shifts the frequency of light. The dynamical phase factor 
for light propagating for a time t is given by exp(— iut) where t is the duration 
the light is inside the phase plate©. If the frequency of light inside the phase plate 
is slightly different, say u>', then the effect of the medium would be to generate 

5 I thank Alex Lvovsky for pointing this out to me. 

6 Actually the frequency remains constant but the speed of light is changed so that the time t 
needed to propagate through the dielectric, but for the sake of developing a unified phenomeno- 
logical picture of linear optical elements we deal with the medium as if it it was changing the 
frequency 



60 



Mode theory and linear optical elements 



a relative phase shift of exp(— iAut) with Auj = u' — u. To construct the uni- 
tary transformation Ups corresponding to a phase shifter we therefore assume an 
effective interaction Hamiltonian of the form 

Hp S = hAua ] a . (3.67) 

The unitary transformation U PS = exp(—itH PS /h) then becomes 

Up S = exp(i0a*a) (3.68) 

with (f) = —Aujt. 

Let us explore the effect of this operator on a general quantum state \ip) = 
X^o^™l n ) of the single-mode harmonic oscillator. Because \n) = (a^) n \0) / ' \fn\ 
we can express this state also as some function f(cJ) = fn(^) n applied to the 
vacuum state. Therefore 

oo 

71=0 

oo 

= £/»(EWE&) n tfps|0) • (3.69) 



n=0 



Now ?7 PS |0) = |0), the problem of finding the effect of U P s on a general state is 
therefore reduced to finding its effect on a) . This is given by the 



Phase shifter rule: 



'PS 

Proof: for any two operators A, B one has 



U PS aUl s = e'^a (3.70) 
UpsatUl = e i(t> a ] (3.71) 



°° K 

e A Be~ A = Y;JV> withA '° = B ' Kl+1 = ^ (3 " 72) 



, 11 



(this statement can be proven by calculating the Taylor expansion of e xA Be~ xA 
around x = 0). In our case we have B = a and A = icj)a)a, so K = a and 
K\ = i(j)[a^a, a] = —icpa. Continuing this leads to K[ = (—i(j)) l a from which the 
first rule easily follows. The rule for a) is simply the adjoint of the rule for a. 




Graphical representation of the phase shifter operation U PS 
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We have to check if this definition of the action of a phase shifter is in agree- 
ment with the observations for a classical beam of light. This means that, if the 
light is in a coherent state |a), the mean (complex) amplitude of the electric field 
should exhibit a phase shift of <\>. This is shown as follows, 



Therefore the complex amplitude a of the coherent state is changed to ae 1 ^ as 
expected. We remark that the step in which Up S is moved into the argument of 
the exponential can easily be proven when one writes the exponential as a Taylor 
series. 

3.4.3 Beam splitter 

A beam splitter is a device that splits a classical light beam into two beams with a 
certain ratio of intensities. In quantum optics the analogous action would be that 
an incoming photon in mode a; n should be prepared in a coherent superposition 
of the two outgoing modes a out and 6 out . Because a beam splitter also does not 
change the total number of photons, unitarity dictates that we have to include 
a second input channel b m to describe a beam splitter (see figure). The reason 
is that otherwise a single state |l ain ) would be mapped to some superposition of 
two output states, which cannot be accomplished with a unitary transformation if 
different superpositions are possible. 



U PS \a) 



Ups exp(aa' — a*a)\0) 

C/ps exp((W — a* a)Up S Ups\0) 

exp([/p S (<W — a*a)£/p S )£/ps|0) 

exp(ae^a t -a*ae-^)|0) 



(3.73) 




Graphical representation of the beam splitter operation U BS 
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There are two different conventions in the literature with respect to denoting 
operators for input and output channels. The notation used in the previous para- 
graph assigns different names to input and output channnels. However, one can 
also avoid introducing a name for the output operators and express everything in 
terms of the input operators. This is the convention we will be using here. 

To construct the unitary operator £7 B s we devise a Hamiltonian that describes 
the transfer from a photon in mode a to mode b and vice versa. A phenomenolog- 
ical ansatz that preserves the total number of photons is given by 

H BS = (tfa + (*a j b (3.74) 

Similarly as for the phase shifter we find the corresponding unitary operator to be 

U BS = exp(6tfa - 9*a ] b) . (3.75) 

In complete analogy to the phase shifter we can again infer that £/bs|0) = |0) and 
that the transformation of a general state can be deduced from the transformation 
law of the creation operator. We use again Eq. (13.721) to find this. 

U BS aUl s = e ^a-e*a% ae -e^ a+ e'a% (3y6) 

We then have K = a, K 1 = 6*b, and K 2 = -\6\ 2 a = -\6\ 2 K . It is then easy to 
see that K 2 i = (— \6\ 2 ) l a and K 2 i+i = (— \9\ 2 ) l 6*b. If we sum the terms for even 
and odd / separately we can see that they correspond to the Taylor expansion of 
the cosine and sine function, respectively: 



U BS aU} 



BS 



oo 

E 

1=0 


l\ 




oo 

E 

n=0 


(2n)! 


oo 

^E 

n=0 


oo 

E 

n=0 


(-1)1 


6\ 2n 


(2n 





K> 



2n+l 



(272+ 1)! 

i) n \e\ 2n e* 



E 



|0 

For all operators involved we find the following 



(2n + l)\ 

n=0 V ' 

9* 

a cos 101 + b — sin \6\ (3.77) 



7 If polarization is important these rules should be applied for each polarization direction, see 
the remarks given for the polarizing beam splitter below. 
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Beam splitter rules: 

9* 

UbsuUIs = a cos \9\ + b — sin \9\ (3.78) 

fI 

U BS bU^ s = b cos \9\ -a — sin \9\ (3.79) 

fI 

[/ Bsa tf/t s = a t cos |0| + &t —sin |0| (3.80) 

\9\ 

U BS b ] Ul s = cos \9\ -a f ^- sin \9\ (3.81) 

fI 

These equations describe a general lossless beam splitter. If the input-output no- 
tation is used one simply identifies a out = U BS aUl s anc * so on ' ^ we a PP^ these 
rules to a superposition of two classical light fields, 

|^in) = D a (a)D b (/3) |0) = exp(aa f - a* a) exp(0 - (3*b)\0) (3.82) 

we find, in the same way as for the phase shifter, that the classical complex am- 
plitudes transform as 

9* 

«out = ol cos \9\ — f3— sin \9\ (3.83) 

fI 

Q 

Pout = j3 cos \0\ +a— sin |^| (3.84) 

fI 

which exactly corresponds to the classical behaviour we would expect. 

The modulus and phase of 9 determine the reflection and transmission coeffi- 
cients of the beam splitter and the phase shift between a reflected and transmitted 
wave, repectively. Clearly the transmission coefficient for the electric field am- 
plitude is given by T = cos \ 9\ and the reflection coefficient by R = sin \ 9\. The 
phase shift is not treated consistently in the literature because in a real beam split- 
ter the relative phase between a reflected and a transmitted wave is random (e.g., 
making the propagation length of these two beams inside the beam splitter equal 
to within one wavelength would be extremely costly). One therefore can choose 
the phase of 9 at will. Two conventions are often used: Arg 9 = n/2 and Arg 
9 = 0. We will use the latter. For a % 50 -% 50 beam splitter with \9\ = 7r/4 these 
leads to the simple transformation law 

U BS aUl s = A={a + b) (3.85) 
U BS bUl s = A=(b-a). (3.86) 
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We remark that this transformation law corresponds to a specific physical realiza- 
tion of a beam splitter which, however, is not realized in practise: imagine the 
beam splitter to be composed of two prisms with refractive indices n>, n < (where 
n> > n<) which are glued together. It is well known that light that is reflected 
from an interface between two dielectric media acquires a phase shift of ir if it 
approaches the interface from the side with the lower refractive index. Imagine 
that beam a enters the beam splitter through the n > side. It then does not aquire a 
phase shift. However, beam b then enters the beam splitter through the n < medium 
so that it aquires a phase shift of 7r when it is reflected and leaves the beam splitter 
through the a output mode. This n phase factor just corresponds to the minus sign 
appearing in Eq. (I3.86I) . 



3.4.4 Polarization rotator 

This is a device that allows to generate any superposition of polarization vectors 
in a beam of light. It can be realized by a suitable arrangement of half-wave and 
quarter- wave plates. The latter are pieces of birefringent material that induce a 
phase shift of ix (half-wave) or n/2 (quarter-wave) between two orthogonal lin- 
ear polarization states. Here we will follow a very common convention and take 
horizontal and vertical polarization as basis states. The annihilation operators for 
photons with horizontal or vertical polarization in beam a are then denoted by a# 
and a v , respectively. 

Since a polarization rotator takes a photon with one polarization and trand- 
forms it into a photon with another polarization (conserving the number of pho- 
tons) the phenomenological Hamiltonian for this process should look like 

#pr = (a v a H + Ca ] H a v ■ (3.87) 

While the physics is very different, this Hamiltonian is mathematically equivalent 
to that of a beam splitter. We therefore can immediately write down the action of 
the corresponding unitary operator 

f/p R = exp(x a v a H ~ X*^H a v) , (3.88) 

which leads to the 
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Polarization Rotator rules: 
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Graphical representation of the polarization Rotator operation Upr 

The half-wave plate itself is actually not a polarization rotator but a phase 
shifter for a particular linear polarization. If it shifts V relative to H its action 
is summarized by a|l#) + /5|ly) — ► «|1#) — /3|ly). This can be accomplished 
by Up S y(n) = exp(z7ra| / ay) or a v — > —a v . Likewise, a quarter-wave plate 
corresponds to Upsy(n) = exp(2|a y ay) or a v — > ia v . 

If the half wave plate is rotated by an angle i3 it generates a ir phase shift for 
a photon in the state cos $ a v — sin a# and leaves cos t? a# + sin ay invariant. 
For the mode a# this behaviour corresponds to the transformation law 

an — > a# cos 2d + ay sin 2i9 . (3.93) 

Similarly, application of a rotated quarter-wave plate then allow us to generate an 
arbitrary relative phase between the two modes in this supoerposition. In this way 
Upr can be generated using phase plates. 



3.4.5 Polarizing beam splitter 

A polarizing beam splitter is a device that transmits horizontal but reflects vertical 
polarization. It obviously corresponds to a unitary operator that acts on ay only. 
We can use the beam splitter operation with 9 = n/2 to describe this: 



Polarizing beam splitter rule: 



UpbscihUI qs = a H (3.94) 

^pbsOv^pbs = b v (3.95) 

CWWpbs = 4? (3-96) 

EWtr*4s = h\ (3.97) 
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It should be mentioned here that an ordinary beam splitter ideally acts on both 
polarizations in the same way. Thus, if polarization plays a role in an experiment, 
the beam splitter rules should be applied to the two pairs of differently polarized 
input modes a#, bu an d ay, by in the same way, e.g. a# — > {o-h + / etc. 
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Graphical representation of the polarizing beam splitter operation U PB s 

3.4.6 Schwinger representation of angular momentum algebra 

There is a nice interpretation of linear optical elements in terms of rotations. To 
see how this works we introduce the Schwinger representation of angular momen- 
tum algebra, 

Ji = -(a^b+tfa) (3.98) 

J 2 = % -{b ] a-a ] b) (3.99) 

J 3 = -(0*0-6*6). (3.100) 

It is easy to check that these operators fulfill the angular momentum algebra 
[Ji, Jj] = iEijkJk- As usual we can find common eigenstates of J 2 and J 3 . Some 
algebra yields 

J 2 = -(a*a + 6*6) 2 + -(a*a + b ] b) (3.101) 

Hence states with a total number N of of photons correspond to a representa- 
tion with J 2 = N 2 /4 + N/2. Equating this with the usual angular momentum 
eigenvalue /(/ + 1) of J 2 we find that states with a total number N of of photons 
correspond to a spin / = N/2 (irreducible) representation of angular momentum 
algebra. The magnetic quantum number, i.e., the eigenvalues of J3, are just one 
half of the photon number difference between the two modes. 

Angular momentum theory tells us that an operator of the form U = exp (2i(n- 
J) corresponds to a rotation around the axis n by an angle where n is a unit 
vector. Comparing this to the operators describing linear optical elements above 
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shows that a beam splitter (and a polarization rotator) just correspond to a rotation 
by an angle £ = \9\ around the axis n = sin xe x — cos x e y, where x is the phase 
of 9. Furthermore, if we generate a relative phase shift between two modes a and 
b by applying opposite phase shifters to each mode, U = £/p Sia (0)^ps,b( _ <f>), the 
result corresponds to a rotation around the z-axis. 



3.4.7 Application: Mach-Zehnder interferometer 

As a simple application of the linear optics formalism we discuss a Mach-Zehnder 
experiment as shown in the following figure. 




Mach-Zehnder interferometer 



It consists of a single photon that enters in mode a (red beam). It passes a %50- 
%50 beam splitter and aquires a phase shift if it passes through path b (green 
beam). After passing a second %50-%50 beam splitter we detect whether the 
photon went through output port b of the second beam splitter. 

Mathematically, the input state of the Mach-Zehnder interferometer is given 
by iV'in) = a^|0). After beam splitter one it is changed to 

|^i> = M'r/4)W = -^(o t + 6 t )|0>. 
The application of the phase shifter yields 

hA 2 > = £/p S ,bhM = ^(at + &V*)|0). 

Finally, the second beam splitter transorms the state to 

l^out) = ETbsOt/4) |^ 2 > 

= i ((a* + 6+) +6^-0^)10) 

= (^(l-e^)+^(l + e^|0> 



(3.102) 



(3.103) 



(3.104) 
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The probability to detect the photon in mode b is then given by Pb = \ | 1 + e 1 ^ \ 2 = 
|(1 + cos</>). It obviously displays the phenomenon of interference between the 
two paths of the interferometer because it depends on the relative phase between 
the two paths. 

As a side remark we mention that this calculation can also be applied to dif- 
ferent types of interferometers. It only requires that the photon has the possibility 
to travel along two different paths, but nowhere have we made use of the spe- 
cific form of these paths. We therefore can use it in the same way to describe a 
Michelson, Sagnac, or Fabry-Perot interferometer (see figure below). Of course 
the calculations would differ in a more detailed description which takes into ac- 
count the evolution of the spatial mode of the photon. 
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3.5 The Hong-Ou-Mandel dip and non-orthogonal 
modes 

The Hong-Ou-Mandel dip (HOM dip) is a simple and famous demonstra- 
tion that a naive application of the single harmonic oscillator model is not always 
appropriate. The basic idea of the HOM dip is that the output state of two pho- 
tons that simultaneously enter a beam splitter depends strongly on whether they 
are in the same mode or not. Consider and input state with two photons in dif- 
ferent arms of the beam splitter, but in one case they have the same polarization, 
|VVs) = Oh^hIO)' wne rcas in the second case they have different polarization, 
|?/W>) = Oh^vIO)- The corresponding output states after the (lossless) beam split- 
ter are given by 



l^out.s) = ^BslVVs) 

= ^(4 + ^)(4-4)|o> 

= \{{bif-{aifm (3.105) 

IVWd) = ^bsIVVd) 

= ^(4 + 4)1°) ■ ( 3 - 106 ) 



This means that if the photons are in exactly the same mode they will interfere in 
such a way that they always exit the beam splitter together in the same output arm. 
The probability to find one photon in each of the output arms is zero. On the other 
hand, if the photons are in orthogonal modes they exit the beam splitter randomly 
in the same or in different arms. 

However, in reality the photons will never be exactly in the same mode, and if 
one has less than perfect control about when each photon enters the beam splitter 
their wavepackets might only partially overlap. The input state is therefore a ran- 
dom mix of two photons in modes that are generally non-orthogonal. Therefore, 
the probablility to find one photon in each of the output arms (so-called coin- 
cidence measurements) will generally vary with the delay between the photons, 
which can be controlled by varying the length of the path of one photon before it 
enters the beam splitter. One therefore will see a dip in the coincidence probability 
as a function of the delay. 

The rest of this Section assumes that you are already familiar with the concept 
of a density matrix. If this is not the case you may first want to read the first pages 
of Chapter|7] We assume that the initial state of the photons is given by a mixture 
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of non-orthogonal Gaussian modes, g in = g a ; m <g> £ b ,in with 
Q&M = / dr/ (r)|l ^( T ))(l ^( T )| 



^B,in = J dr7 B (r / )|l B/0(r / ) )(l B ^ (r / ) | , (3.107) 

with J drf a (r) = J dr' fb(r') = 1. The action of the beam splitter on one of the 
states |1 ,# T ), l6,v(r')) results in| 



^Bs|la,0,W = l^i-i) + l^2-o) (3.108) 
1 

- 2 



a{bl, - a\b\) |0> 



+\ (b% - a\a\) |0> 
The probability to detect one photon in each detector is given by 



PuCar^r); 6:^(0) = <*i-il*i-i> 
1 1 

In the derivation we have used that wavepackets traveling in different beams A,B 
have zero overlap. The commutators have been calculated in Sec. 13.11 

We thus have shown that for an arbitrary pair of states in the two beams the 
probability to detect one photon at each output port is given by 

P M (a : 0(r);6 : ^(r')) = \ (l - |<0(t)|V>(t')>| 2 ) , (3.109) 

where the scalar product is given by Eq. (13.91 ) and the modes are normalized 
positive-frequency modes of the electromagnetic field. Of course one has P2-o( a : 
4>(t); b : ip(r')) = 1 — P M (a : 0(r); b : ip(r')) by conservation of probability (can 
also be derived explicitly). The total probability to detect a 1-1 output in state 
is then given by 

P-i = y'drdr'/ a (r)/6(r')Pi.i(a:0(r);6:V(r')) 

= \-\{^*t Ur)h{r') Mr)mr'))\ 2 (3.110) 



8 To keep the notation short we omit here the dependence on t and r' . 
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If the integral on the right-hand side is negligible, we have recovered the result of 
the simple HOM dip explanation presented at the beginning of this section. 

To show that this is the case for realistic pulses we consider the special case of 
Gaussian wavepackets of the form 



e-We^le, (3.111) 



2cke cT^ 

Here T represents the duration of the pulse, in Alex's experiment about 170 fs. 
These wavepackets are only different in their polarization vectors and e^, 
which determine the qubit state of the photons and at the same time determine 
the scalar product for perfect mode matching, 

(0|^)o=(0(r)|^(r))=e;-e^. (3.113) 

It is not hard to derive the mode overlap for different times, 

K0(r)|^(r'))| 2 = £^ Le "^K<^>o| 2 (3-114) 

with u = ck. If we assume that the mode mismatch probability is given by a 
Gaussian of width At, 

_ r 2 

f a (r)=f b (r) = ^- (3.115) 



we find 



1 1 |2 Ar 2 + 4T 2 (T 2 + Ar> 2 

"1-1 = n _ n # o 3 (3.116) 

4T(T 2 + Ar 2 )!^ 2 



2 2' 



In the limit that T and Ar are much larger than 1/oj this becomes 

11 T 
Pi-i = |U|WA | 2 — (3.117) 

11 2 2 l\VIV/0| - V ) 

This means we are allowed to use the simple HOM analysis of the previous section 
whenever the width Ar of the mode mismatch distribution is much larger than the 
pulse duration. If the mode mismatch distribution is narrow, we recover the ideal 
case (i.e., P M = for identical polarization states). 

We can relate this probability to the HOM dip as follows. The HOM dip is 
a decrease from a constant coincidence (1-1) rate for non-overlapping states to a 
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smaller (zero, if the states are completely overlapping) coincidence rate if they are 
pretty much overlapping. If we denote the ratio between these rates as d, then we 
have 



P\-\ | At— »oo 

= l-|(</#)o| 2 Q (3.118) 

with Q = T/\/T 2 + At 2 . This allows us to infer Q = 1 — d where d corresponds 
to the HOM dip for identical polarization states. Assuming that d does not depend 
on the polarization we therefore find 

Pui = \-\\{mo\ 2 {l-d) (3.119) 



3.6 Aspects of nonlinear optical elements 

Nonlinear optical elements, such as a nonlinear crystal like BBO, differ substan- 
tially from linear elements and are generally much harder to study theoretically. 
The polarization P(x) of a nonlinear medium depends nonlinearly on the electric 
field of a light beam and can, to some degree, be expanded in a Taylor series of E. 
Alternatively, one may assume that the nonlinear dielectric medium generates an 
effective interaction between photons. A simple example of a physical process in 
which such an interaction is generated is the case of a three-level atom in cascade 
configuration (see figure; for an introduction into few-level models for atoms see 
chapter^. In absence of light the atoms are in the lowest energy level \g). A 
photon os frequency u\ can pump the atoms into the first excited state \e\). Only 
if this process took place can a photon of frequency 0J2 pump the atoms from |ei) 
to 1 12). Obviously the optical properties for photon 2 depend on the presence of 
photon 1, i.e., an effective interaction between these photons is generated. 



CO, 



Ie 2 > 



CO. 



CO, 



|g> 



A simple cascade model for atoms in a nonlinear medium. 
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This interaction may result in different effects. For instance, the atom may 
decay through the same channel through which it has been excited, thus creating 
photons of the same frequency which it had previously absorbed. Denoting the 
annihilation operator for a photon of frequency Ui by a, we then can expect to 
have an effective interaction Hamiltonian of the form H{ nt = Ka\ (r 2 ttitt\, where 
k is a coupling constant that determines the strength of the effective interaction. 
Of particular interest is the case when toi = u 2 so that these to photon modes are 
equal, ai = a 2 = a. We then find a Kerr-type nonlinearity with 

H Ken = n(a j ) 2 a 2 (3.120) 

which is a very common model in quantum optics. Another important process 
up- and down- conversion. If the atom decays through mode H a 3 two photons 
are absorbed and one other is created. The corresponding effective interaction 
Hamiltonian would look like H{ nt = na\a 2 ax + K*a^a\a\. The first part describes 
up-conversion in which two photons of lower energy are converted into one photon 
of higher energy. Energy conservation implies that u 3 = lui + to 2 must hold. The 
second part describes down-conversion whereby two lower-energy photons are 
created from one photon of higher energy. For oj\ = cu 2 one finds 

-f^conv = na\a 2 + n*a 3 (a^) 2 (3.121) 

which is a very important model for the creation of entangled photons and for 
squeezing. 

3.6.1 Squeezed states of light 

Let us consider the situation in the figure below: a nonlinear crystal whose effect 
can be approximated by H cow is subject to an intense classical field of frequency 
u 3 . 



a a, 
=t > 3 



Basic principle of a squeezing experiment 

If the change in the input beam can be neglected we can replace the operator 03 
by the coherent state amplitude a 3 of the classical field so that we get an effective 
Hamiltonian for mode a alone, 

H sq = ga 2 + g*^) 2 , (3.122) 

9 Selection rules exclude this possibility for atoms in free space. However, the process is possi- 
ble in a crystal or if external fields are applied. 
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In the same way as for linear optical elements we can deduce that the unitary 
squeezing operator, which is usually denoted as S{(3), is given by 

S(f3) = exp(/3(a t ) 2 - 8 * a 2 ) , (3.123) 

Again we can use Eq. (13.721) to derive the transformation properties for the anni- 
hilation and creation operator. As in the case of the beam splitter the series can be 
separated into even and odd I and the corresponding sums can be identified with 
the Taylor expansion of cosh and sinh. This leads to the 

Squeezing rule: 

SaS* = a cosh(2|/3|) - a*^- sinh(2|/3|) (3.124) 

IP I 

Sa j S j = a f cosh (2 1/3 1) - a^- smh(2\f3\) (3.125) 

IP I 

This transformation law, in which annihilation and creation operators are mixed, 
is very famous and also called a Bogoliubov transformation. It appears frequently 
in theoretical physics when a system is governed by a nonlinear Hamiltonian and 
one linearizes the dynamical equations for small perturbations around a given 
mean field solution. An example would be the description of quasi-particles in 
a Bose-Einstein condensate. In our case, the mean field solution is given by the 
intense pump field and the small perturbations are given by photons in mode a. 
Another, more exotic application of Bogoliubov transformations which does not 
rely on linearization is particle creation in a curved space time. For instance, 
Hawking radiation can also be described using a Bogoliubov transformation. 

Before we explore the features of squeezed states of light we remark that 
squeezing is not restricted to a single mode. In general the pump field may gen- 
erate two photons in different modes a and b. The interaction Hamiltonian is then 
replaced by H sq2 = gab + g*cfi$. It is not hard to see that this leads to the 

Two-mode squeezing rule: 



U sq2 aU} q2 


= acosh(2|/3|) - 


-6t£ 

101 


sinh(2|/3|) 


(3.126) 




= a t cosh(2|/3|) 


B* 

~ b W\ 


sinh(2|/3|) 


(3.127) 




= 6cosh(2|/3|) - 


t p 
~ a ¥\ 


sinh(2|/5|) 


(3.128) 




= 6 t cosh(2|/3|) 


(3* 


sinh(2|/3|) 


(3.129) 
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Let us now explore the physical content of the single-mode squeezing trans- 
formation (I3.125I) . We start with calculating the mean number of particles in a 
squeezed vacuum state |0 sq ) = S((3)\0). Using the fact that S^((3) = S(—(3) one 
finds 

n = (0 sq |a t a|0 sq ) 
= (0\S j a^aS\0) 
= (0\(&<fS)(S*aS)\0) 

= (0| (a t cosh(2|/5|) + a^sinh(2|/3|) ) f a cosh(2|/3|) + a) ^- sinh(2|/3|) ) |0) 



= (0|a^sinh(2|/3|)atA sinh(2 | / 5|)|0> 

= sinh 2 (2|/5|) (3.130) 

Hence a squeezed vacuum has a non-zero mean number of particles. The physical 
reason is simple: down conversion in the nonlinear crystal populates the initially 
empty mode a by creating pairs of photons. The energy for this process is pro- 
vided by the pump beam. This is an example of an active optical element. An 
especially interesting feature of a squeezed vacuum becomes obvious in the num- 
ber representation, 

|0 sq ) ~ |0) + 2|/3||2) + -.- , (3.131) 

i.e., a squeezed vacuum contains only even number states. In particular, a weakly 
squeezed vacuum state can be considered as a superposition of a two-photon state 
with the vacuum which can be used to create entangled two-photon states. 

However, the great interest in squeezed states is not related to photon produc- 
tion but to the fact that the mean deviation for certain observables can be sup- 
pressed below their vacuum value. To shed light on this consider the case that 
mode a is initially in a coherent state \a). Down-conversion then transforms this 
state into a squeezed coherent state |a sq ) = U sq \a) = f/ sq -D(a)|0). To experimen- 
tally analyze its properties one generally use the quadrature operators 

1 i 
X{9) = —j=(a + a^) cos6 H — -=(a) — a)sin# 

V 2 v 2 

= -= (ae~ ie + aV) (3.132) 

i 1 
Y{9) = —j={a) — a)cos6 ^(a + a^)sin^ 

V 2 v 2 

= -5= (aV e - ae~ ie ) (3.133) 
v2 
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For 9 = these operators coincide with the position and momentum operator of a 
harmonic oscillator. For each choice of 9 they form a pair of conjugate observables 
because [X(9), Y{9)] = i. We have seen in our discussion of coherent states that 
for 9 = they are proportional to the electric and magnetic field operator in 
specific situations. In general situations X and Y are proportional to E and B for 
some value of 9, but they also test other features of the quantum state of light as 
we will now see. 

The mean value of X and Y for a squeezed coherent state are given by 

(X) = (0\D^S j XSD\0) 

= (0\D ] ^={ae~' w + a ] e i0 )D\O) 
V2 

= (0|^((a + a)e-' ie + (at + a *)e ie ) |0> 

= 4= (ae- ie + a*e ie ) (3.134) 
V2 v ' 

(Y) = -j= (a*e ie - ae~ ie ) (3.135) 

The calculation of the variances follows the same principles but is tedious and left 
to the reader as an exercise. The result is 

AX 2 = - (ch 2 + sh 2 + 2ch sh cos(2# - ip)) (3.136) 

AY 2 = - (ch 2 + sh 2 - 2ch sh cos(2fl - <p)) (3.137) 

with P/\P\ = exp(iip), ch = cosh(2|/3|) and sh = sinh(2|/?|). For the special 
choice quadrature operators that corresponds to 9 = ip/2 we find 

AX = -^-=e m (3.138) 
v2 

AY = ^=e~ m (3.139) 
v2 

Hence the noise in Y is reduced as compared to the shot noise (i.e., the noise of Y 
in the vacuum state) while the noise of X is increased. Note that for this case we 
have AX AY = 1/2 which indicates that squeezed states are states of minimum 
uncertainty, i.e., they fulfill Heisenberg's uncertainty relation for X and Y with 
an equality sign. 

This feature is very useful for a variety of reasons. To see this we represent 
a squeezed coherent state graphically in a similar way as the ordinary coherent 
state: 
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Graphical representation of squeezed coherent states. 



Examples a,b,c represent three different squeezed coherent states of light. Each 
state is represented by an ellipse whose center is at < X > +i < Y > with respect 
to the X, Y coordinate system. The corresponding unsqueezed coherent state is 
represented by a (grey) circle whose area indicates the uncertainty of a quadrature 
measurement. A squeezed coherent state. The ellipse of the squeezed coherent 
states indicates that the noise for X is increased and for Y is reduced. Example 
a) corresponds to a general squeezed coherent state. If the ellipse is arranged as 
in example b) we have reduced noise when measuring the phase of the complex 
amplitude and increased noise when measuring in modulus. This is of practical 
importance in interferometry where improved knowledge of the phase (at the cost 
of increased fluctuations in the photon number) is desirable. Example c) would 
allow to have improved knowledge of the number of photons in mode a at the cost 
of increasing the noise in the phase. In a sense one can consider a photon number 
state \n) as the limit of infinite squeezing of the photon number. 



Homodyne detection 



The question remains how to measure the quadrature operators X{6) and Y{6). 
The standard technique to accomplish this is homodyne detection. 
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l¥> 




b |a) 



Homodyne detection to measure quadrature operators. 



We want to measure the mean value of the quadrature operators for a quantum 
state of a signal light beam in mode a. To do so we use a %50-%50 beam 
splitter to combine this light field with a classical light beam in mode b that is 
called local oscillator and is described by a coherent state \a). This field is much 
more intense than the signal beam, i.e., the mean number \a\ 2 of photons in mode b 
is much larger than (tpla^ a\tp) . The resulting input state can be written as [tj,> in ) = 
Dt,(a)f (g^)|0) where we have represented the state |^) of mode a in the form 
= f(a^)\0) with / some Taylor-expandable function of the creation operator. 
In practise, / is a polynomial of low order. D b (a) denotes the shift operator for 
mode b. 

We measure the photo current at each output channel of the beam splitter and 
subtract the results. Because the photo current is approximately proportional to 
the mean photon number measured by the detector the output signal AI is pro- 
portional to 



Hence, if the phase 6 of the local oscillator (a = \a\e l9 ) is fixed, A J is propor- 
tional to (X(9)). 



AI 



(^n\Ui s (a^a-tfb)U BS \ip in ) 

<^in|| ((^ - fo f )(a - b) - (&t + a t)(6 + a)) |</v> 

-(^|(at& + &V))|^in) 
-(0\pDt(a*b + tfa))D b f\0) 
-(0\f(a\b + a) + (tf + a*)a))f\Q) 
— (i/j\(cra + a*a)\ip) 



(3.140) 
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3.6.2 Nonlinear media of Kerr type 

Together with parametric up- and down-conversion, the Kerr effect is one of the 
most important nonlinear effects in quantum optics. In classical electrodynam- 
ics the Kerr effect can simply be described as an intensity-dependend index of 
refraction, 



where I is the intensity of the field and n 2 the non-linear refractive index. Typ- 
ically n 2 is very small so that large intensities are needed to create a significant 
effect; in (fused silica) optical fibres, for instance, one has n 2 ~ 3 x 1CT 16 cm 2 /W. 
However, in media exhibiting electromagnetically induced transparency n 2 can be 
as large as 0.1 cm 2 /W ll2"Tl . although in this case there are practical restrictions on 
the intensity of the beam. Because a plane electromagnetic wave of frequency uj 
varies in a dielectric as exp(inkz — ojt), where k = lo/c is the wavenumber in 
vacuum, we see that a variation of the intensity / changes the phase of the wave 
as exp(in 2 Ikz). If this effect happens with a single light beam it is called self 
phase modulation. In some media the intensity J 2 of a second beam (with differ- 
ent polarization and/or frequency) changes the phase of the first beam. In this case 
the refractive index n\ of the first beam varies as ri\ = n + n 2 I 2 and one speaks 
of cross phase modulation. 

To construct a quantum- version of the Kerr effect we start by considering the 
intensity operator. The classical intensity of a beam of light can be derived from 
the Poynting vector 



This vector describes the change of energy of the electromagnetic field inside 
some volume. The energy density is given by 



n = n + n 2 I , 



(3.141) 



S = —E x B 



(3.142) 



u = 




(3.143) 
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so that, with V an arbitrary volume in free space and dV its boundary, 

dE f 3 du 



<Fx — 
dt J v dt 



d 3 x ( —B ■ B + e E ■ E 
— I d 3 x (-B ■ V x E + E • V x B 



Ak) Jv 

-— [ d 3 iV-(£x B) 
Ho Jv 

[ d 3 x VS 
Jv 



d z A-S (3.144) 

dV 

where d 2 A is the surface element on the boundary dV which is a vector whose 
direction is normal to the surface. In the derivation we have used Gauss' law and 
Maxwell's equations in free space. 

An operator for the field intensity can be found by replacing the fields in the 
Poynting vector by the respective operators. However, there are some subtleties 
that have to be adressed. The complete operator would look like 

S = ±(E {+) x B {+) + E { ~ ] x B i+) + E {+) x B^ + E H x B^) . (3.145) 

The first and the last term are very rapidly oscillating (as exp(±2itut)) and there- 
fore would average out over the time scale on which a detector reacts. They are 
therefore dropped. The third term would give a contribution even when the state is 
the vacuum, which is undesirable. For this reason, and as it is often done in quan- 
tum field theories to avoid non-zero vacuum expectation values, one exchanges 
the sequence of the operators so that all annihilation operators are to the right of 
all creation operators. This is called normal ordering. We therefore find for the 
Poynting operator the expression^ 

S=i(£ H xB (+, -B H x£ (+) ). (3.146) 

I most cases in quantum optics one is interested in the behaviour of a well 
focused beam or pulse of light. This was the fundamental assumption behind the 
paraxial approximation, for instance. If all photons in the system are composed 



10 The minus sign stems from rewriting the normally ordered form in term of the vector product 
of B and E. 
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of superpositions of plane waves which are close to a central wavevector k = kn, 
where n is a unit vector, then we can use the plane-wave relation B = (k x 
E)/uj k = n x E/c along with E ■ k « 0. This reduces the Poynting operator to 

S « 2ce -B (_) ■ E {+) n . (3.147) 

This suggests to define the intensity operator, which describes the power of the 
pulse flowing through a unit area perpendicular to the radiation direction, as 

I(x) = 2ce E { ~\x) ■ E i+ \x) . (3.148) 

Let us now develop an intuition of what we would expect for a single mode 
theory of a Kerr medium. When classical light propagates along, say, the z-axis, it 
generally aquires a phase factor of the form exp(iknz). The nonlinear phase shift 
is therefore given by exp(ikn 2 Iz), where z is the length of the nonlinear medium. 
Ignoring the difference between group and phase velocity for this estimation, the 
time the light travels through the medium is given by t = z/d = zn /c, where 
d is the speed of light in the medium. Hence, the nonlinear phase shift can be 
written as exp(ikn 2 Ict/n ). 

Turning to quantum mechanics we replace the intensity by the corresponding 
operator. In single-mode approximation the electric field operator becomes E m 
E a (x) a, where E a (x) is the mode function associated with a. The intensity 
operator therefore becomes 

I(x) w 2ce \E a (x)\ 2 a^a . (3.149) 

Since a classical field amplitude acquires a phase shift exp(ikn 2 Ict/n ), the cor- 
responding quantum mechanical observable should be modified by the operator 
exp(ikn 2 Ict/n ). From our discussion of coherent states we know that this ob- 
servable is the annihilation operator a. We therefore expect a phase change of the 
quantized radiation that is proportional to a^a. The Kerr Hamiltonian 

#spm = y (a^a 2 (3.150) 

accomplishes this for self phase modulation (SPM). For cross phase modulation 
(XPM) one would have 

Hxpm = hKa'b'ba . (3.151) 

For the case of self phase modulation, let us consider the Heisenberg equation of 
motion for a, 

ihd t a = [a, H[ nt ] 

= kcJcl 2 , (3.152) 
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which has the solution 

a(t) = e- iKta ° ao a , (3.153) 

with a = a(0) the annihilation operator in Schrodinger picture. Such a Hamil- 
tonian can for instance be achieved by a intensity-dependend energy level shift in 
resonant atoms. We will see how such level shifts are generated when we discuss 
adiabatic elimination (see Sec. 16.4.21 . 

We can now compare this solution to the classically derived expression for 
self phase modulation. In Sec. 13.4.21 we have seen that the phase shift acquired 
by the annihilation operator is just the negative of the phase shift acquired by the 
complex amplitude of a coherent state. We therefore can make the identification 

K = 2kc 2 e —\E a (x)\ 2 , (3.154) 
n 

which obviously depends on the mode. For the Gaussian mode (13.421 ) the maxi- 
mum amplitude is given by Eq. (13.441) . This leads to a nonlinear coefficient of the 
form 

K G =U -7-7=, (3.155) 

HQ TAy/lX 

where A = tvWq is the transverse area of the mode and r the duration of the pulse. 
This is easy to interprete. Pq = Jtlu/t corresponds to the power needed to emit 
a single photon during time r. I = tiw/ (tA) is the corresponding intensity of 
a single-photon pulse if it is confined to a transverse area A. For a r =130 fs 
pulse of a single telecom wavelength photon (A = 1.3/xm) in an optical fibre with 
a core diameter of 2iv = 5/im this amounts to P = 1/uW, J = 6 W/cm 2 , and 
kq = 10~ 15 u;. In other words, to achieve a phase shift of n G t = 1 rad, a single 
photon would need to travel for z = ct ~ 2 • 10 5 km through an optical fibre, 
which clearly is impossible because of losses. To achieve a 1 rad phase shift over 
a distance of 1 km the pulse peak power would need to be P = nP = 0.25 W, 
corresponding to a pulse with 2 ■ 10 5 photons in it. 

Eq. (13.1531) allows us to derive the time evolution of a general quantum state 
of light under the influence of SPM, but it easier to directly calculate the time 
evolution using the Hamiltonian. With an expansion in the number state basis, 
IVK^)) = Yln^nitfln), we get the Schrodinger equation 

idt^n = 2 n ( n ~ !) V>n , (3.156) 

so that 

\if,(t)) = MQ)e~ iKtn{n ~ 1)/2 \n) . (3.157) 
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In particular, a coherent state evolves into 

\^{t)) = e-H 2 / 2 V ^-e-^-^n) , (3.158) 

which is not a coherent state anymore (because the coefficients are not of the form 
a n (t) I y/n).). Hence, in SPM a classical state evolves into a non-classical state. 

Application: controlled phase gate 

An important application of XPM would be the optical implementation of a con- 
trolled phase-gate for quantum information Like the controlled not-gate, 
the c-phase gate together with all possible single-qubit operations would enable 
universal quantum computation. For the logical qubit states |0), |1) and the corre- 
sponding two-qubit state 1 00) , • • • , the c-phase gate has the following truth table, 

, i,j = 0,1. (3.159) 

That is, the logical state |11) acquires a n phase shift and all other basis states are 
unchanged. 

In principle it is easy to achieve this using XPM. Let us consider two-photon 
states where each photon is in a different spatial mode a, b and can have horizon- 
tal or vertical polarization. The encoding of a qubit is given by the polarization 
of each photon, say <-> H and 1 <-* V. For instance, the logical state |01) then 
corresponds to the physical state (v H by | vac) , where we denote the vacuum state by 
| vac) to avoid confusion with the logical states. We now assume we have a non- 
linear crystal that generates XPM only if the photons have vertical polarization. 
This is described by the Hamiltonian 

Hxpm = a v avb v bv ■ (3.160) 

Then the initial state 

IVv) = a o a ^l/|vac) + a ia^fe| / |vac) +ai a|/^|vac) + aiia| / fe| / |vac) (3.161) 
evolves into 

|VWt) = aoo4Al vac ) + a m a} H b v \\&c) + a w a v b j H \vac) + e~ lKt a n a v b v \\<ic} . 

(3.162) 

If the photons propagate through the crystal for a time t = tv/k we have accom- 
plished an optical c-phase gate. 

However, there are several problems with this scheme. The most important 
one is that for nonlinear crystals the constant k is so small that nt <^ 1 in all 
practical cases. In recent years the possibility to achieve large XPM nonlinearities 
in atomic gases (see, e.g., Ref. [50]) has generated new hope in constructing an 
optical c-phase gate, but as we will see in the next two subsections even this 
approach may face some severe problems. 
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3.6.3 Shortcomings of the single-mode approximation I: Mode 
distortion 

While it is appealingly simple for a nonlinear theory, the single-mode approxi- 
mation for Kerr nonlinearities is only of restricted value. Like in all nonlinear 
processes, the assumption that all photons occupy a single mode is often inappro- 
priate, even if this mode function is adapted to the nonlinear problem. 

Let us evaluate this problem in a simple one-dimensional multi-mode theory, 
where we replace the positive-frequency part E^ + \x) of the vector operator of 
the electric field by a scalar operator £(z) that fulfills the commutation relation 

PC*) . ^CO] = - z ') = VS(Z - z>) , (3.163) 

ZEqP 

This model, which is often called scalar QED, neglects the transverse disperion 
and the polarization degree-of-freedom of a light beam, but is otherwise similar to 
ordinary QED. The area F can be interpreted as the area of the light pulse trans- 
verse to the z-axis along which it propagates. For instance, it could be identified 
with the area of the core of a single-mode optical fibre. £(z) creates a (scalar) 
photon at position z. To create a photon in a particular mode, e.g., a Gaussian 
wave packet 4>{z), we can intruduce the corresponding annihilation operator as 

a, = - dz (f)*(z)S(z) , (3.164) 

which ensures that [S(z), at] = <p(z) so that we can extract the mode function 
from the commutator. 

We assume that the nonlinear Hamiltonian is of Kerr-type at each point inside 
the medium. We therefore can write is as 

H = | J dz' {£\z')) 2 £ 2 (z') . (3.165) 

Let us consider the dynamics of an initial state which is a superposition of a one- 
photon state and a two-photon state, 

|^(0)) = ±al\0) + ±(a\) 2 \0) (3.166) 

Because the Kerr Hamiltonian does not change the number of photons we can 
write the state at any time as a superposition of a general one- and two-photon 
state, i.e., 

\il>(t)) = J dz P(z,0)£^z)\0) + J dzdz' Q(z,z',0)£ j (z)£\z')\0) 
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with the time dependent functions P(z, t) and Q(z, z', t) fulfilling the initial con- 
ditions 

P(z,0) = (3-167) 

Q{z,z',0) = ^4>(z)^z') (3.168) 

The Schrodinger equation id t \tp) = H\tp) then yields the differential equations 

id t P{z,t) = (3.169) 
id t Q(z,z',t) = kt] 2 S(z- z')Q(z,z,t) (3.170) 

This means that the single-photon part is unperturbed, P(z, t) = P(z, 0). How- 
ever, the two-photon part will evolve into a state that cannot be written as a product 
of two mode functions anymore, Q(z, z', t) ^ <fii(z, t) 4>2(z', t). This means that 
the two photons become entangled, it then wrong to assume that there is a mode 
function that is occupied by both photons. 

While it would be instructive to go in detail through this calculation, it is 
somewhat plagued by the appearance of the 5 distribution in Eq. (13. 170b . This is 
a consequence of our over-simplified Hamiltonian where we tacitly assumed that 
we have a point-like interaction, 

H = | J dzdz' £\z')£\z)5{z - z')£{z) £{z') , (3.171) 

rather than having an interaction potential of finite range, 

# = I J dzdz' £\z')£\z)V{z - z')£(z) £{z') . (3.172) 

A treatment based on an extended interaction potential V(z — z') is possible but 
beyond the scope of this course. 



3.6.4 Shortcomings of the single-mode approximation II: clas- 
sical SPM 

Potential problems of the single-mode approximation can also be seen in the evo- 
lution of a classical state. Let us check if the Kerr Hamiltonian does correctly 
reproduce self -phase modulation of a classical beam of light, i.e., let us calculate 
the evolution of the mean electric field amplitude if the initial state is a coherent 
state. One easily finds in Heisenberg picture 

E{t) = (a\a(t)\a) 

= (a\e- lKta o ao a \a) 

= a(a\e- iKta ° ao \a) . (3.173) 
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This can be evaluated using 
Proposition 1: 

(a\e lKt ^ a °\a) = exp (|a| 2 (e ifti - 1)) . (3.174) 

Proof: Let I(t) = (a\e lKta o a °\a) . Clearly 7(0) = 1. The time derivative of I is 
given by 

d t I = iK(a\ala e iKta ° a °\a) 

= ina* (a\a e iKta ° ao \a) . (3.175) 

To move a to the right we note that a e tKta o ao = e i^<4>^(e~ iKta l ao aoe iKta ^ ao ). 
The term inside the parentheses is equivalent to the free time evolution of a har- 
monic oscillator of frequency — k. It is well known that the latter is given by 
e tKt a . We therefore find 

d t I = ikq* (a\e iKta ° a °e iftt a \a) 

= iK\a\ 2 e iKt I(t) . (3.176) 

This is a differential equation for which has the solution (13.1741) □ 

As a consequence, the mean electric field amplitude of a coherent state evolves 

as 

E(t) = aexp (\a\ 2 (e iKt - 1)) (3.177) 

This result has some intriguing consequences. Let us first consider the case of a 
usual nonlinear crystal in which n 2 , and therefore k, is small. In this case we can 
assume that nt 1 so that 

E(t) ^ ae iKtH2 . (3.178) 

This is just the usual result for SPM. Note that the smallness of nt does not imply 
that the total phase shift A0 SPM = Im(|a| 2 (e lK * — 1)) ~ Kt|a| 2 needs to be 
small; for sufficiently large intensity \a\ 2 of the coherent signal field we can easily 
achieve phase shifts in the order of 1 rad or more. 

On the other hand, if we consider a giant nonlinearity for a weak signal field 
the situation becomes very different. In this case we have A0 SPM = | or | 2 sin(/ci) 
which does not grow linearly in time. For a small value of |a| the SPM effect is 
small regardless of the size of nt; neither a giant nonlinearity nor increasing the 
length of the nonlinear medium will result in an increased SPM phase shift. Even 
worse, the phase shift is unavoidably connected with a loss of fidelity because 

| exp(H 2 (e iK ' - 1))| = exp(|a| 2 (cos(fi;f) - 1)) < 1 . (3.179) 
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This equation essentially states that in the presence of a large nonlinearity a co- 
herent state will evolve into a very different state rather than just picking up an 
SPM phase shift. 

This is seemingly in disagreement with experiments and has led to a number 
of theoretical investigations (until about the mid nineties) who are trying to adress 
this problem using a multi-mode description, i.e., a description that includes the 
dynamics of the quantum field operator. However, there seems to be no general 
agreement yet on the best way to describe Kerr media in QED. 

As one example for a model, I will briefly discuss the approach taken by 
Boivin, Kartner, and Haus iBH . The central ingredient of their theoretical model 
is to include a time-dependent nonlinear response rather than one which is con- 
stant in time. The corresponding Heisenberg equation of motion for S(z, t) takes 
the form 

d z S = inn NL (z, t) £(z, t) (3.180) 

with 

/oo 
dt' f(t') E ] (z, t - t') S(z, t-t') + m(z, t) , (3.181) 
-oo 

where /(£') describes a delayed Kerr- type nonlinearity and the operator m incor- 
porates the noise associated with the nonlinear medium. It can be shown that 
this model predicts the usual SPM phase shift for a coherent state, regardless of 
the strength of the interaction. However, the associated noise inevitably leads to 
a loss of coherence which would make Kerr nonlinearities unsuitable to construct 
a controlled phase gate for quantum information ll52l . So far this theoretical pre- 
diction has not been experimentally verified. Furthermore, the lack of agreement 
about which model properly describes Kerr media would call for applying alter- 
native models to study this problem. 



The reason why it is necessary to include a noise operator can be motivated by considering a 
damped harmonic oscillator with Heisenberg equation of motion dta(t) = —70(f). The solution 
a(f) = exp(— 7t)a(0) would not fulfill the equal-time commutation relations [a(t),a^ (t)} = 1 
anymore. To cope with this one introduces a noise operator 771(f) which commutes with a, <r 
and fulfills [771(f), rrv(t')] = 2j. The modified dynamical equation, d t a(t) = —70(f) + m, 
has the solution a(f) = exp(— 7f)a(0) + L dt ' cxp(— 7(f — t'))m{t'). It is then easy to check 
that [o(f),o t (f)] = 1. While this may look like a mathematical trick, the noise operators have 
a concrete physical meaning: in any lossy medium the losses are generated by the interaction 
between the atoms of the medium and the radiation field. The noise operators are an effective 
description of the noise in the radiation field induced by the atoms. 
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Chapter 4 



Quantum Mechanics of Atoms: a 
short Review 



4.1 The Hydrogen Atom, Parity 

Hamilton operator for a proton and an electron (V= Coulomb potential): 

H =S; + £: +v(lx '- x ' l) (41) 

R = M 1 +m (MpX p + m e x e ) is the center-of-mass coordinate, r = x e — x p the 

relative coordinate between the two particles, M = M p + m e , - = ^- + The 
Hamiltonian then can be rewritten as 

P 2 p 2 

H = + — + V(\r\) (4.2) 

In traditional atomic and molecular physics the relative motion is the central 
topic, in atom optics the center-of-mass motion. 

Introduce spherical coordinates for the relative motion: 



P 2 K 2 n ^ 



2/i 2/x \ r r 



,2 



(4.3) 



with p = —ihV and X being the orbital angular momentum. One finds [L, \r\] 
0. The eigenstates of H are therefore also eigenstates of L: 



i> = Y lm (d,cf>)fi(r) 
L 2 Y lm = h 2 l(l + l)Y lm 



(4.4) 
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The radial part fi(r) of the wavefunction fullfills 



fc 2 i , h 2 1(1 + 1) 

- — d r rd r + 

z/ir 2/j, r z 



E fi = \ - — - d rrdr + — + V ( r ) ) fi ( 4 - 5 ) 



Solution: f nl ~ r'exp ^— • Laguerre polynomials, a = Bohr radius 
« 0,5 A. 

^nim = fniYi m (9, (p) \nlm) (electronic state) (4.6) 

Here n denotes the principal quantum number, / the orbital abgular momen- 
tum, and m the magnetic quantum number. There is one additional quantum num- 
ber: Spin s of the electron, s = ±|, = \nlms) . 

Parity: 

r — ► —r 

\r\ — ► |r| (4.7) 

Y lm (e, <p) — ► y /m (7r - e, <p + tt) = (-1)^40, 0) 

=>• The parity of \nlms) is (—1)'. Parity is very important for the dipole selection 
rules of atomic radiation. 



4.2 Fine structure and spin 

Of high practical relevance is the interaction of an atom with a static magnetic 
field, 

H = -fiB (4.8) 

where n is the magnetic dipole moment of the atom. A classical particle with an- 
gular momentum L (e.g., on a circular orbit) and charge q possesses the magnetic 
moment 

V = 7TT7 L - ( 4 -9) 

The orbital contribution of the electron is 

li e = -^—l. (4.10) 



2m 



e 



For the proton we have because of M p 3> m 



p 
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Spin of the electron: Q 

2m e 2m e ^ ^ 

Fine structure is a consequence of the magnetic field which is created in the 
rest frame of the electron because of its motion around the nucleus. In the elec- 
trons rest frame the nucleus represents a moving charged particle, corresponding 
to an electric current that induces a magnetic field. 

Rest frame of the nucleus: electric field E = — VVcw> B = (the nucleus 
creates no magnetic field as it is not moving in this frame). 

Lorentz transformation to the electron's rest frame H: 



B' = v x E 

= -V X Wcoul 
= -JjP X VVcoul 



(4.13) 



where v is the velocity. 

Vcoui = V Cou i{\r\) Wcw = e r d r V cou [ = r-d r V Co ui 

i/i \ i n r \ (4 - 14) 

=>& = -- {^drVcoul) P x r = - (- r d r V Cou i) 

Thus, there exists a coupling between spin and orbital angular momentum of the 
electron: 

H S B = —^eC B 

t( 1 w ^ (4 - 15) 



> r 

An exact calculation of the matrix elements {n'l'm! s'\Hs Bairns), including 
relativistic corrections, leads to 

AE FS = --m eC 2 (aZ) 4 ^- ( — ^—r - — ] 

2 1 > W + i An J (4.16) 

AE FS /h w lOOGFz 



'The factor of 2, which appears here in contrast to the orbital angular momentum, is called 
gyromagnetic factor g. It follws from Dirac's equation for a relativistic electron that g = 2. 
However, vacuum fluctuations of the electromagnetic field cause small deviations from this value 
(ffw 2,002). 

2 Actually the electron is moving on a non-inertial (circular) orbit and it is necessary to proceed 
differently. However, the result presented here agrees with the correct one up to a factor of 1/2 
(Thomas factor) and explains the relevant physical ideas. 
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where j — I + ~ denotes the total angular momentum of the electron. 

Without relativistic corrections AEps would depend on I and s instead on 

3=1 + 8. 

4.3 Hyperfine structure 

Hyperfine structure is caused by two different effects: 

1 . Interaction between the electron's angular momentum and the nuclear spin. 

2. QED corrections (vacuum fluctuations) 

ad 1.: A magnetic dipole moment xreates a magnetic field which other magnetic 
momenta can interact with. For the hydrogen atom the total angular momentum 
j — s + I of the electron creates a magnetic field at the position of the nucleus 
with which the latter' s spin interacts. 
Hyperfine spin: 

F=S k + j (4.17) 

Is has to be remarked that the addition of Sk and j to F is not fully equivalent 
to the addition of the total spin Sq = Sk + S e and the electronic orbital angular 
momentum L to F. This happens although the operators themselves are exactly 
equal. The reason is that with each of the two additions one associates a certain 
set of basis vectors diagonalizing F. The two bases associated with F = Sk + j 
and F = Sg + L are different, see Appendix ll0.2l 



Chapter 5 



The Interaction between Atoms and 
Light 

5.1 Minimal Coupling 

The classical Maxwell equations are gauge invariant, 

A^A' = A + V X , <P^<P' = <P-d tX . (5.1) 

This invariance holds also in the presence of charge densities g and currents j and 
should also be present in the quantum description of electromagnetic fields. 

5.1.1 Derivation of Minimal Coupling 

The equation of motion for charged particles, 

mx = qE + qx x B electric + Lorentz force 
can be derived from the Lagrangean , 

L = -mx 2 + qxA(x, t) — q<f)(x, t) , (5.2) 

in the following way. The canonically conjugated momentum p is defined by 

dL 

Pi:= — =>Pi = m±i + qAi(x, t) (5.3) 

Therefore, in the presence of an electromahnetic field the canonically conjugated 
momentum is different from the kinetic momentum II = mx. The Euler-Lagrange 
equations are given by 

d dL 
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Insert <E21> and (IP) : 

9 9 
—Pi = mxi + <7t^A {x{t),t) + ?7^rA 

=^ mi 8 = - qd t Ai + qx k (diA k - d k Ai) (5.5) 

verwenden —di<p — d t Ai = Ei (by definition) 

=>■ mxj = + qx k (diA k - d k Ai) (5.6) 
However, Bj = (curl A)j = e jmn d m A n 

'* ^ikjBj £ikj€jmndmA n 
^jik^jmndm-^-n 
i&im&kn ^in^km)d m A n 

= d t A k - d k Ai (5.7) 
mxj = qEi + qx k e ikj Bj 

= qEi + g(i x B)j qed. 

Using the Lagrange function one can derive the Hamiltonian and Schrodinger's 
equation: The classical Hamiltonian is defined by 

H = piii - L 

Using x — — (p — qA) we want to write it as a function of p and x =>- 

1 11 1 

H = Pi — {pi - qAi) - -m—(p - qA) 2 - q—{p - qA)A + qcf) 
m 2 m z m 

H = ^{p-qA) 2 + q<t> 

The corresponding quantum mechanical Hamiltonian can be found by canonical 
quantisation: 

x -> x ; p -> p ; [Xi,pj] = 6ij 
In position representation: 

x — »• x ; p — > —iHV 

and therefore 



1 o 

H= — (-ihV - qA(x,t)) + q<j)(x,t) , ihd t ip(x,t) = Hip(x,t) 
2m 



(5.8) 

This is the Schrodinger equation for a non-relativistic particle in an electromag- 
netic field. The form (p — qA) is called minimal coupling. 
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5.1.2 Gauge Invariance of Minimal Coupling 

Applying the transformation (15.11 ) to the Hamiltonian (15.81) one could be led to the 
conclusion that the Schrodinger equation is not gauge invariant: 



H' 



1 2 d\ 

— (-ihV-qA(x,t)-qVx) + q<f)(x,t) - q—if) 



ihd t ip(x,t) 7^ Hif)(x,t) 

To restore the invariance under gauge transformations one also has to unitarily 
transform the wavefunction: 



This implies 

ihd t ip' 



if>' — if; exp(iqx/fr) 



dt 



J Hif) - q^i) 
at 



(5.9) 



— (-iW -qA- qV X f + 
2m 



(e*«/ty) - q^' 



H'if)' 



(5.10) 



The minimal coupling scheme is of outstanding importance for high energy 
physics. Apart from gravity, all fundamental forces (that is, electromagnetic, 
weak, and strong interaction) do have the same basic structure. The simplest 
form of a gauge invariant theory is realized for the electromagnetic field. It is also 
called U(l) gauge theory because the transformation of the state (15.91 ) corresponds 
to a multiplication with a Unitary 1x1 matrix (= complex number of modulus 1). 
The other fundamental forces couple more fields. For example, the Special unitary 
SU(2) gauge theory of (electro-)weak interaction transforms the fields according 
to 

i ) = cxp(,,T ■ x) ( t 



(5.11) 



Here ip e is the wavefunction of the electron, if) Ue that of the electron-neutrino, and 
<j are the Pauli matrices. The word special implies that the matrix exp(icr ■ x) h as 
determinant 1. The gauge potential for this theory can be written as A = A^'ai 
and therefore conists of Three vector fields A^ . These three fields do roughly 
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correspond to the W± and the Z Boson which mediate the weak interaction be- 
tween different particles, in the same way as photons do for the electromagnetic 
force. Q 



5.2 The Power-Zienau-Woolley Transformation 

Although minimal coupling is of highest importance for the description of the 
fundamental forces, it has several disadvantages when one wants to describe the 
interaction between atoms and light: 

• For electrically neutral atoms which are composed of electrons, neutrons, 
and protons, minimal coupling is somewhat inconvenient. It would be better 
to describe the interaction using the atomic multipole moments. 

• In minimal coupling the gauge potentials instead of the electric and mag- 
netic field describe the coupling. Although this is not a principal problem 
a direct coupling to the physical fields would be manifestly gauge invariant 
and more intuitive. 

• In Coulomb gauge the scalar potential (ft can be identified with the unre- 
tarded Coulomb potential, since the latter solves Poisson's equation (which 
the equation for 0). The fact that the Coulomb potential is unretarded makes 
it appear as if minimal coupling would produce an instantaneous long-range 
force between the atoms. 

The Power-Zienau-Woolley transformation (Refs. lESI and lE9l . see also QUI 
and [31]) is a unitary transformation of the minimal coupling. It leads to an 
equaivalent description of the interaction, but with the above problems removed. 



5.2.1 The Basic Idea 

To derive the so-called dipole coupling Hone can start from the fact that the Euler- 
Lagrange equations are not changed when one adds a total time derivative to the 
Lagrnangean, 

, Tl , d „, . ddL' dL' ddL dL 
L (*,*,*) = L(z,z,t) + -F( X ,t) =►__-_ = __-_ 



'Actually the weak and the electromagnetic forces are entangled. Both together are described 
by a U(l)xSU(2) gauge theory. U(l) thereby refers to hypercharge and not to the electromagnetic 
charge. 

2 A more precise name which is also used in the literature is multipolar coupling since the 
full Power-Zienau-Woolley transformation leads to a coupling of the electromagnetic field to all 
multipole moments of the atoms. However, in the applications one very often only takes the dipole 
moment into account. This is also the case for the simple example in this subsection. 
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For simplicity we consider here only constant fields E and B (for the general 
proof see below). 

E = E ; B = B <p = -xE 

1 

A = --(x x B ) 

Choosing F = —qxA one arrives at 

jF = -qxA+^(xx x)B (5.12) 
L' = -mx 2 + qxA + qxE —qxA + -ix x x)B (5.13) 

V v ' 

L 

1 q 

= -mx 2 + qxE H IB mit I = (a; x mi) (5.14) 

2 2m 

Setting d := qx electric dipole moment of the particle one finds for the new 

fi :— y~L magnetic moment of the particle 
Lagrangean 

L' = \^mx 2 + dE + /iB 

Advantage: Now the canonical momentum agrees with the kinetic momentum, 
p = mx. In addition, the particle now couples to the physical fields and therefore 
is manifestly gauge invariant. It also involves the dipole momenta which are better 
suited for electrically neutral particles. 



5.2.2 The Complete Transformation 

We consider an ensemble of point particles with charges q a and positions x a . 
Their charge density is given by 

g{x) = s ^q a 5{x - x a ) (5.15) 

a 

Since we are particularly interested in the behaviour of a gas of neutral atoms, we 
assume that the total charge is zero, J2 a q a = 0. The polarization of macroscopic 
Electrodynamics fullfills divP = — (f? ge b) (see Sec. ??). For the present case its 
solution is given by 

P(cc) = V/ duq a (x a - R)S(x- R-u(x a - R)) (5.16) 
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Proof: 



diPi(x) = / du q a (x a — R)idiS(x — R — u(x a — R)) 

= ^ / duq a (x a - R)i5'(xi- Ri-u(x a - R)i) x 

S(xj - Rj -u(x a - R)j)5(x k - R k - u(x a - H>)0) 

where we sum over i and k) = (1, 2, 3), (2, 3, 1), (3, 1, 2). 5' is the derivative 
of the 5-distribution. We have 

{x a - R)iS'(xi - Ri- u(x a - R)i) = -d u 5{xi R, u(x a - R)i) (5.18) 

and hence 

d,iPi{x) = — / du q a d u S(x — R — u(x a — R)) 
= -^2q a {S(x - x a ) -8(x-R)} 

a 

= Q5(x -R)~Y^ qJ(x - x a ) (5.19) 

Q 

The total charge Q vanishes for neutral atoms, and the second term corresponds 
to the charge density q.e.d. 

To understand all consequences of the Power-Zienau-Woolley transformation, 
one needs to consider a gas with more than a single atom or molecule. We there- 
fore consider a charge distribution in which the charge carriers are grouped into 
electrically neutral atoms or molecules. The nucleus of the mth atom is placed at 
a position R m , and the particle coordinates (electrons + nucleus) are denoted by 

Xrv^ . 



tie 

• 




m=l „ m = 3 




. m=2 

• 


• 


I* 

R S 'Af 
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In minimal coupling, the complete Hamiltonian of this system is given by 

5^ Vfcoui(m, m') + ^ Vcoui(m) . (5.20) 

m<m' m 

Here, E L is the operator of the transverse electric field. It is the same opera- 
tor which we up to now (in absence of charges) simply have denoted by E. In 
the presence of charges, the complete electric field takes the form E = E L — 
Wcou\{x). The Coulomb interaction between the particles of one atom is given 
by 

V Ct Jjn) = A 1" ( 5 - 21 ) 

On the other hand, the interaction between the particles of two different atoms can 
be written as 

V Coal (m,m')= V - 1 a - % - - (5.22) 

To construct the unitary transformation that connects minimal coupling to 
multipolar coupling, we employ the total polarization of the gas, which is a sum 
over the polariztion of the individual atoms: 

P(x) = J2P m (v) (5-23) 

m 

P m (x) = Y] du q am (x am ~ Rm) S(x - R m - u(x am - R m )X5.24) 
The Power-Zienau-Woolley -Transformation is then given by 



U = exp[- I d 3 x P{x) ■ A{x) ) (5.25) 



Explicit expressions for the transformation 

Our task is to find explicit expressions for the transformed Hamiltonian and other 
operators of relevance. Since U is a function of x am and A(x), it is obvious 
that because [A^x), Aj(x')} = [Ai(x), Bj(x')} = the vector potential, the 
magnetic field and the particle coordinates as well as the polarization transform 
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trivially, 

A = U ] AU = A (5.26) 

B = B (5.27) 

x am = x am (5.28) 

P = P. (5.29) 

All other fields can be treated by using the theorem exp(— iS)E exp(zS') = E — 
i[S, E] + ■ ■ ■ . We find for the electric field 

Ef(x) = Ei{x)- l -j "dV P^A^x'lEUx)} (5.30) 

= Ei{x) - l - [ dV Pj{x')—5iAx - x') . (5.31) 

fi J £o 

Recalling that the transverse (source-free) part of a vector field P is given by 
(TP)i(x) = J <fV Pj(x')5±(x — x'), we can derive 

E L (x) = E ± (x) - — P L (x) . (5.32) 
This equation has some interesting consequences. If we write it in the form 

e E ± (x) = e Q E L {x) + P ± (x) = (5.33) 
it becomes obvious that the old operator EoE 1 of the electric field in the new 
picture plays the role of the dielectric displacement field D . In other words: if 
one uses minimal coupling, the electric field is given by E L and the dielectric 
displacement by EqE 1 - + P ± . In dipole coupling, the electric field is given by E 

and the dielectric displacement by e E + P = EqE l . Since the operator E L 
has a simple form, it is also used in multipolar coupling to express the interaction. 

it remains to calculate the transformation of the particle momenta. One easily 
finds 

= Pam + hV am S. (5.34) 

The calculation of the gradient will be done below. In summary we obtain for the 
transformed Hamiltonian the expression 



# = T,J2^r(P<> m + W a7ri S-q am A(x a J) 2 



+ £ 4 fd 3 x h^-P^ + SB 2 )- 



2 J " W 

V Coui(m, m!) + ^coui(m) (5.35) 

m<m' m 
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Calculation of the polarization integrals 

This expression can be considerably simplified. We start with the second line: 

/ d 3 x (t) 1 - P^f = J- f d 3 x ((D 1 ) 2 - 2D 1 • P x + (pM (5.36) 

= -L [d 3 x ((D ± ) 2 -2D ± -P ± + (P-P^) 2 ) , 

2eq J \ / 

where P" denotes the longitudinal part of the vector field (i.e., the Fourier trans- 
form of the field is parallel to the wave vector, not transverse). We can now make 
use of the following relation, 

J d 3 x R x (x) ■ S(x) = J d 3 x R ± (x) ■ S ± (x) . (5.37) 

To prove it we first show that the transversalization operator is idempotent. By 
definition, we have 

(T 2 R)i(x) = J d 3 x'S^(x - x') J d 3 x"5f k (x' - x")R k {x") . (5.38) 



With 



ij k 2 J \ jk (fc') 2 



5i{x-x") (5.39) 



it follows that 



(T 2 R)i(x) = j d 3 x"5i(x-x")R k (x") 

= [TR)i{x) (5.40) 
In addition, T is Hermitean with respect to the standard scalar product: 

J d 3 x R{x) ■ (TS)(x) = J d 3 x Ri(x) J d 3 x'5j j {x - x^S^x') 

= J d 3 x'{TR)(x') ■ S{x') (5.41) 
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Relation (15.371) then follows from 

J d 3 x R x (x) ■ S(x) = J d 3 x (TR)(x) ■ S(x) 

d 3 x {T 2 R){x) ■ S{x) 
d 3 x (TR)(x) ■ (TS)(x) 
d 3 x R ± {x)-S ± (x) . 
Eq. (15.371) is helpful to simplify the integral over the polarization: 



(5.42) 



1 

2To 



d 3 x P 2 -2P-P !I 



= -L [d 3 x fp 2 -2(pl') 2 + (pl') 2 
2^o J ^ 

= ^~Jd 3 x (P 2 -(PH) 2 ) (5.43) 



The integration over the longitudinal part can be calculated exactly. Since the 
source of this term is the charge distribution, divP = divP" = —q, it is reason- 
able that it is related to the Coulomb potential. To demonstrate this we first rewrite 
the latter in the form 



1 



Atxeq 

-1 

4vre 
1 



Aire 



d 3 x' 


q(x') 


\x — x'\ 


d 3 x 


divP'^aj') 


\x — x'\ 


d 3 x' 





\x — X 



(5.44) 



The gradient of the Coulomb potential, which simply is the negative longitudinal 
electric field, we find 



d 3 V Cou{ {x) = J-^VPf^)^^^- 

fd 3 x' pU^djdi-. — - — r (5.45) 
J \x — x'\ 



-1 



4ire 



To calculate the double derivative of the integrand, we use the Fourier trans- 
form 

e -ikx ^ 

d 6 x — — — = 

\x\ k 2 



(5.46) 
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To prove this one can use spherical coordinates, with the z-axis along k. One then 
finds! 

ikx roo — ikr cost!) 

d 3 x- 

r 

-ikr Akr\ 



poo 


r 2 dr dcos'd dip 


Jo 




= 2ir j 


f°° 1 
f rdr — — (e 
o — ikr 


2m 
~k 


f v 

{ n5{k) — i- — 
I k 











V 

Ti5{k) + i— 
k 



(5.47) 



We thus find for the double derivative 



1 „ „ f d 3 k 47re ik - {x - x,) 



9 A- t 77 = d A 



x-x'\ 3 J (2tt) 3 k 2 

Att I rf3fc r ik.( x - x ') k i k j 
1 (27T)3 e k 2 



-4tt (Sij6{x - x') - 8±{x - x')) (5.48) 



so that 



djVcoAx) = J d 3 x' P\{x')(-An) (5ij8{x - x') - 5±(x - x')) 

= -P]{x) (5.49) 

The longitudinal electric field is therefore completely characterized by the polar- 
ization, so that 



^-Jd\{pW) 2 = -| J d 3 X (V^cou,) 2 

d 3 x VcouiAVcoui 



£o 
2 



-- / d 3 x VcouiQ , (5.50) 



3 The omission of the principal part V in the last line is based on the following considerations. 
Like every distribution, V/k is only defined as an integrand and needs to be multiplied with a 
test function (which has to be infinitely often differentiable, normalizable, without poles). In the 
second-to-last row, an additional factor of 1/fc appears, which is not a test function. Therefore, 
the result is only properly defined if one considers test functions which go at least linearly against 
zero for k — > 0. In theoretical physics, however, it is customary to use this equation without caring 
about test functions. The principal value is then meaningless, since it does not result in a finite 
value if the divergence is of the form 1 /k 2 . It is therefore simply omitted. 
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where the last line follows from the Coulomb potential obeying Poisson's equa- 
tion. Inserting Eq. (15.151) for the charge density, one arrives at 

- — fsx (pi') 2 = -- y y — — r (5.5D 



In this expression, we have to remove the terms with indentical particles since they 
diverge. The rest is symmetric under exchange of primed and unprimed indices, 
so that 



__L/ A(pl)2 = _£ £ 

m<m' a m ,c 

-E E 



4ire Q \x am - av J 

m<rn' a m ,a , m' 
m 

Qatm Qa' m 



„ , A-KEq \x am - x a > \ 

m a m <a' m 1 m m 

= - y V(m,m') - yV(m) (5.52) 

m<m' m 

The integral over the longitudinal part of the polarization therefore exactly 
cancels the Coulomb interaction between the atoms. This is the reason why we 
decomposed the polarization into the complete and the longitudinal polarization 
in Eq. (I5.36I) . The remaining integral over the complete polarization field can now 
be decomposed into the individual contribution of each atom, 



In Eq. (15.241) one can see that the polarization of the mth atom is strongly localized 
around the atom itself. Since in an atomic gas the distance between the atoms is 
much larger then their size, one can neglect terms with m ^ vnl and thus deduces 



in 



^yjd'x ((Pi? + (Pif 

2^E/ d 3 x(Pi) 2 + yV(m) (5.54) 
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Hamiltonian and multipolar expansion 

Inserting all calculations in the transformed Hamiltonian (15.351) one obtains 

H = ^ H m + H md + H m{ (5.55) 

m 

H ™ = E Sr- + v ^ + 2T / rf3x ^ (5 - 56) 

^rad = ^~ I d z x (£> x ) 2 + [d 3 x B 2 (5.57) 
# int = -— / d 3 x P 

+ S ( nam ' Pam + Pam ' n °" n + (nam)2 ) (5 - 58) 

m,a m m 



with 



^J^-^A^J (5.59) 

ft(V Qm )< / d« <?a m Oa m - H m ) ■ A(R m + u(x am - R m )) - q am A t (x am ) 
Jo 

J du q am (^Ai(R m + u(x am - R m )) 

+(x am - R m )jud i A j (R m + u(x am - Rm))) - q am Ai(x am ) 



The vector n am can be brought into a form which is explicitly independent from 
the potential. To do so one uses 

duud u A(R m +u(x am -R m )) =A(x am )- / du A(R m + u(x am -R m )) 

Jo 

(5.60) 

which can be proven by partial integration. Introducing the notation r am := 
x am — R m it follows that 



[fir 



I du q am (Ai(R m + ur a J + {r am ) judiA^Rm + ur a J) 
Jo 

- du q am (ud u Ai(R m + ur am ) + Ai(R m + ur am )) 
Jo 

/ duq am ((r am ) j ud i A j (R 

m ~l~ UT Qm ) u { r a m + Ur a m )) 

Jo 

/ du q am u{ r * m )j {diAj(R m + ur am ) - djAi(R m + ur am )) (5.61) 
Jo 
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and therefore, because of diAj — djAi = Ei^B k , 

n a m = q am / u du (x am - R m ) x B(R m + u(x am - R m )) (5.62) 
Jo 

This is the final form of the multipolar Hamiltonian. We emphasize again that 
no approximation has been made aprt from the assumption that the total polariza- 
tions of two different molecules have little overlap. It is therefore fully justified to 
use this coupling instead of minimal coupling. The most important features of the 
Hamiltonian (15.551) are as follows. 

• It is explicitly gauge invariant since it only contains the electric and the 
magnetic field. 

• The Coulomb potential between different atoms has disappeared. The inter- 
action between the atoms is now exclusively mediated through the retarded 

dynamical fields D and B. 



5.3 Dipole approximation and dipole coupling 

The typical size r of an atom is in the order of Bohr's radius (?» 0.5 x 10^ 10 m). 
Since optical wavelengths A = 2ir/k are in the range of several hundreds of nm, 
it is justified to restrict the multipolar expansion of the interaction Hamiltonian 
after the first few terms. The reason is that multipole moments of order I are 
proportional to (kr) 1 and kr <C 1 is valid, see Eq. (I1.18I ). We therefore keep only 
the lowest order terms and obtain for the polarization 

p(x) = Y Pm ^ 

m 

= Yl Yl qa ™ ~ Rm "> / duS ( X ~ R m ~ U(x am - R m )) 

~ Yl Yl qa ™ ~ Rm ^ I du5(x-R m )+0(\x am -R m )\ 2 ) 

= Y 5 ( X - R m) Y q «™ ( Xa ™ ~~ ^ 

= Y S ( X ~ R m)d m , (5.63) 

m 

where d m is the dipole moment of the mth atom/molecule. Analogously one finds 
for the second coupling term 

n a m ~ q am I u du (x am - R m ) x B(R m ) + 0{\x am - R m )\ 2 ) (5.64) 
Jo 
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and thus 



u du (x am - R m ) x B(R, 



-q am L am B(R r 



(5.65) 



This is the dipole approximation. The associated interaction Hamiltonian is called 
dipole coupling and given by 




(5.66) 



where \± m = J2 am q am L am /(2M am ) is the magnetic dipole moment that is asso- 
ciated with the orbital angular momentum of the electrons. The spin contribution 
simply doesn't appear here because it is omitted in our treatment of the Power- 
Zienau-Woolley transformation. It would of course be present if we start from 
particles with spin. 

We remark that only the complete multipolar Hamiltonian is unitarily equiv- 
alent to the minimal coupling scheme. As soon as approximations are made, the 
use of the different interaction forms leads to slightly different results. We will 
consider an example for this in the next chapter. In quantum optics, using dipole 
coupling usually results in slightly better results. 



5.4 Selection rules for atoms 

One of the most important questions in atomic physics and quantum optics is 
which levels are coupled by a given electromagnetic field. To examine this we 
consider a single atom whose nucleus rests at the origin (R m = 0, it is a very good 
approximation to neglect the center-of-mass motion in the derivation of selection 
rules). Since all the other particles are electrons, the electric dipole moment can 
be written as 

d m = ^ ^ lotm i X ce m Rm) 

I > m 

= ~ e ^ ] X a m 

=: -er (5.67) 

r is the sum over the position of the electrons relative to the nucleus, in the 
following we will omit the index m, since we consider a single atom only, d m — >• 
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d, a m — > a etc. Here, the sum includes only the electrons. If there are any 
contributions from the nucleus, they will be added explicitly. 

We have seen in chapter[4]that atomic states are, to a very good approximation, 
eigenstates of the total angular momentum 

F = s K + ^(X Q + s a ) , (5.68) 

Q 

with the atomic state = \F,m F , ■ ■ ■). Since the interaction is mediated 
through the dipole moment, the transition amplitudes between \F, m F ) and \F', m' F ) 
are proprtional to 

d F , mF ;F',m' F ■= (F,m F \ - er\F' \m' F ) . (5.69) 

The derivation of these matrix elements is surprisingly simple and elegant if one 
employs the concept of Tensor operators, which is summarized in appendix ll0.2l 
In particular, a Tensor operator of order 1 is also called a Vector operator R q 
(q — — 1, 0, 1). In the notation of App. 110.21 the spherical vector components R q 
are given by R q = Ti )(? . These operators are defined by their commutator relations 
with a given angular momentum operator L: 



[L + , R + ) 


= 


[L_,i2_] 


= 


[L+,R-\ 


= V2hR 


[L-,R+] 


= V2hR 


[Lg, R±\ 


= ±hR± 


[L z , Ro] 


= 


[L±, Rq] 


= V2hR 



Here we have defined L± = L x ± iL y . Equivalent to these relations is 

[L a ,R b ] = ihe abc R c (5.71) 

if the Cartesian components R c are related to the spherical components by R x = 
— i? + )/v / 2 and R y = i(R_ + R + )/y/2 as well as R z = R . An example 
for a vector operator is the position operator of a Schrodinger particle, with L the 
orbital angular momentum operator. However, it is an easy task to prove that the 
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dipole moment d is a vector operator, too: 

[F a ,d b ] = [(s K ) a + ^(L a + s a ) a , -e^2(xp) b ] 

a f3 

a,/3 

= ~e S a pihE abc (xp) c 
= -iefe a6c y^(a; 0! ) c 

a 

= ihe abc d c (5.72) 

To a large degree, the matrix elements of vector operators can be calculated 
algebraically. 

±h(F,m F \d ± \F',m' F ) = (F,m F \[F ' z ,d±]\F' \m' F ) 

= (F,m F \F z d±-d±F z \F',m F ) 

= h(m F — m' F )(F,m F \d±\F' ,m' F ) (5.73) 

It follows that (F, m F \d±\F', m' F ) can only be non-zero if Am := m F — m' F = 
±1. Since d z commutes with F z one can show that 

= (F,m F \[F z ,d z }\F',m' F ) 

= h(m F - m' F )(F,m F \d z \F',m' F } . (5.74) 

In summary, we obtain the selection rules for the magnetic quantum number 

(F,m F \d z \F',m' F ) only if m F = m' F (5.75) 
(F, m F \d+\F', m' F ) ^ only if m F = m' F + 1 (5.76) 
(F, m F \d-\F', m' F ) ^ only if m F = m' F - 1 (5.77) 

Further statements about the matrix elements of d can be derived with the aid of 
the 

Wigner-Eckart Theorem (here only for vector operators): The spherical matrix 
elements 

{eL, do, d + } := {(d x — id y )/y/2, d z , —(d x + id y )/y/2} of a vector operator are 
given by 



(F,m F \d q \F',m' F ) = (F', m' F ; 1, q\F, m F , (1, F'))^ 



V2F + 1 



(5.78) 
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The factor T>(F, F') is called reduced matrix element. 

This very important theorem states that the matrix elements of vector operators 
(and tensor operators in general) are proportional to Clebsch-Gordan coefficients 
(see Appendix 110.21 where also the proof of the theorem is presented). 

The Wigner-Eckart theorem has an intuitive interpretation. Electric dipole 
coupling is associated with absorption or emission of a photon, which has spin 
one. Since the total angular momentum is conserved, one has to add the atomic 
and photonic angular momenta in the way described in Sec. ll0.2l 

The Wigner-Eckart theorem does not make any statement about the form of 
the reduced matrix element X>(F, F'). For a given system one can derive it by 
explicit calculation for a given transition m,m'. However, this is usually very 
difficult and demands considerable numerical efforts. One therfore either uses 
the experimental value or the following estimate. The dipole moment operator is 
defined as electron charge times electron distance from the nucleus. Hence, the 
reduced matrix element should be of the order of eao, where ao is Bohr's radius. 
This rule of thumb gives often very reasonable estimates. 

Another important question is the physical meaning of the spherical compo- 
nents d q of the dipole moment. To shed light on this we enumerate the polarization 
vectors for left/right-circularly polarized light and linear orthogonally polarized 
light as 



1 

3 V 



e±i = £± = =F— (e x ± ie y ) , £q = e z (5.79) 



It should be noted that these vectors are orthogonal in the sense of e* • e q > = 5 qtq i. 
With these vectors we can express the dipole vector in its spherical components, 

i 

d = d x e x + d z e z + d z e z = ^ d i £ l ■ ( 5 - 8 °) 

<?=-! 

We thus can rewrite the interaction Hamiltonian for an atom with a ground state 
multiplett F g and excited state multiplett F e by using the Wigner-Eckhart theorem, 

H m = -d D 

- ^ \F e ,m e )(F e ,m e \d q \Fg,m g )(Fg,m g \ e* q ■ D+H.c. 

m g ,m e ,q 

= "7p=Y ^ \F e ,m e )(F g ,m g \(F g ,m g ;l,q\F e ,m e ,(l,F g ))e* q -D+ 

v e m g ,m e ,q 



Hence, left-circularly polarized light of a laser beam running along the quantiza- 
tion axis (= z-axis), i.e., D = e + exp(ikz) exp(-iujt), couples to d + and therefore 
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requires m e 



m, 



+ 1. Likewise, e = e z couples to do and implies m e 



in 



9- 



The corresponding laser beam needs to propagate perpendicular to the quantiza- 
tion axis, since its wave vector is orthogonal to e z . □ 



F z Eigenzustande 
m=-l m=0 m=l 




F=l 



m=0 



F = 



m g ±l 



We have to emphasize that circularly polarized light only implies m, 
if the laser beam propagates along the quantization axis. For another direction the 
laser light can generally couple to all spherical components d q . For instance, a left- 
circularly polarized beam along the x-axis has polarization e = (— e y — ie z ) / \/2. 
Because d ■ e = —id /\^2 — i(d + + cL)/2 it couples to all spherical components, 
so that m e = m g + i, % — — 1, 0, 1 are possible. 

To conclude this section we summarize the general dipole selection rules, 
which can be derived from the properties of the Clebsch-Gordan coefficients: 
Am = 0,±1 
AJ = 0,±1 

For J' = J the transition m = — > m! = is forbidden. 



Strictly speaking, this is too strong a statement. Since real laser beams are always focused, 
they correspond to a wavepacket to which also wave vectors contribute which are not exactly 
parallel to the axis of propagation. A real laser beam therefore always has a small longitudinally 
polarized part. This part, however, can usually be neglected. 



The Interaction between Atoms and Light 



Chapter 6 

The two-level model for atoms 



If atoms in their ground state \g) with energy E g are exposed to a laser beam 
composed of photons with energy ftw, they can absorb the energy of the photons 
and make a transition to an excited state |e) with energy E e . Energy conservation 
implies E e « E g + ho. If no other atomic energy levels are close to E g + hu, one 
can neglect all states but \g) and |e). One then can think of the atom as a neutral 



e') 



Huj 



Figure 6.1: Absorbtion of a photon by an atom 



particle with just two internal levels \g) and |e); it then resembles for instance 
a neutron with spin ~. Models with two levels appear frequently in quantum 
mechanics, for instance for 

• real spin \ particles 

• photons with polarization 



the band edges of valence band and conductor band 

atoms and molecules 
ions 



L 

V 
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• Josephson effect in superconductors 

• interferometer 

• qbits 

The mathematical content of this chapter can also be applied to the other physical 
systems, but we will focus on atoms and molecules. 

Two-level models provide an intuitive picture of complex physical phenom- 
ena. Instead of their simplicity they often lead also to excellent quantitative re- 
sults. 

6.1 Derivation of two-level systems 

Example: hydrogen atom The eigenstates (without spin) have the form \nlm). 
In absence of a magnetic field they are degenerated with respect to m. 

— — — 2P 

2S 

— IS 
Figure 6.2: Sketch of hydrogen energy levels 

If laser light of frequency uj ps (E 2 p — E\s)/ h is used, all multiplets except for 
IS" and 2P can be neglected. We then obtain a system with four states. Sometimes 
systems consisting of two multiplets are already called 2-level systems, since only 
two energy levels are involved. In the presence of hyperfine degeneracy such 
systems can include a considerable number of states. To achieve a real two-level 
system (with two states instead of just two energy levels), one has to use laser 
with a suitable polarisation and to keep an eye on spontaneous emission. 

Examples for two-level systems 

1. Example: 

Absorbtion und Emission 

If one pumps with circularly polarized light, only the two states to the right 
may form a two-level system. 
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^ \e,F',m F ,);F' = 0,1, 



F = 2 1 , 

_ _ _ F = i J \9,F,m f ) 

Figure 6.3: Hyperfine splitting of Rb 87 . Even if light is used that is resonant with 
two energy levels (for instance F — 1 — > F' — 2), one generally has to consider 
many states. 



2. Example: 

r - 3 

e 2 

Absorbtion und Emission 

r _ 1 

9 2 

Mechanism: 

• |m 9 = |) is only coupled to |m e = |) and forms a two-level system 

• \m g = — |) is coupled to \m e = |). This state decays spontaneously 
into \m g = —\) and \m g = |). Hence, for each spontaneous emis- 
sion population is transferred to the real two-level system \m g = \), 

\m e = |). 

• After many spontaneous emissions almost all atoms are transferred to 
this system. 




6.1.1 Representation of operators 

In the rest of this chapter we will only consider a two-level system with states 
\g), \e). The representation of any operator acting on this system can be found as 
follows: 

oo oo 

O = \n) (n\ O \m) (m\ ~ \n) O nm (m\ 

n=0 m=0 n,m=e,g 



mit: O nm = (n\d\m) = (° ee °A 

\ u gz u ggJ 



The most important operators are: 
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Identity operator 1 

' E, 



1 

1 



Hamiltonian H, 



o 



£ j since \e), \g) are eigenstates of H . 



Electronic position and dipole operator x oc d = —ex with matrix elements 

(e\d\g) := d eg = d* ge 

(e\d\e) = (g\d\g) = 

The vanishing of the siagonal elements is a consequence of |e) and \g) being 
eigenstates of the parity operator, and of d having parity —1. One therefore 



has 



•es \ 0/ e » \1 01 \d' I \x' 

(6.1) 

where x eg is the matrix element of the relative coordinate of the electrons. 

Electronic momentum operator We first proceed as with the position operator 
to find (g\p\g) = and (e|p|e) = 0. To derive the off-diagonal elements 
we use the following trick: 

[x, H ] = —p (6.2) 
This follows from H = J^a^Pl + J2 a ,/3 V ( a ^) and P = E Q Pa> 



[x at H ] = 

a' 6 
ill 



m 



-p, 



a 



e 



Hence 



(e\p\g) = -Tjr (e\[x,H ]\g) 

= ^{E e -E g )(e\x\g) 

so that we find for the matrix representation of the electronic momentum 
operator 

m e u ( -id eg \ 

^ = . (6.4) 

n 
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6.1.2 Coupling between a two-level atom and the electromag- 
netic field 

Consider an electromagnetic wave with a positive frequency part of the vector 
potential of the form 

A {+) (t, R = 0) = A Q e~ lu}Lt R. Center-of-mass of the atom (6.5) 

The interaction Hamiltonian in dipole approximation reads (E = —A) 

Hint = —d ■ E 

= -d(E {+) + 

= -iu L (d eg \e) (g\ + d* eg \g) (e|) {A e^ - A* e^) (6.6) 
For comparison we present the minimal coupling Hamiltonian: 

H' mt = —p-A = -zuj (d e9 |e> (g\ - d* eg \g) (e|)(A e— ^ + Atf**) (6.7) 

Tfl e 

Apart from various signs, which can be transformed away by a different choice of 
the phase of the basis vectors and a shift in the center-of-mass coordinate, the main 
difference is the occurence of the factors ujl and u , respectively. The reason for 
this difference is that we only consider a part of the full Hamiltonian (that includes 
all levels), and only the full Hamiltonians are unitarily equivalent. In Fig. l6.4l the 
unitary transformation is represented as a rotation in Hilbert space. This has the 

Rest') 

<n 

Figure 6.4: Unitary transformation = U\il>) 

following meaning: the unitary Power-Zieman-Woolley transformation "rotates" 
the states in Hilbert space. The two-level approximation is a projection onto a two- 
dimensional subspace. Depending on whether one performs this projection before 
(|e) and \g)) or after (|e ; ) and \g')) the Power-Zieman-Woolley transformation, one 
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obtains different results since one considers different subspaces. 
In practical calculations this difference is unimportant since we are only interested 
in the case lu ps uj l . Q In the following we will use the — d ■ E coupling. We 
define H Rabi frequency 



Q :— jdeg ■ Eq — jdeg ' A UJ L 



(6.8) 



with Eq := Aquo l and the less important abbreviation 

1 

h 

We arrive at 



il = -d eg El (6.9) 



H int = -ih |e) (g\ (fie - ^* - tte iuJLt ) + ih \g) (e| (tt*e iuJLt - tt*e~ luJLt ) (6.10) 
The unperturbed Hamiltonian is given by 

E e \ 1 , „ 1 — i \ .i 1 



By choosing a particular value for E g (the zero-point of energy can be chosen at 
will) one can represent H in various ways: 

E e + E g = # = ^^ 2 , w = Ee I Eg (6.12) 



o oy = '"~ u ^0 o 

We will frequently make use of the rotating-wave approximation. This is a widely 
used and well-justified approximation which considerably simplifies calculations. 
To derive it, we first switch to the interaction picture. Thereby the temporal evo- 
lution associated with H is already incorporated in the wave function: 



e 



-*ffot/*|^\ = C /-u,\ (614) 



The transformed Schrodinger equation becomes ih \ip) = H \ tp) and the Hamilto- 
nian is transformed according to 

H = -ihrftl + U\H Q + H mt )U = U^H mt U 

~ , n (6.15) 
= -ih \e) (g\ (ne^ {WL ~ wo)t - Sle l ^ L+Ua)t ) + H.c. 



1 ojq f» ujl is a requirement for the two-level model to be valid. 

2 Very often the definitions used in the literature differ from our convention by a factor of two. 
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The term e^-c+^o)* is fastly oscillating and therefore averages to zero during the 
evolution. This leads to the following simple form of the Hamiltonian in interac- 
tion picture, 



H ww (t) = -ihe- z/xt n \e) (g\ + iKD?e lIxt \g) (e 



(6.16) 



where we have used the laser's detuning A = ujl — cu . 
Transforming back to the Schrodinger picture we find 



H 



h 











(6.17) 



By another unitary transformation we can achieve that H only depends on the 
absolute value of VL. This transformation is achieved by multiplying |e) with a 
phase factor e lLp , where — iVt = \Vt\e lLp . 



U : \e) 



U 







1 

This leads to the final and simple form for the Hamiltonian, 

|fi|e _<Wi * 
n\e iuJLt 



h 



(6.18) 



(6.19) 



6.2 Rabi oscillations and Landau-Zener transitions 

We will now deal with the time evolution of a state with Hamiltonian H as dis- 
cussed in the previous section. 

Without the time dependent factor e Wht this equation would be easily solvable. 
Luckily, we can remoce this factor by a unitary transformation: 



The transformed equation becomes 



ihd, 



~2~ a " + -2~ a * + h W ) ) U I (6 " 22) 
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We therefore find for H: 



H=h 



-f 



(6.23) 



Such a time-independent Hamiltonian is only possible in rotating-wave approxi- 
mation. The counter-rotating terms would still depend on time after transforma- 
tion Eq. (I6.21I ). We now can easily solve the Schrodinger equation: 



r Mt)\ =e -its/n Ayo) 



(6.24) 



Let us now rewrite H. Since the Hamiltonian is Hermitean and traceless, there 
are only three independent parameters. 



H = hn ■ (T. 



/lmQ\ 
n = I Refi J 



The expression e itn cr can be expanded into a Taylor series, 

oo 1 



-itn-er 



l=0 



The product of two a matrices is given by 



OiOj = ie ijk a k + (%J 
[n ■ <x) 2 = riiaiUjaj = ninj{iSij k a k + 
f = \n\ 2l i 



n 



[n ■ (7 

We thus can bring e~ %tn cr into the form 

\2l 



-itncr 



l=0 



(2()! 



i=o 



-it) 



2I+1 



(2Z + 1)! 



: n ■ a 



-itn-cr 



cos(|n|t)J —i—, — j-sin(|n|t) 



n 



We introduce the notation 



w •= n 



mi 



(6.25) 



(6.26) 



(6.27) 
(6.28) 
(6.29) 



(6.30) 
(6.31) 



(6.32) 
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which only depends on the modulus of Vt. For an atom that initially is in its ground 
state, the solution becomes 





— - sin(wt) 
cos(wt) — i£- sin(wi) 



(6.33) 



The probability amplitude ip e (t) to excite an atom oscillates with frequency w, 
this phenomenon is called Rabi oscillations. We have 

V e (t) = \Mt)\ 2 = ^sin 2 (wt) (6.34) 
The maximum value of V e is reached for wt = n/2 + nn and given by 

mis 

cc , ' — (6.35) 





n 


2 


(f) 2 + 


n 


2 



If the laser is not detuned against the resonance frequency, A = 0, complete 
excitation (V e = 1) is possible. V e depends strongly on the detuning, 

A»|ft| T 3 e«^-<1 (6.36) 

The excitation probability therefore depends on the ratio of Rabi frequency and 
detuning. H These results are of very high relevance for atom optics, quantum 
optics, electron spin resonance, nuclear spin resonance, Interferometry, and many 
other areas. 

With laser pulses of a particular duration one can prepare well-defined atomic 
states. The so-called ideal | pulse El creates a 50-50 superposition of ground and 
excited state: 

A = , = ^ tf = e — = -^(} - 1 ) (6.37) 

The ideal n pulse transfers all electrons to the excited state if they were prepared 
in the ground state, or vice versa: 

A = t = ^^ U =(l~o) (63S > 



3 It is often claimed that the excitation probability depends on the ratio of Rabi frequency and 
spontaneous emission rate. We will see later that this ratio is also important for small detunings. 
The ratio of £1/ A is usually more important, though. 

4 In practise the detuning of a laser against the resonance frequency of an atom depends through 
the Doppler effect also on its center-of-mass motion. At room temperature this makes an ideal ir/2 
pulse impossible. For cold atoms and with some tricks one can achieve such transitions with very 
good efficiency. 
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To visualize the states of a two-level system one can use the Bloch vector, 
whose components are given by (cr) = (v, w, u). One then has 

u = (cr z ) = |^ e | 2 - \ip g \ 2 

v = (a x ) = + Vg^e (6.39) 

The Bloch vector of a pure state is of unit length, it always moves on the Bloch 
sphere of radius 1 around the origin. 




Figure 6.5: The Bloch sphere 
For a non-detuned laser beam Rabi oscillations have the following form, 

ft) - KS?) ^ 

The corresponding Bloch vector becomes 

v\ /-sin(2fit)\ 

w = (6.41) 

u) \-cos(2fii)/ 

For a 7r pulse, t — the Bloch vector moves on a semi-circle, see Fig. 16.51 

Interferometry with two internal states: an example for the application of n 
and it/ 2 pulses two-level systems and Rabi oscillations can be used to realize 
an interferometer. The two states thereby play the role of the arms of the interfer- 
ometer, Rabi oscillations provide beam splitters and mirrors. Sketch of set-up: 



6.2 Rabi oscillations and Landau-Zener transitions 




l^o) = \g) 



IV> 2 ) 





























-( 






y/2 



(\g) 



e* |e» 



1^3) 





'2 |^2) 


1 










1 - 


U 


1 + e <v 



1 

71 



3«V 



^-(|0)(l + e *) + |e) (l-e*)) 



Detect the number of atoms in |e): 
P e = K^|e>| 2 



1 



(1 - e itp ) 



(1 — cosy?) 
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6.2.1 Landau-Zener transitions 



Landau-Zener transitions are not as important as Rabi oscillations, but they nev- 
ertheless explain many relevant phenomena. The difference to Rabi oscillations is 
that the detuning is allowed to vary in time, A = A(t) = A t 



H = h 



A 

2 

n 



2 1 



where we use the rotating- wave approximation and assume G K. for siplicity. 
Equations of motion: 



iljjg 



2 

-A tl/;g + Qtpe 



Differentiate ip g with respect to time to obtain 

ii>g 

Inserting (16.421) results in 



1 • 1 • ■ 

-A 1p g + -A tVp g + tttpe 



fa 



^ + t -^+ l -A 



(6.42) 



Introduce the new variable y = t y A 



y 2 fi 2 i , 



This is the standard form of the differential equation for parabolic cylinder func- 
tions (see El). 

The solution can be expressed in terms of confluent hypergeometric functions 
and, for tp g (— oo) = 1, typically look like this: 



\%\\ 




0.8 


r 


0.6 




Va~ « n M| 

J 





to. 



- 2 
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-4-2 2 



The asymptotic excitation probability P (+oo) can be calculated analytically 
and is given by 

P (+oo) = 1 - |^(+oo)| 2 = 1 - e " 2A o 

Landau-Zener transitions are of importance in many different areas. They have 
first been discussed in moleular physics: 




• Light of frequency lul is only resonant for a particular distance R between 
the nuclei. 

• The time evolution of the electrons is described by states \e(R)), \g(R)) and 
is much faster than the nuclear motion (since M p 3> M e ) 

• In the equation of motion for the electrons, the distcance between the nuclei 
then only enters as a (slowly varying) parameter: 

in (i>e\_( -AE(R(t))/2 n \f ^ e \ 

UJ V n AE(R{t))/2 J\ip g J 

• Around Rq one has AE(R(t)) = E e - E g w AE(R ) +E' (R - R ) 

hjJ L KiV Q t 

A(t) =oo L - AE(R(t))/H = -E'vot/H 
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Landau-Zener transitions are also important for STIRAP and Bloch oscillations. 



6.3 Dressed States 



Dressed States are generally eigenstates of the full Hamiltonian H, including an 
interaction with one or few light modes. The simplest example is the interaction 
with a laser beam as described above, 



Hip = Eip 



The corresponding eigenvalues are given by 



+ M 



(±) 



2w V i^/wi A/2 
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This phenomenon is called "avoided crossing": For f2 = the eigenstates are 
ip e and ip g ; their eigenvalues are crossing at A = 0. For f2 7^ 0, the coupling 
of the states by the laser beam transforms this crossing into two non-intersecting 
hyperbolas with minimum enerhy distance 2fr|fi|. 

It is remarkable that, if one follows one of the hyperbolas, one evolves for E_ 
from ip g to ip e , and for E + from i\) e to ip g . This means that, for time-dependent 
A(t), one can make a transition from ij) g to ipe while always staying in an instan- 
taneous eigenstate of H(t) (see also the section on the adiabatic theorem below). 

This behaviour is closely related to Landau-Zener transitions, tunneling be- 
tween E- and E + corresponds to not exciting the atoms. For a slow rate of change 
an atom tries to stay on one of the hyperbolas for a Landau-Zener transition. 

The asymptotic transition probability can also be understood by estimating 
the tunneling probability between the two hyperbolas. According to the WKB 
approximation (see, e.g., Ref. (ITST). the tunneling probability between the hyper- 
bolas can be expressed as P = e ~ AEt / h , with ^ = 2 \ The typical time scale 
t for tunneling can be estimated by the time when the distance between the two 
energy levels has grown to 4 \ 

A = A t « 4 |0| 

4 101 
t = -4-^- 

A 

2jn[4|nj 8|o| 2 
=>. p = e A o = e A o 

6.3.1 Dressed states in a cavity 

The notion of a dressed states is more frequently used if the atoms are coupled 
to the quantized radiation field. The form of the corresponding two-level Hamil- 
tonian can be derived as before. The only difference is that the Rabi frequency 
becomes an operator: 

h = \d eg -E {+ \x) (6.43) 






A - 


-> 00 


1 


A - 


-> —00 


1 


A - 


-> 00 





A - 


-> —00 
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This results in the Hamilton operator 



H 



hjJi] 



-a z — 



d eg -E(x)\e) (g\ + d*-E(x) \g) (e 



rad 



Using 



T t —■} — li rr i — >i rax 

U = e 1 2 g 1 a 
we can transform this to the interaction picture, 



Hi 



int 



d eg \e) (g\ e*" ' + d* \g) (e\ ■ E(x, t) (6.46) 



iojot 



We have 



e(x, t) = iJd 3 kJ2 y 2(^ye/ kXe ~ iuJkt£kaaka + h - c - 



2(27T)3£ 

Visualization of the terms |e) (g\ a ka and |e) (g| aj^: 



|e><g|aka 



e><g|aka 



(6.44) 



(6.45) 



(6.47) 



Figure 6.6: The term on the right-hand side is somehow unphysical and is ne- 
glected in the rotating-wave approximation (RWA) 



In RWA, one therefore neglects the processes of absorbtion & de-excitation 
and emission & excitation. 

Having made the RWA, we transform back to the Schrodinger picture, 



H 



fojJn 



Schroed 



-o 7 



■J d 3 k ^hw k a} ku a ka + h J d 3 k ^ J|e) (g\ a ka g ka +U. c. 

a a 

(6.48) 
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where g ka is given by 



g ka = -id eg ■ e fcg J ofc/0 _, 3 . e lkx (6.49) 



2^(2vr) 3 e 

The Jaynes-Cummings model, which describes the interaction with a single pho- 
ton mode, has been particularly well studied. From Eq. (16.481) one finds 

H = h (jY a z + VlcJci + gcr + a + g*a~a ( j . (6.50) 

where a ± are defined by 

^(!!) (6.5D 
This model is meaningful if 



0/ V 1 



• only one photon mode is occupied (n a ^> n ka ), 

• only one mode of the cavity is resonant. 



During the last 15 years huge progress has been made in cavity QED. The fun- 
damental idea can be explained with the aid of an idealized cavity whose walls are 
perfect mirrors. The mirrors impose boundary conditions on the electrodynamic 
(quantum) field, which guarantee that only modes with wavelengths A n = 2L/n 
are allowed. The corresponding frequencies cu n = 2iic/\ n oc n are discrete and 
are only resonant for cu n ps (E e — E g )/h. Therefore, only modes with a certain 
value of n need to be taken into account. Typical values for microwaves are cav- 
ities with a length L ~ 1cm and n ~ 1; for optical frequencies, one usually 
employs cavities of length L ps 1 — 10 cm and n ~ 10 6 1. 

Real cavities are never ideal. A measure for their quality is their finesse Q, 
which represents the mean number of round-trips which a photon makes before it 
is lost. Very good cavities have Q ~ 10 6 . Because the mirror material of the cavity 
walls is usually slightly birefringent, a cavity is only resonant for one direction of 
the photon's polarization. 

If an atom is placed into a 3D cavity, its spontaneous emission is modified 
since it can only emit photons in certain modes (of which none may be resonant), 
and because it can reabsorb the photons since they stay close to the atom. 

Dressed states in a cavity can generally be represented by 

l^> = ^2 W n ) l e ' n ) + ^a( n ) \9, n) ) (6.52) 
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We are looking for the eigenstates that obey H \ip) = E\ip). This leads to the 
following system of equations, 

(E + ^ - nhw L )ip g (n) = hg*y/nij) e (n-l) (6.53) 



(S_^o_( n _l)^ £ )^ e ( n _l) = hg*V^iP g (n) (6.54) 



where we choose ip e (n — 1) for the second equation because it is this component 
which couples to ip g (n) (exciting the atom is accompagnied by the absorbtion of 
a photon!). Setting Q n := gy/n and A = lu l — lu we can write this as 



This is the same equation as for a classical e.m. field, so that the eigenvalues are 
easily found to be 




E n , ± =(n--)fcu L ±h\ (-) +|fi n | 2 (6.56) 



The energy also depends on n which means that the phases of states with different 
n evolve differently ( e -j£; ™.±*/fr). This leads to destructive interference between 
their amplitudes and to collapse and revival of the excitation probability. 
References: 0, 



6.4 Models with few levels, dark states, Raman tran- 
sitions 

The most common extension of the two-level model is of course the three-level 
model. Such a system can be achieved by choosing lasers in a suitable way. An 
example are transitions from J = 1 to J' = 1, where the transition from J' to J 
with m = m! = is forbidden by a selection rule. Shining two laser beams with 
polarization a + and ct_ on the atom therefore depopulates the state m = 0, J — 1. 
After this has happened, one obtains a so-called A system consisting of the states 
\J=l,m = ±1) and \J' = l,m' = 0). 

The Hamilton operator of this system in RWA reads 

/ E e hQ + e~ iuJ+t WL-e-^- 1 \ 

H = nn* + e iuj+t E + 

\ nnie™-* o E_ ) 

= EJ+l W)*e iw + t E + — E e (6.57) 

\ m*_ e-^- 1 E_-E e j 
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Figure 6.7: The arrows represent the transitions which are driven by the cr± light. 
The eavy lines represent de-excitation due to spontaneous emission. 



Again we perform a unitary transformation which removes any time depen- 
dence from H: 

$ ± = e _iw± V± (6-58) 
This leads to the time-independent Hamiltonian 



EJ (6.59) 






n + 


Q_ 


n* + 







n*_ 





A_ 



where we have set A± = u± — — h E± . 



6.4.1 Dark States 

Consider the case of two photon resonance: A + = A_. This means that for 
A ± 7^ the individual photons are not in resonance with the respective transitions 
|±) — > |e). However, we have h(u + — cu_) = EL — E + which follows from 
A + = A_ because of wegen 



f\w+ — E P + E A 



E e + E_ 



If the condition for two photonen resonance is fulfilled, there exists an eigen- 
state of the Hamiltonian of the form 



1+) 



(6.60) 



which is called a dark state. Here we have set Q 



\Q i I . The action 
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of the Hamiltonian on this state is given by 

/ Vt + \ / \ 
H\rj} D ) = h I n* + A+ tt^/tt 

A / \ J 

for A + = A_ = A. Hence, in this case \ipD) is an eigenstate of H with no excited 
component. The excitation probability is therefore zero, the atomic gas is for this 
(!) particular combination of light field transparent. H The transition probabilities 
interfere destructively. 

Simple special case: Q + = =>- \^ D ) = |+) 



le> 




l-> l+> 



If the Rabi frequencies are changed, \i/) D ) is coupled to the excited state again. 
For any set of Rabi frequencies, there is a bright state 

l«=Tr l+>+ if!-} <6 - 61) 

which is orthogonal to \^ D ) and is coupled to the excited state, H \ip B ) — Q |e) + 
A|Vb) 

l^u) and \iPb) form an alternative basis for the subspace spanned by |+) and 
|— ). Spontaneous emission from |e) can generally end in \^d) or \ipB) since the 
dark state is only dark for the two laser modes and is coupled to other light modes: 

\M \1> B ) 

This can be used to pump the atoms into \iI>d) (see the chapter on laser cool- 
ing). For A + 7^ A_ the dark state \i(> D ) is not an eigenstate of the Hamiltonian 



5 This is also the case when we have no two photon resonance. 
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anymore and will evolve into the bright state \ipb), which can be excited. The 
pump cycle is then broken. 

The physical meaning of dark states can be understood by considering a tran- 
sition from J g = L g = 1 to J e = L e = 0. Although this example is not very re- 
alistic (since L„ usually vanishes in the ground state) it displays the main features 
of dark states. Consider the form of the electronic states |e) = \L e = m e = 0) 
and \m) = \L g = 1, m g = m) in position space. Using spherical coordinates, 

(r sin $ cos ip 
r sin "d sin (p 
r cosi!) 

we find for the states 

(a;|-l) = y 1 ,_ 1 (i?,^)r/ fl (r) 

= re~ tlp sm$fg(r) = (x - iy)f g {r) 

(x\0) = r cos'^?/ 9 (r■) 
= zf g {r) 

(x\l) = (x + iy)f g (r) 

(x\e) = Y 0t0 fe{r) 
= fe(r) 

This results in the following matrix elements of the dipole moment -ex, 
d e -\ = (e| — ex | — 1) 

(x 
V 
z 

/fx 2 — ixy 
d 3 x f*f-i(r) xy - iy 2 
\ xz — iyz 

By symmetry under x — > —x' or y — ^ —y' one has 

d e ,-i = -e J &/:(r)/_! J -iy 2 



6 In realistic atoms one for instance has L g = 0, S — 1, L e = 1 so that J g = 1 and J e = for 
a A configuration. Though the angular momenta have just the reversed values as in the example 
above, selection rules imply that the arguments which will be presented in the following can also 
be applied to this case. 



134 



The two-level model for atoms 



Symmetry also leads to J f(r)x 2 = J f(r)y 2 , so that d e -i '- 



- 2 - -■ (l,-z,0). Analo- 

gously one can show that d e ,i ~ (1, i, 0) and d e )0 ~ (0, 0, 1). 

We now focus on the example of light that is liniearly polarized along the x 
direction, so that E ~ e x . Because 



V2 



V2t 



this is equivalent to using left- and right-circularly polarized light of same inten- 
sity. This leads to the Rabi frequencies 



n 



h 











i 


II 





I- 1 


1 




\ ) 





1 ^ 








—i 


II 





1- 


J 









The Rabi frequencies are therefore equal, so that the dark state (16.601) reads 



|^) = _(|+1) 



-1)) 



In position space this becomes 

1 



= V2iyf g (r) 



(6.62) 



is antisymmetric in y direction and corresponds to a p state along the y axis. 

I x P-Hantel 




With the knowledge of the spatial form of the dark state we reconsider the 
interaction between the atoms and the laser beams. Classically, an electric fields 
interacts with an atom by displacing positive and negative charges against each 
other along the direction of E: 
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The initially symmetric charge distribution thus becomes asymmetric. In quan- 
tum mechanics, one has a similar picture. E creates a superposition of |e) and \g). 
In our example, |e) is symmetric in position space, and all ground states are anti- 
symmetric. Lj A general superposition is not symmetric. 




The charge distribution of the electrons is given by q(x) = — e\il)(x)\ 2 . It is 

symmetric both for the excited state |e) and for each of the ground states. Only a 

superposition corresponds to an asymmetric, displaced charge distribution. 

▲ 

► 

v 

The superposition of symmetric and antisymmetric states is therefore the quan- 
tum mechanical description of the displacement of charges in an electric field. 

If we apply a laser beam (whose electric field is in the direction of its polar- 
ization) to normal atoms, it is possible to find a combination of an s and a p state 
which corresponds to a charge displacement along the direction of the electric 
field. For atoms in a dark state this is impossible: in our example the electric field 
is along the x axis, but the dark state (16.621) is antisymmetric along the y-axis. This 
means that for atoms in a dark state, the charge distribution cannot be displaced 
along the direction of the electric field. 



7 Just as a reminder: for real atoms the situation is usually inverted. 
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6.4.2 Raman transitions 

Consider the time evolution generated by the Hamiltonian 



VL + 

H = h ( Q; A+ (6.63) 

n* o a_ 



iip± = A±ip± + Q*±ijj e 



Ansatz: ip± = e lA±t ip± 



iip e = n + e~ iA+t ip + + fi_ e - iA -*^_ (6.64) 
# ± = n* ± e iA±t ijj e (6.65) 

Since \tp e \ < 1 one finds that \ip±\ < \Q±\ holds. The formal solution for ip e is 
given by 

^ e (t) = ^ e (0) -i [ dt' \n + e- iA+t '4> + (t') + n_e _iA - t '^_(t')) (6.66) 



We now focus on the case of large detunings |A±| ^> \£l±\. In this case, ip±(t) 



varies slowly compared to e iA± * since 



< IfLJ <c |A±|. We therefore can 



set ip± (t') ~ ^± (t) in the integral, which results in 

Ut) « ^ (e" 4A+ * - 1) + ^ (e-' A -* - 1) £_(t) 

Inserting this into Eq. (16.651) leads to 

4 ± = n* ± e iA ±t |5± (e""M - 1) Mt) + ^- (e-^ - l) 
or, explicitly for one component, 
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where fastly oscillating terms average to zero and therefore can be neglected 
against the slow terms. U hi summary, we thus can write down an effective two- 
level description for the two ground states, 

»£(£)-*■(£) 



with 



1"+| 2 n *+ n ~ c i(A+-A-)t 



The excited state is never significantly populated and could be omitted in these 
equations. This procedure is called adiabatic elemination of the excited state. 

The matrix elements of H e s are proportional to f2 2 /A. It is of fundamental 
importance that this expression can be large even when the excitation probability 
P e = f2 2 /A 2 is small. This makes Raman transitions an excellent tool to coher- 
ently manipulate atomic states, since a small excitation probability implies that 
loss of coherence due to spontaneous emission is strongly suppressed. The time 
scale of loss of coherence (decoherence) due to spontaneous emission is given by 
the excitation probability times the rate of spontaneous emission, ^yP e . Typical 
numbers for Raman transitions are f2 = 10 9 Hz, 7 = 10 7 Hz, A = 10 10 Hz. This 
results in fi 2 /A = 10 8 Hz and -fP e = 10 5 Hz. 

The effective Rabi frequency 

fieff := (6.69) 

is a measure for the strength of the coupling between the two ground states. We 
have omitted the index ± for the detunings, since for A_ ^ A + one runs into the 
problem that H ei{ is not Hermitean anymore. However, H ef f is nevertheless a very 
useful operator if one knows why this happens. The reason is that the adiabatic 
evolution of the excited state which was discussed above, 

V> e = ^ + (e-^ - 1) + (e-' iA -* - 1) (6.70) 

is actually the first order term of a Taylor expansion of ipe with respect to 
Thus, ipe is only correct to first order in if all other frequencies in the Hamil- 
tonian, Vt±, A + — A_ A±, are much smaller than the detunings. This implies 



8 In contrast to this, in Eq. (I6.65i we had fast terms only. Hence they produce the dominant 
effect and cannot be neglected. 



138 



The two-level model for atoms 



that we have to do the following expansion: 
1 1 



A_ (A_-A + ) + A 
1 1 
a! ~ A 
1 



_-_(A_-A + ) + ... (6.71) 



A++ 0((^)) (6.72) 

This means that to first order the expressions 1/A± are the same, if eff is therefore 
Hermitean to first order in Introducing the effective Rabi frequency as in 
( 16.691 ). one can represent _ff eff as 

H eS = h[ n , e _ i(A+ -A_)t ( 6 - 73 ) 



'eft" 



* 



This Hamiltonian is analogous to that that of a two-level system. 



6.5 Adiabatc theorem, STIRAP 

The method of STImulated Rapid Adiabatic Passage is a highly efficient tool of 
coherent manipulation which combines dark states and Raman transitions. It is 
based on the 

Adiabatic theorem: Let H(t) be a time-dependent Hamiltonian with instanta- 
neous eigenstates \n(t)), H(t) \n(t)} = E n {t) \n{t)). Then the evolution of a 
state \ip(t)), which initially coincides with the eigenstate \n(0)), (approximately) 
follows adiabatically \n(t)) if H(t) varies sufficiently slowly, \(n(t)\m(t))\ 
\E n (t) — E m (t)\/h). This means that 

\ip(t)) w e im \n(t)) (6.74) 

This is the "folk version" of the theorem. More elaborate versions which also 
present the correction terms can be found in [@J|40l. The theorem is based on the 
following observation. 

We want to solve ihd t \^(t)} = H(t) \ip(t)) with the initial condition |^(0)) = 
|n(0)). We expand the state in the instantaneous eigenbasis of H(t), 



W*)> = Ht)) 

It 



(6.75) 
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Inserting this into Schrodinger's equation results in 

ihd t \i/j) = ih^2{c n \n) + c n \h)} 



H\i(>) 

Y,E n {t)c n {t)\n{t)) (6.76) 



Multiplying both sides from the left with (m(t)\, and using (m(t)\n(t)} = 5 nm 
one arrives at 

ic m + i ^2 ° n ( m \n) = ~^° m (6-77) 



Ansatz for the coefficients c m (t): 

c m {t) = exp f-i 1^ dt' E m {t')/h \ ■ ~c m (t) (6.78) 
This leads to the follwing equation for c m , 

c m = -Y j ~c n (m\n) exp (-i jf dt' (E n (f) - E m (t'))/hj (6.79) 

If H(t) varies slowly, then (m\n) is small compared to the phase factor J (E m — 
E n ) for m n. We then can neglect these terms and arrive at 

c m « - (m\fn) c m (6.80) 

Because 

(m\rh) = — (rh\m) = —((m\rh))* (6.81) 
we find that (m|m) = ia(t) is an purely imaginary number, so that 

c m (t) = exp f-i a(t')dt'^j (6.82) 

Remark: the phase J a(t)dt is called Berry phase PTI and has a geometric mean- 
ing. On a closed path in Hilbert space the Berry phase always has the same value, 
regardless of whether one traverses the path fast or slow, since Berry's phase is 
invariant under reparametrizations of time. 
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6.5.1 Application of the adiabatic theorem to dark states 

hw=£i+>-£i-> < 6 - 83 > 

with 



n = V|fi+| 2 + |^-| 2 (6.84) 

Let us consider time-dependent Rabi frequencies which are slowly varying, VL± = 
Q± (t) . The dark state then also becomes a function of time, 

H'd) = H'd(T)) . (6.85) 

We assume that the initial state of the atoms is the initial dark state, \ij)(0)) = 
|^d(0)). Since the dark state is an eigenstate of the initial Hamiltonian we then 
find that the state follows adiabatically the evolution of the dark state \ip(t)) = 
|^d(T)). This is the principle of STIRAP. With this technique it is possible to 
transfer atoms from one ground state to another with very high efficiency. More 
specifically, one does the following. 



-e- 





+> +> 



Figure 6.8: The initial state is on the left side, the final state on the right side. The 
circle indicates the dark state. 



. |V(0)> = |-> = with n_(o) = 

• Change Q + and f2_ such that fi+(ti) = 

• Then = 1+) = IMh)) withfi+fa) = 0. 

This procedure corresponds to a counter-intuitive sequence of two laser pulses 
where one first applies on that laser pulse which does not couple to the atoms. 

In practise, STIRAP has a maximum efficiency of up to 99%. It is remarkable 
that, from the perspective of the adiabatic theorem, the atoms always remain in 
a superposition of |— ) and |+) (\4>d)) even though these states are only coupled 
through |e). In the ideal case (infinitely slow variation) the transfer from | — ) to 
|+) happens without ever exciting the atoms. However, one then also has to wait 
for an infinitely long time. 
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Figure 6.9: The first pulse corresponds to the second to fi_ 
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Chapter 7 

Incoherent Interaction and Density 
Matrix 

7.1 Density Matrix Formalism 

The density matrix g allows us to describe the time evolution of incoherent sys- 
tems (opposed to pure quantum states). In particular, it also makes it possible to 
describe open quantum systems. 

An example: A (classical) apparatus in a laboratory produces a quantum state, but 
only with a certain classical probability. Each time the experiment is performed, 
the apparatus may produce a different state with a certain probability. This could 
for instance be an oven which produces an atom beam sometimes with velocity Vi 
and sometimes with v 2 . 

One therefore only knows with some probability i = 1,2,... whether a state 
has been prepared. The states j^j) need not to be orthogonal. Of course the 
classical probabilities must fullfill J2i'Pi = 1' smce some state has been prepared 
with certainty. In a closed quantum system (not coupled to a reservoir), each state 
\ipi) then evolves unitarily according to \ipi(t)) = e~% Ht |^j(0)). If in one run 
of the experiment the state \ipi) has been prepared, the expectation value of any 
observable for this state is given by 

(JfyHViWWiW) (7-1) 

In each run of the experiment, the state j^j) is prepared with probability p^. One 
can calculate the full expectation value by averaging over all states which can be 
prepared, 

(A) = J> (^h = J> (^(t)|a|^(t)) (7.2) 

i i 
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If {|n)} is an orthonormal basis of Hilbert space, one can write the identity oper- 
ator 1 as J2 n \ n )( n \- We therefore find for the expectation value 

in n i n 

(7.3) 

The Density matrix £ is defined by 



(7.4) 



The expectation value of an operator then can be expressed as a trace over the 
operator and the density matrix, 



(A) = Tr (U) 



(7.5) 



Properties of the trace 



(i) The trace does not depend on the choice of basis. 



(ii) Tr (AB) = Tr (BA) 



Proof: 



ad (i) Let {|n')} be another basis of Hilbert space. Then 

Tt(A) = (n\A\n) = H Yl l n '> ( n '\^\ n ) 

n n n' 

= y ( n '\ i n )^ n i n '> = Yl ( n '\A\n') 



ad (ii) 



Tr (AB) = Y (n\AB\n) = ^ (n\A\n) (n'\B\n) 

n n,n' 

= (n'\B\n) (n\A\n) = (n'\BA\n') 



7.1 Density Matrix Formalism 
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Properties of the density matrix 

(i) £ is Hermitean. 

(ii) Tr^= 1 

(iii) Tr (g 2 ) < 1 (= 1 for pure states) 
Proof: 

ad (i) trivial because of E [0, 1] 
ad (ii) 

Trf? = ( n b#i) (^l n ) = Y Pi \n) (n\ipi) 

n,i i n 

= Y Pi (V'ilV'i) = Y Pi = 1 
i j i 

ad (iii) 

Tr(^ 2 ) = Y ( n \ Y Pi ^ WM^i) ^M) = Y PiP i i^il^i) (V'il^i) 

n ij i,j 

= Y PiP i I (V'ilV'i) I 2 < = 1 



For a pwre state one has: 



<i 



A mixture has at least two different states: 

Q = Pl +P2 \lp2){lp2\ + • • • , Pi 7^ 

The physical difference between a mixture of two states and a superposition 
becomes very clear by considering the following example: Let qg describe a mix- 
ture: 

Qc = \i\e){e\ + \g){g\) 

Half of the atoms are excited, and half are in the ground state. Let us calculate the 
expectation value of the dipole operator, d oc \e)(g\ + \g)(e\. 

(d) oc Tr(M\e)(g\ + \g)(e\)) = ^Tr(|y)<e| + \e)(g\) 
= 1 -(e\...\e)+ 1 -(g\...\g} = 
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Let q s describe the following superposition, 

|V) = ^(|e> + |<7» 
Qs=\^)^\ = \{\e) + \g)){(g\ + \e)) 

= \{\9)(g\ + \e){e\)+\{\e){g\ + \g)(e\) 
v v ' 

QG 

q s has non-diagonal matrix elements, which are called coherences. The expecta- 
tion value of the dipole operator is given by 

(d) cx Tv{U\e)(g\ + \g)(e\)) = Tv{{g G + l -{\e)(g\ + \g) (e\)){\e) (g\ + \g)(e\)) 

= Q + l Tx{\e){e\ + \g){g\) = l 

2 * v ' 

i 

The non-diagonal elements of the density matrix describe the degree of quantum 
mechanical coherence of the system, i.e., its ability to show intertference effects. 
This is the reason why they are called coherences. The diagonal elements are 
equaivalent to the population of each state. 



7.2 Liouville Equation, Superoperators 

We consider the time evolution of g in a closed quantum system. A pure state of 
course fullfills the Schrodinger equation: 

This can be used to derive an equation for the density matrix: 

ihd t g = ih^2pid t (\ipi)(ipi\) = ih^pi(d t \ipi) (ipi\ + |^>i) d t (^|) 

i i 
= J2P^( H \^P^)(^\-\^)(^Pi\H) 



ihd t g = [H, g] (7.6) 



This is the Liouville equation for closed systems. 
The trace of the density matrix is time independent: 



ihd t Trg = Tr([H, g}) = Tr (Hg - gH) = 



7.2 Liouville Equation, Superoperators 
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[H, ■] is a special cases of so-called superoperators, which act on operators instead 
of states. 

Cg := [H, g) 

is called Liouville" =Superoperator. Let us now solve the Liouville equation 
ihd t g = Cg. Formally this is quickly done: 

g{t) = e"^(0) = (1 - % -tC - ^CC + ■■■ )g(0) 

= Q{0) - f g(0)} ~ ^f- 2 [H, [H, g(0)]] + ■■■ = e-% m g(0)ei m 

Proof: 

f{t) := e^ c g(0) - iW t f = Cf = [H, f] 
git) : = e-i Ht g(0)e* Ht 
to show: fit) = git) 

ihd t g = ihi- l -Hg + l -gH) = [H,g] 

fit) and git) fullfill the same differential equation. 

fiO) = giO) = g(0) 

fit) and git) fullfill the same initial conditions. 

= g(t) 

Entropy, Equilibrium Density Matrix (This is a short deviation to statistical 
mechanics) 

The entropy of a quantum system is defined by 

S = -k B Tr (gin g) = -k B (In g) (7.7) 
S describes the degree of order in a system. 

In ( 17.71 ) we take the logarithm of the density matrix. Generally, a function fiA) 
of an operator is defined by refering to its diagonal representation. 

A \a n ) = a n \a n ) — > A = a n \a n )(a n \ 

For Hermitean operators we have (a n \a m ) = 5 nm 
f should be expandable into a Taylor series, 

i 
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This can be generalized to arbitrary powers: 

A = \a n )(a n \ 



i\ 

l n 



f(A) = ^ if ^2 a n \a n )(a n \ = ^ /(fln) K)(a n | (7.8) 



Z! 

We can employ this to calculate the entropy of a pure state: 



Q= l-\l/>) (lf>\ 

S = -k B Tr(g\ng) = -k B Tv(l ■ \t/,)(if,\ • ln(l) \if,)(if,\) 

=o 

The entropy of a pure state is 0. 

For a system in equilibrium we have 



H,g 







Q = <£4> 

=>- g has the same eigenvectors as H 

H \E n ) = E n \E n ) 

=> Q = ^Qn \E n )(E n \ 
Thermal equilibrium: maximise the entropy S under the constraints 
• Tr (g) = 1 (Conservation of total probability ) 

H) = E (For canonical or grand canonical ensemble) 



To do so we introduce two Lagrange parameters Ai, A2. Maximising results in 





dg n 

Since g is diagonal in this basis we find 

_d 

dg 



[s - Ax(Tr (g) - 1) - A 2 ((#) - e) = 



- I -k B ^2 Qn> ln Qn> ~ ^ \ ^2 Qn ' ~ 1 ) ~~ ^ 2 ( e n ' En ' ~ E ) f 

n I V n' / \ n' /J 



= (In g n + 1) - Ai - \ 2 E n 
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g n = exp j-^T" (^1 + ^K) | 

~ g K B 

Ai is fixed by ]T n £n = 1 

n 

Consider this as a function of H 

-p-H 
g ~ e fc s 

Let A 2 = J =► e = ie~^, /3 = ^, Z = Tr (e"^) 

Liouville's Theorem The entropy of a system is not changed if a unitary trans- 
formation is applied to the density matrix operator. 

g' = UgU ] 

S' = -k B Tr (g'lng') = S(g') = S {UgU ] ) 
Taylor expansion of the operator inside the trace: 

K(g') := g'lng' = Y, K nQ' n = Y. K ™ (WT 

n n 

Evaluate the term in parentheses, 

(UgUt) n = UglfU^gU^UgU^ . . . 
i 

= Ug n U ] 
^K(g') = uY,K n g n rf 

-k B Tr (UglngU^) 
= -k B Tr (U j Ug\ng) = S 

The entropy cannot be changed by the time evolution U = exp(—iHt/h) 
in a closed system. If one wants to change it, the quantum system must be open. 
This means that the system can exchange energy, momentum, or particles by a 
coupling to the environment or reservoir. This is exploited in cooling of quantum 
systems. 
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7.3 Reduced density matrix, Zwanzig's Master Equa- 
tion 

Direct product of Hilbert spaces Consider two Hilbert spaces Hs, Hr with 
bases {(S^)}, {\R a )}. Their direct product is given by 



^2*Pna \3 n ) g> \R a ) 



This space does not only contain product states of the form \S n ) <S> \R a ), but also 
entangled states, z.B. 

|Si> ® \Ri) + \S 2 ) ® \R 2 ) 
Entangled states play a very important role in many modern areas of Physics 
(Quantum information, Bell inequalities, Quantum computing, ..)• They are spe- 
cial with respect to the correlations between the two subsystems R and S. Cor- 
relations in entangled states cannot be explained by any cassical model. More on 
this topic can be found in Ref. [ 19 1. 

Let q be the density matrix in the space Hs <8> Hr. The expectation value of 

operators which only act on states in Hs is given by (^As^ — Tr (g (^A s (g> 
Here, the operator As acts only on Hs'- 

A s :Hs^H S :\S n )^A s \S n ) 

As CS> II r is the representation of this operator on the product Hilbert space. It acts 
like 

A s ®l R :H s ®H R ^Hs®H R 

\S n ) ® \Ra) A S \S n ) ® l R \R a ) 

7.3.1 Reduced Density Matrix, Projection Superoperators 

The trace of an operator in product space can therefore be decomposed as 

Tr {6) Hs ®h r = (^"l ® & (l 5 »> ® 1^)) 



n ci 

= ^<S , n |Tr*(£)A s |S' n ) 



7.3 Reduced density matrix, Zwanzig's Master Equation 
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In the last step we have used that the operator acts only on one subspace. 

Tr R (q) EE ^ ( R a\e\ R a) 



Define 

Qs = Tr i? (Q) 



reduced density matrix. 

q s act on Hs ■ In particular we have 

A s ) = Tr 5 (q s A s ) = Tr 5 (Tv R (q)A s )) 

Obviously, the reduced density matrix contains the full information about the ex- 
pectation values of all observables on the system Hilbert space Tis- 
Aim: Find an equation for the evolution of q s alone. 

Projection Superoperators In quantum mechanics, ordinary projection opera- 
tors P, like P = |e)(e| with P \ip) = \e) (e\ip), fullfill P 2 = P, pt = p 

Analogously we find for the action of a projection superoperator on an opera- 
tor 

V 2 A = VA 

Examples: 

1 . Reduced density matrix 

Vq = Q R (0)Tr R (q) , (7.9) 
where Qr(0) corresponds to the density matrix of the reservoir at time t = 

V 2 Q = M0)Tr fl \V q) 

= Qr{0)Tv r {q r {0) Tv r {q) ) 

V 2 q = qrWYtr (q r (0)) qs = Vq 

2. Projection on the populations 

Vq ee |n) (n| q |n) (n\ 

n 

V projects q on its diagonal elements, i.e., all coherences are removed. 
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7.3.2 Zwanzig's Master Equation 

Zwanzig's Master equation is an equation for Vg alone, which is derived from 
Liouville's equation 

ihg = Cg 

for g. To do so we first introduce a second projection superoperator Q := 1 — V, 
which fullfills Q 2 = The projection superoperator is time-independent, V = 

0, 

=> ih^-Vg = VCg = VC(V + Q)g 
ot 

d 

=> ih—Vg = VCVg + VCQg 
ot 

analogously with Qg 

th^-Qg= QCVg+QCQg (7.10) 
ot 

First we solve ihQg = QCQg. When the superoperators are time-independent, 
the formal solution is easyly given: 

Qg(t) = exp(-iQ£Qt/h)Qg(0) . 

The general case is considerably more complicated and leads to a so-called time- 
ordered exponential. Since this appears very often in the interaction picture and 
in quantum field theory, we willconsider this in more detail here. 

The Time- Ordered Exponential 

Formal solution: 

i% f Qgdt= [ QC(t')Qg(t')dt' (7.11) 
Jo Jo 

=> Qgit) = Qg(0) -'-J QC(t')QQg(t')dt' (7.12) 

where in the last step we used Q 2 = Q. This solution is only formal since it 
contains itself: Qg(t) depends on Qg(t') with t! < t. If Qg(t') is rewritten using 
again Eq. ( 17.121) . one arrives at 

Qg(t) = Qg(0) - % - J* dt'(Q£(t')Q){Q e (0) - % - J* dt" QC(t")QQg(t")} 

(7.13) 



There is an important difference between projection operators and projection superoperators. 



/ P lib) 

For projection operators one has \ip) = I ^ ^ ) . However, for superoperators we have A 

f VAV vAq 
{ QAV QAQ 
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Iteration leads to 



=> QQ(t) = Qg(0) - % -J^ dt'(QC(t')Q)Qg(0) 

+ l dt ' f df(Q£Q(t'))(QjCQ(t"))Qg(0) + ... (7.14) 

We introduce the time-ordered product : 

(7.15) 



T(A(t)B(t')) :-- 



A(t)B(t') for t > t' 
B(t')A(t) for t' > t 



This allows us to rewrite the multiple integrals in the series in the following way: 

dhQCQih) / dt 2 QCQ(t 2 ) = dt 2 dt 1 T(Q£Q(t 1 ) QCQ(t 2 )) 
Jo Jo J to. 



(7.16) 

T(QCQ(t 1 ) QCQ(t 2 )) has also the correct ordering of the operators (the operator 
at a later time is to the left of that at an earlier time) if t\ <t 2 . Consider therefore 
the following integral with modified boundary values for t 2 : 



dii / dt 2 T(Q£Q(h) QCQ(t 2 )) 
Jo 



dh / dhQCQih) QCQ(t 2 ) + / &bx / dt 2 Q£Q(t 2 ) QCQ(h) 
Jo Jo Jti 

(7.17) 

The first term on the right-hand side (rhs) of Eq. (17.171) appears in the series 
(Eq. ( 17.141) ) vor. The second term on the rhs of Seite von Eq. (17.171) can be written 

as 

/ dh / dt 2 QLQ(t 2 ) QCQ(h) = dt 2 dhQCQih) QCQ(h) . 

Jo Jt x Jo Jo 

(7.18) 

With t[ = t 2 und t' 2 = U it follows that 

2ndTerm in Eq. (|7~T71) = [ dt[ [ * dt' 2 QCQ{t[) QCQ{t' 2 ) (7.19) 

Jo Jo 

Eq. ( 17.191 ) is again a term which appears in the series Reihe (Eq. (I7.14I) ). 

t t-ti 

dti / dhQCQih) QCQ(t 2 ) 



1 fdh f dhTiQCQih) QCQ(t 2 )) (7.20) 
^ ! Jo Jo 
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We thus have found an expression for the first term in the series in which the 
integrals have independent boundaries. It can be shown that generally 

rt rti rtn-i 

dt, dt 2 ... dt n Q£Q{t 1 )...QCQ{t n ) 
Jo Jo Jo 

= - I dh... [ dt n T(QCQ(t 1 )...QCQ(t n )) (7.21) 
ni Jo Jo 

hJo 

+^ f l dtl - J dt nT(QCQ(h) . . . QCQ(t n )}Qg(0) 

= Texpf^-J dt'Q£Q(t')]Qg(0) (7.22) 

= :U QQ (t)Q e (0) (7.23) 

This is the solution of the homgeneous differential equation 

ihQg = QCQ(Qg) (7.24) 

The time-ordered exponential T exp is defined by its Taylor series, in which in- 
tegrals over different variables ti,t 2 , ■ ■ ■ do appear. The time-ordering operator 
T guarantees that the operators appear in the same order as in the series (17.141) 
erscheinen. 

Derivation of Master equation continued 

To derive the Master equation it is now useful to make the ansatz Qg(t) = UQQ(t)g~Q(t) 

=>■ ihd t Qg = ihU QQ g Q + ihU QQ g Q (7.25) 
With Uqq = -j:QCQ(t)U QQ (t) and using Eq. CUB we find 

ih QQ = U Q l Q (t)QCVg(t) (7.26) 
= ~QQ^)- l j L J^t'U Q 1 Q {t')QCVg{t l ) (7.27) 

=>Qd(t) = Qq(0) - j.U QQ (t) £ dt'U^WQCQPQit') (7.28) 
Insert Eq. ( 17.281 ) into the equation for Vg: 



iKPg = VCVg + VC(t)Q^Qg(0) - l -U QQ (t) jf dt'U Q l Q (t')QCVg(t') 



(7.29) 
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Eq. (17.291 ) is Zwanzig's Master equation. 

This is a exact differential equation for Vq alone, but it generally cannot be solved 
without approximations. To simplify it we now consider the special case of a 
projetion superoperator of the form 17 .91 with the following features: 

1. g(0) = Qr(0) <8> f?s(0), i.e., there are no correlation between the system and 
the reservoir at time t — 

2- [qr(0), Hr] = 0, i.e., the reservoir is in a stationary state of Hr (we set 

H = Hr + H S + Him). 

3. R(g R (0)H int ) = e.g. H mt = —d ■ E, here the reservoir corresponds to 
the radiation field 

=> Tr fi (g R (0)H int ) = -d (E) (t = 0) = 0. | 
When these conditions are valid one finds 

Qq(0) = q(0)-Vq(0) 

= Mo) ® Qs(0) - Mo) ® Tr i? (to(°) ® ^(0)) 
= gji(0) (8) es(0) - qr(0) ® es(0) 

= (7.30) 

The following relations for an arbitrary operator X are left to the reader as an 
excercise: 

VCVX 
VCQX 
QCVX 

Furthermore, 

QCQX = (l-V)CQX 

= CQX-VCQX 

= C(l-V)X-VC mt X 

= CX-CVX -V&mX 

= [H s + H R + Hi*, X] - [H s + H R + Hi*, Qr(0)Tt r (X)] 

-Qr(0)Tt r ([H m ,X\) 
= [H s + H R + H inU X}-[H s , 6R {0)TtR(X)} 

-[H R , q r (0)]Tt r (X) - [H inu Qr(0)Ttr (X)} 

-QR(0)Tr R ([H mt ,X]) 

2 In a stationary state in free space one has (E) = 0, since such a state corresponds to a mixture 
(i.e., not a superposition) of number states. Because E ~ a + a + one then finds (n| a + a + |n) = 

(n| (V«|n - 1) + \AH 7 T|"+ 1» = 



= r(O)[Hs,Ttr(X)]=VCsX (7.31) 
= to(0)Tr^([if int ,X])=P£ int X (7.32) 
= [H iBt ,Q R (0)Tr R (X)] = £ iot VX (7.33) 
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Using [H R , Qr(0)] = one deduces 

QCQX = [H s , X-VX} + [H R , X] + [H Mt X-VX}- V[H mU X] 
= [H s , QX] + C R X + [H m , QX] - VC mt X 
= (C s Q + jC R + C int Q-VC im )X (7.34) 



Additional properties: 



C S Q = QC S (7.35) 
VC S = jC s V (7.36) 



Insertion into the master equation (Eq. (17.291) ) results in 

ihVg = VC s Vg - -VC int J dt'UQQ{t)UQ X Q {t')C int VQ{t') 

U QQ (t) = Texp (^1 j\t"QC{t")Q\ 

C(t)X = [H(t),X] = [H s + H R + H mt , X] 

Approximation: take only terms up to 2nd order in H int into acount (assumption: 

Hint <C H R , Hs) 

=> Uqq « exp [~jjQ(£s + C R )Qt 

(if H s and H R are time independent). 

QC S = C S Q QC S Q = C S Q= QC S 

QC R = C R Q = C R 

=> Uqq = e -iQ£st/h e -iC R t/h 

In the last step C S C R = C R C S was used, which can be proven in the following 
way: 

(C R C S - C S C R )X = [H R , [H s , X}} - [H s , [H R , X}} 
= [H R ,[Hs,X]] + [H s ,[X,H R ]] 
= -[X,[H R ,H S ]} =0 
=o 

where Jakobi's identity was employed. One then finds 

e -iQC s t/n = y^{ — l J_L(QC s y (QC 8 ) n = C n s Q n = C n s Q n>0 

Z — / n \ I 

n=0 
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_^ e -iQCst/hQ£<p _ e -i£ s t/HQjr<p 

so that the master equation can be cast into the form 

iKPq = VC s Vq - l -VC mt f dt' exp { ^ (t - t') {C s + C R ) } C int Vg(t') 
h Jo . ft 



free evolution 



Vq = q r (0)Tt r (q) = Q R (0)Qs(t) 
ihg R (0)g s = g R (0) Tr R ([H s , Qr(0)q,]) - 



[H s ,8s] 
-t 



- j:Qr(0)T*r ([H M , £ dt!e^^l £ ' +c *V h [H ilU , te(0)&(f)]] 



ihg.(t) = [H s ,g s (t)} - dt'TrR r[fr iBt , e - i ( t -* , )^ + ^/ R [fr int ,^(o)fo(0]] 



This is the general master equation for a system fulfilling H int <C H R , H s and 
H s = H R = 0. 

Assumption about the form of the interaction: H int = —diEi, i = 1, . . . , N 
where operators di act on H s and Ei on 7i R . This is a very general form which 
applies to most systems. 

Tv R ([H mt , ...])= Tr R ([diEi, e^^X^*)/*^,, ^(0)^')]] 



with e~ iCt/h X = e -ivt/h Xe im/h ^ see beginning of this chapter) 

=► e-^-^+^d^QRiO^it') 



with 

d t (t - t') = e -i(t-t')Hs/Hj ie i(t-t>)H s /H 

E^t - t') = e -i(t-t')H R /n Eie i(t-t')H R /H 

QR(t) = Qr(0) since [H R , q r (0)] = 
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Remark: di(—t) denotes the operator di in interaction picture 

Tr K ([H in t, . . .) = Tv R ( dtEi, [diit-t'^it-t'), g R {0) U s (t-t')g s (t')Ul(t-t'l ] ) 

= -Qs 

with f/" s (t) := e - iH ^/ h 

^Tr R (...)= Tr R (cZ;£^(t - - H)qr(0)q. ~ d^g^O^d^t - - t 7 ) 

- - t')£z(* " t')g R (0)g s diEi + Q R {0)g a di{t - t')Ei{t - 0<*i#») 

= ddiit - t')g s Tr R (EiEttt - t')g R (0)) - d^d^t - t')Tr R (E ieR (0)Ei(t - t'))- 

-dt(t- t')g s d{Yt R (Et(t - t)QR(0)Ei) + g s dt(t - t')diTr R (g^E^t - t)E % ) 

Define F mim2 (ti,t 2 ) := Tv R (E mi {t x )E m2 {t 2 )g R {Q)) = {E mi {t x )E m2 {t 2 )) 

Tr R (...) = did^t - t')g s F u (0, t-t')- d^d^t - t')F u (t - t' , 0)- 
-di(t- lf)Q a diF u (0, t-t') + g s d t (t - t')d t F H {t - t 1 , 0) 

ihd tQs (t) = [H s , g a (t)] -ij o dt ' 1 ~ K Mt - t')g s ] + 

F K (t-t',0)|> a d { (t 

We now employ the widely used Markov approximation. Idea: The reservoir's 
correlation function Fu(0, t — t') decays rapidly with t — t' . (The reservoir "has 
no memory" for its own state). For small t — t' one then can approximate g s (t') 
by 

g s (t')^U s (t'-t)g s (t)Ul(t'-t) 

U s describes free propagation. One thus assumes that the system evolves freely 
for short a time. 

Qs = U s (t- t')g s (t')Ul(t - t') = g s (t) 

ihd tQs (t)= [H s ,g s {t)} - dt , (F il (0,t-t')[d i ,d l (t-t')g s (t)] + 

F li (t-t',0)[d l (t-t')g s (t),d i \ 
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Define r = t-t' dr = -dt' 

^lhd tQs (t) = [H S ,Qs(t)]~J^ dT^F il (0,T)[d i ,d l (T)Q s (t)]+F ll (T,0)[Q s (t)d l (T),d i } 

Since Fu(0, r) decays rapidly one can extend the integration to oo 
=>• Final result for the master equation in Markov approximation: 



ihdtQ.it) = [H a ,Qs{tj\ - ^J dr^F a (0,r)[di,^(r)^(0] + F u (r, 0) [ Ss {t)di{r), dj 



The master equation in Markov approximation 

• is a linear differential equation and, in contrast to the general master equa- 
tion, no integro differential equation. 

• has constant coefficients (given by the integral over r). 



7.4 Spontaneous Emission 

To describe spontaneous emission in free space we treat the atom as the system 
and the quantized radiation field as a reservoir. From the perspective of the atom, 
the state of the reservoir is essentially constant since spontaneously emitted pho- 
tons do not return to the atom. 

H s = ^ E n \n) {n\ , \n) — atomic energy levels 

n 

H R = / d 3 k^2huj k al^a k , a 

H int = -dE{x) = -d i+) E {+) - d^E^ 

where x denotes the position of the atom and we have made the rotating wave 
approximation. 

E^ + ' = positive-frequency part (oc a k a ) 
= creation operator 

• For a two-level system: d^ = d eg \e) (g\ 

• For two multiplets: = ]T \J e , m e ) (J g , m g \ d m ^ mg 
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The initial state of the radiation field is the vacuum: qr(0) = |0) (0|. In Markov 
approximation, the master equation the becomes 

ihg s =[H 3 , Q S ] --J dr j^, dj(r)Q S ] (E^t)) (7.37) 



Here we setH 

d a = d (+) ; E x = 

d 2 = d("> ; E 2 = E^ 

It is easy to see that for the vacuum state for the photons one has 

(0\EP4 +) \0) = (O^H^IO) = (0\EP4 +) \0) = , 
resulting in the master equation 

'[*,..-]<..■>■■■} = \(E { a +) Et\r))[di + \di-\r) Qs ] 



(E { b + \r)E^)[g s di + \r),di 



The correlation function has the following form (the derivation is left as an excer- 
cise for the reader) 

(0\Ei + \x,0)Et\x,r)\0) = TT ~^6 ab / du k u 3 k e'^ (7.38) 

With the aid of the free evolution operator U (r), the second correlation function, 
in which r appears in the first operator, can be expressed as 

(0\Ei + \r)Et\0)\0) = (0\U(r)Ei + \0)U\r)Et\0)\0) 
Usef/|0) = |0) : 

mi + \r)Et\0)\0) = (0\El + \0)U\T)Et\0)U(T)\0) 

= (0\E { a +) (0)Et\-r)\0) 

The time dependence of the atomic operators for free evolution is given by 

d {+ \r)= \Jem e )(J g m g \d metmg e^ E ^- E ^ h 



3 Strictly speaking each di in the sum above denotes a vector component of either or S \ 
To keep the notation concise we omit this here. 
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'o 



fc roo poo 

Y,( d m e ,m g ) b \J e ,m e )(J g ,m g \ u 3 du dr e -(-— e,« 9 ) 
„ ™ Jo Jo 



67r 2 eoc 3 



m e ,m 



i'"9 



witho; me)ma := {E me -E mg )/h. Using / °°(lre laJr = 7r5(u;) - zT/a; one arrives 



at 



/■OO »r 

/ dr 4 + )(r) (^ +) (r)^- } (0)) =— — £ | J e , m e )(J 9 , m 9 | (d. 



3 ■ / j ^ 3 



^m e mg) , 

^ ^ | t/ e , ?72 e ) ( J^, 777.^ | (d me , m9 )ft(i^ memg -\- iK m 



m e ,m g 

The operators ^ < " + ' ) and X ^ + only act on the atomic degrees of freedom and are 
time independent. Our result for the master equation for atomic multiplets then 
becomes 



ihgs = [H S ,Q S ] - \ {[d£\(K^-ik>t ) )es] + [Qs{K^ + ik'^),S; 



where we have introduced ^ := ( K^*) and K ' ' := ( K^ + ' 



(7.39) 

(***)'. 

For a two-level system these operators have the simple form 

j_ ^ POO 

K {+) = \e)(g\ d eg -p^; = \e)(g\ d eg ... (7.40) 

- '(+) 

The numerical value of the integral in K is irrelevant. It not only diverges (and 
therefore has to be renormalized), but it also leads to an incorrect value when it is 
calculated in two-level approximation and rotating wave approximation. Inserting 
these operators into Eq. (17.391) we find the Master equation for two-level atoms 

1 1 

q = __[^j | e }( e | _ -[^A Lamb |e)(e|,^] (7.41) 

1 

— — (a + (T_Q + g<j + a_ — 2<j_q<j + ) 
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with 
and 

A 

'Lamb 



\r1 I 2 , ,3 

7 := (7.42) 



K g | 2 r 00 Pa; 3 
Ai, mh := - _ " — t / au;- 



QTT 2 he c 3 J (lu — lu ) 



Consequences Obviously, the real part K^ + > and the imaginary part K /( - + > de- 
scribe very different phenomena. 

- '(+) 

• K contains an integral over lu with a principal value and is therefore 
related to the commutator oc A L . It describes a shift in the resonance fre- 
quency: lu' = lu + Ax. A L is called Lamb shift. 

• Ax is divergent and needs to be renormalized. A derivation which goes 
beyond the two-level approximation can be found in Ref. [8]. A^ does not 
change the unitarity of the time evolution. 

• The real part K + is the origin for the term proportional to 7 which de- 
scribes a non-unitary time evolution (see below); 

7 denotes the spontaneous emission rate . 

• The expression — ^(R^Rg + gR) R — 2RgR)) corresponds to the most gen- 
eral form of an incoherent time evolution which conserves the probability 
and is local in time. It is called Lindblad form. R is a generalized raising or 
lowering operator. 

Example: a thermal radiation reservoir at temperature T instead of |0) =>- 
thermal excitation of the atom: 

—r](T)(a^a + g + ga^a + — 2a + ga^) . 

More on this can be found in Ref. @. 

In the following we will assume that the resonance frequency ujq already con- 
tains the Lamb shift, ujq = cu' + A^. The master equation is given by 

% 7 
q = --[hu) \e)(e\ ,g\ - -{a + cr_g + ga+O- - 2a^ga + } (7.43) 

We write the density matrix as 

f £ee Qeg \ 
\Qge Qgg J 
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and first calculate the commutator with the system Hamiltonian: 

p, 



• Q 



[\e)(e\,g}= (I ° 

For the Lindblad form we find 

a +(T _ = \e)(g\g)(e\ = |e)(e| 







so that 



2^ee Qeg 
Qge ~2g 



\e)(e\ g + g\e)(e\ - 2\g)^g\e)(^-- 

The master equation then can be brought into the form 

Gee g eg \ _ / g eg \ _ 7 (2g ee g eg . (7,44) 



J 2\g ge -2q. 
which results in the following equations for the matrix components: 

^ee = ~igee => Qee{t) = £ ee (0)e~ 7 * (7.45) 

Qeg = (-^0 - ~) geg Qeg(t) = e -<Wot e _ 2 *p eg (0) 

g 9 g = igee g 9 g (*) = g 9 g (0) + (1 - e~ 7 * ) £ ee (0) 

This solution is characteristic for spontaneous emission: 

• The excited state decays incoherently since during the transition from e to g 
no coherences g eg , g ge do appear. In the case of coherent Rabi oscillations 
one finds that g eg , g ge first start to grow. 

• If there are initial coherences g eg , g ge they decay with rate |. This is because 
\g)(e\ contains an excited component which can decay. |e)(e|, on the other 
hand, contains two excited components. 

• The population lost in g ee is transferred to g gg , so that probability is con- 
served. 

Tr g = 1 -> Tr g = = g ee + g gg 

This is a general feature of the Lindblad form: 

Tr g = -1tt{I$Rq + gR ] R - 2RgR?} 

= -iTr {I$Rq + R?Rg - 2R ] Rg} = 
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• The master equation for the density matrix is superior to a phenomenolog- 
ical description of spontaneous emission based on the wavefunction only. 
With ip e (t) = e~ a t ip e (0) We find for the trace of g: 

\ij e \ 2 + \ij g \ 2 = e-^ e (0)| 2 + |^(0)| 2 ^ const. 

Probability is therefore not conserved. One may try to repair this by some 
ad hov ansatz, but such an approach leads to the creation of coherences, 
which should not occur in an incoherent process like spontaneous emission. 

Mt) oc ^(0)| 2 + (l-e-^)|^ e (0)| 2 g eg = W g ^ 

For t = and t — > oo the coherence g eg is going to zero, but during the 
evolution this is not the case. 

Visualization on the Bloch sphere: see Fig. 17. II 




Figure 7.1: 7T pulse (left) and spontaneous emission (right) on the Bloch sphere. 

Both for a Rabi oscillation with total angle n and for spontaneous emission the 
initial and final state are the same. The difference becomes evident on the Bloch 
sphere: for Rabi oscillations, the Bloch vector moves on a semi-circle its length 
remains unchanged. The system always remains in a pure state. During sponta- 
neous emission, the Bloch vector becomes smaller and evolves right through the 
center of the Bloch sphere, which corresponds to a maximally mixed state. 

Let us now consider the case when the initial state is an arbitrary superposition. 
Without spontaneous emission the Bloch vector evolves on a circle on the surface 
of the sphere (Fig. 17.21 In the presence of spontaneous emission, this is changed 
to a spiral downwards which ends in the ground state. Again the state is a mixed 
state for intermediate times (Fig. 17.21 right). To demonstrate the transition from a 
pure state to a superposition we calculate the trace of g 2 : 

Trg 2 = g\ e + 2g eg g ge + g 2 gg 
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Figure 7.2: Evolution of an arbitrary superposition for a Rabi oscillation (left) and 
for spontaneous emission (right). 

For g ee (0) = 1 we find from Eq. (17.451) the result 

Tr g 2 (t) = 1 - 2e" 7 * + 2e" 27 * (7.46) 
The time evolution is shown in Fig. Fig. 17.31 The absolute value of Tr g 2 quickly 




Figure 7.3: Tr g 2 during spontaneous emission 

drops to the value of \ (complete mixture) and then asymptotically returns to one 
again. This return to a pure state is only present if one neglects the atomic center- 
of-mass motion. Because of the momentum transfer associated with a sponta- 
neously emitted photon the final state of spontaneous emission corresponds to a 
mixture in momentum space. This process can roughly be described as 

-> J d 3 k\g,p+hk)(g,p + hk\ S(\k\-^-) 

(7.47) 
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7.5 Spontaneous emission cancellation through in- 
terference 

We consider a three-level atom in V-configuration: 



1+) 



\g) 

E e = E + = E_ 

In the general master equation (17.391) for atoms one then finds a term of the form 

poo 

/ drdl + \r)(0\E^(r)Et\r)\0) 
Jo 

\me){g\{dm eg )i{K meg +iK' m 



m e =± 



rn e g) 
& 



Since K„ leg , K' meg only depend on E e — E g they are equal for both transitions, 
H > g and > g. 

K (+) + iK K+) = {K + iK > ){ \ +) {g \ d+g + |_) {g \ d _ g} 

Inserting this into the master equation leads to 

ihQ s = [H s , g s ] 

-l{[d +g \+)(g\ + d_ g \-)(g\ , (K-iK')}(\g)(+\d* +g + \g)(-\d*_ g )Q s 
+(K + iK')[g s (\+){g\d +g + \-)(g\d^ g ) , d* +g \g){+\ + d% \g){-\ ] } 

Under the assumption that0 d + = d _ = d this can be reduced to 



ihg a = [H s , Qs] - j\d\ 2 \ (K - iK') (|+) 4 -)) (g\ , \g) ((+ f (-|) g s 



+(K + iK') ^(|+} + |-})<<?| , \g) ((+! + (-!) 



(7.48) 



4 The choice d + = e M d_ is also possible. 
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If the atoms initially are in the special excited state \ip) = ^(|+) — |— )) we 
find 

6a (P) = ty){rl;\, H s \4)=E e ty) - [H a ,Q.(0)] = Q 

In addition we have ((+| + (— |) \ip) = 0. Because of the term g s (\+) + |— )) = 
that appears in the commutators in Eq. (17.481 ) one then finds 

Q. = Q 

=>- \ip) can not decay into the ground state! 

• Destructive interference between the decay amplitudes of | +) and | — ) causes 
this superposition to be stable ifbH . 

• This effect reminds us of the dark state in a A-system. The major difference 
is that ) is an excited state and represents a dark state for all resonant 
photons. 

• The observation of this effect is difficult since the assumption d + = e ta d_ 
cannot be fullfilled for atoms. This is because atomic energy eigenstates are 
eigenstates of J z 

Example: three-level system in 
V configuration: 

J g = J e = 1 

One then has 

d + = V(e x — ie y )/ \pl couples to a + light 

d = V(e x + ie.y)j\pl couples to cr„ light 

Thus, d + ■ d* =0 and the interference terms in the master equation are 
zero. For instance, we have 

(\+)d + + |-)d_)^|.|^)((+|d; + (-|d:) = p 2 (|+)(+| + |->H). 

By contrast, under the assumption d + = <i_ used above one finds 

2> 2 (l+> + + ("I- 

• For molecules destructive interference is possible in principle. Experimen- 
tal observation of spontaneous emission cancellation was reported in Ref. |[T51l . 
but the result was later questioned ifToll . 
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Chapter 8 



Electromagnetically induced 
transparency 



EIT is one of many effects in quantum optics where light changes the state of the 
atoms, which in turn leads to a change of the index of refraction. The description 
of such effects is based on 



8.1 The Maxwell-Bloch equations 

Equation for atoms: g = Cg as in Chapter |7J with the interaction between light 
and atoms given by 

H = H A + H Rad -- [ P{x) ■ D(x)d 3 x 
e o J 

P(x) = Polarization; e.g. for a single atom at position x : (8-1) 
P(x) = d5(x — Xq) 

Heisenberg's equation of motion: (as in Chapter|2l): 



ihb{x) = [D(x),H] with 
[Di(x),Bj(x')] = -ih€ ijk d k S(x - x) 



from which follows 



(8.2) 



D = — curlB (8.3) 
B = — curl(D-P) (8.4) 
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These are Maxwell's equations U in a dielectric with polarization P. As we have 
learned in the first chapters, one can deduce wave equations from these relations: 

+ curl curl ^ D = curl curl P , (8.5) 

or, using D = e E + P, 

(J^d? + curl curl ^E= -fi P (8.6) 

Together with q = Cq (see Chapter |71) this equation represents the Maxwell-Bloch 
equations. These equations are generally hard to solve and one has to make ap- 
proximations. Frequently used assumptions are the following. 

• Homogeneously distributed atoms 

P(x,t) = P(t) (8.7) 



• The atomic density is very low so that the atoms do not interact directly 
with each other (collisions) and do only interact with the light. 

P(t) = gd (8.8) 

here q denotes the atomic density and d the dipole moment operator for a 
single atom. 

Averaging over the atomic degrees-of-freedom results in 

VoQd? (d) (8.9) 
^^Tr(^d) (8.10) 

The Maxwell-Bloch equation allow to calculate the optical properties of an 
atomic medium. In a gas one can neglect interatomic collisions to a certain degree. 

'We assume that the macroscopic free charge density vanishes, so that divD = 0. It is then 
not necessary to distinguish between the transverse part of the dielectric displacement field and 
the field itself. 



—d? + curl curl 



E 



8.2 Index of refraction for a two-level atom 
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8.2 Index of refraction for a two-level atom 

We consider a two-level atom in a running laser beam along the z-axis, 

E = £(z)ee ikLZ " iWLt + c.c. . (8.11) 

The laser beam is continuously switched on ("cw") so that the envelope £ (x) is 
time independent. For simplicity we assume that the width of the laser beam is 
much larger than the wavelength 2n/kLSO that the envelope depends on z only. If 
the variation of the envelope is small on the scale of the wavelength, we can make 
use of the paraxial approximation (see Sec. ll.8l) . The wave equation (18.91 ) is then 
reduced to 

2ik L e e lkLZ ~ luJLt d z £ - 2ik L e* e ~ lkLZ+luJLt d z £* = — — ^Tr(^Y8.12) 

c 2 e 

Our task is to find a suitable aproximation for the atomic density matrix g so that 
we can solve this equation. 



8.2.1 Stationary solution for a driven two-level atom 

The evolution of the atomic density matrix is governed by the master equation. 
For a two-level atom in free space it is given by (17.431) . In the presence of a laser 
beam, we have to replace the atomic Hamiltonian fkdo \e)(e\ by H = Jtlu |e) (e| — 
d ■ E. Proceeding as in Chapter^we find in the rotating-wave approximation the 
Hamiltonian ( 16.171) . i.e., 

H = huj |e)(e| - hVLe~ iu;Lt \e){g\ - hn*e iu}Lt \g)(e\ (8.13) 

where 

„_ d eg -EM(t = 0) d eg ■ e e^ z £ 

^ = h " h • (8 - 14) 



Similarly as in Sec. 16 .21 we can make a unitary transformation 

U{t) = exp{-iu L t |e)(e|) (8.15) 

With q = UqW we find 

H = -iW ] U + U ] HU 

= -hA\e)(e\-hn\e)(g\-m*\g)(e\ (8.16) 



with A = uj l — u , so that 

% 

H 



t = ~[H,g\- l{\e)(e\Q + Q\e)(e\-2\g)(e\g\e){g\} (8.17) 
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This is a set of linear differential equations with constant coefficients. If g is 
written as 

Q = Qee |e)(e| + Q eg \e)(g\ + g ge \g)(e\ + g gg \g)(g\ (8.18) 
the equations become 

Qee = ~lQee ~~ i(Qge^ ~ Qeg^*) 

Qeg = (l'A - |j f? efl + ifi(£> ee - f? 55 ) • (8-19) 

The equations for g ge and g gg are not needed because of g ge = g* eg and Tr (g) = 

gee g g g 1- 

To consider a cw laser field has a number of advantages. After a sufficiently 
large time (large compared to 7 _1 ) all transient excitations will be damped away 
by spontaneous emission. The atoms then settle in a stationary state in which co- 
herent (deexcitation and ) excitation by the laser field and incoherent deexcitation 
cancel each other. This stationary state is determined by setting g — 0, leading to 

~st I I 

6ee A 2 + ^ + 2|fi| 2 
p st = l - n st 

¥gg x t^ee 

Qea = ~ — ■ (8-20) 

A2 + £ + 2 |fi|2 



This is the solution in the "co-rotating frame". To get the solution in the "labo- 
ratory frame", we need to undo the unitary transformation U = exp {—iui,t | e) (e| ) . 
It follows that 

git) = U(t)~g st U + (t) 

= ~gt \e) (e| + ~g% \g) (g\ + ~g%e-^ \e) (g\ + ~g%e^ \g) (ej (8.21) 

8.2.2 Solution of the paraxial wave equation 

We now employ this solution to calculate the right-hand side of (18.121) by using 

Tr (g(t)d) = d eg g ge {t) + d* eg g eg (t) . (8.22) 
Inserting this into Maxwell's equation results in 

e e ikLZ ~ tuJLt d z S - e* e~ ikLZ+iu>Lt d z 8* = -%^- (d eg g st e e iuJLt + d* ~g%eX^l%) 
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Comparing the rapidly oscillating exponential factors (and polarizations) one thus 
finds 



Ze 



eg^eg 



-% a. 



2e e9 A 2 + 'f + 2\n\ 2 



From the Wigner-Eckart theorem (15.781) we know that 

d eg = (F e ,m e \d\F g ,m g ) 
= (F e ,m e \d q \F g ,m g )e* q 



(F g , m g ; 1, q\F e , m e )e* q ^Jii (8.25) 



MFe,F g ) 
V2F e + 1 

If the laser beam is either linearly or circularly H polarized, we can choose the 
quantization axis in such a way that one vector e qo coincides with the laser polar- 
ization e. Hence, 

d eg = Ve* (8.26) 



with V = (F g , m g ; 1, gol-^e, rn e )V(F e , F g ) / \/2F e + 1. The paraxial equation thus 
reduces to 

Q z g = -iM.\V\ 2 S (8.27) 

2e h l 1 A 2 + ^ + 2|0| 2 

In this chapter we make the additional assumption that the the factor of |f2| 2 in 
the denominator can be kept constant. One then arrives at a simple linear equation 
for E with the solution 

E{z) = expf-z^^P ^3 )S(0) (8.28) 

For the complete positive-frequency part of the electric field this implies 

= e exp I -iu L t + ik L z [ 1 - H-| A ~ %2 ] I £m.29) 

2e h A 2 + ^ + 2|fi | 2 ' 



2 This excludes elliptically polarized light, but it seems not be possible to construct a real two- 
level system for this case anyway. 
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It then becomes obvious that the solution of the Maxwell-Bloch equations cor- 
responds to a (complex) change in the wavevector of the electric field. We can 
compare this change with the definition of the index of refraction : 



k(u L ) = — n(u L ) 



leading to 



n{uj L ) 



2he A 2 + f + 2|ft | 2 



(8.30) 



(8.31) 



The index of refraction is generally complex, n 



n 



in 



E (+) 



E e 
E e 



tikx^-(n +in ) 

P c v ! 



■IU) L 



-iujj^t ik x — k x 



(8.32) 



The real part n determines the phase velocity v p h = ^j- of the light. For 
vanishing detuning the phase velocity becomes t>(A = 0) = c. 



Group velocity: Consider a wavepacket 

E(x,t) = [ Auoe ik{w)x e- iuJt E(uo) 



For a narrow frequency width w « w , we can make the expansion 

dk 



k(oo) ps k(u ) + 



du 



{u-Uq) + ... 



(8.33) 



(8.34) 



E(x,t) 



duE(uj)e- luJt e 



-iu)(t- f£| ^ x) 



o ' (8.35) 



We can infer that the wavepacket travels with a velocity determined by 



x fdkY 1 

7 = const= UjJ 



dui 
dk 



(8.36) 



This corresponds to the group velocity v gr = It varies particularly 

strong around the resonance A = 0. However, at the same time n" , which 
describes the absorption, is maximal for A = 0. One therefore cannot con- 
sider this situation as a real propagation of light since most of it is absorbed. 
The notion of a group velocity is then not very useful. 



8,3 Electromagnetically induced transparency and dark states 
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• "Spectral hole burning": |f2o| 2 is proportional to the light intensity. Hence, 
for high intensities absorption (and all other effects) is strongly suppressed 
even if A = 0. The reason is simply that a high intensity corresponds to 
many photons inside the medium. If all atoms have already been excited by 
a photon, the remaining photons cannot be absorbed anymore. 



8.3 Electromagnetically induced transparency and 
dark states 

The refraction index for atoms with more than two levels can be very different 
from the two-level result. One example is the phenomenon of electromagnetically 
induced transparency for a three-level atom in A configuration. To derive theindex 
of refraction we again have to solve the Maxwell-Bloch equations with 

-IQe 

h = n I a_ oi c sp . Em .g= I -} e . 





( o 


-q_ 








-n* 


A_ 





•^sp. Em. Q 









aw 





2' 



\&e- 


2 Qe- 


7 

2 Qee 








2 Qee 



(8.37) 



Master equation: 



Z ! 

q = — — [H,q] + £ sp . Em.Q = stationary solution (8.38) 

To find the solution of this system of linear equations is simple in principle, but to 
reduce the size of the result we consider the case that f2 + is small and expand the 
solution around f] + = 0. 

For f2 + = 0, the soltution is simply the dark state , i.e., 

Qo = !+}(+! (8-39) 

Since |+) is the dark state for the situation f2 + = and f2_ ^ 0, the entire 
population is pumped into the state |+). 

If we switch on a weak second laser beam (Q + ^ 0) we can calculate its effect 
using perturbation theory: H = H + Hi with 











ft* 


A_ 











A_ 




H Q = hi A_ Hi = h (8.40) 



Inserting the ansatz Q = Qo + Qi + into (18.381 ) we find 

i 

= -j\ H ® + H li Q0 + 6l] + £sp.Em.Q0 + £sp.Em.Ql (8.41) 
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Since g is the solution of the unperturbed equation, we find for the equation to 
first order in Hi 

% % 
- -^[Ho, Qi] + C-sp.Em.Qi = Qo] (8.42) 

This is an inhomogeneous linear system of equations for qi. The solution is given 
by 

e=l+)(+l + W|e><+| + e-+|-)(+|+H.c.} (8.43) 

QO 



with 



Qe- 

Q— 



(A_ - A + )fi H 



(A_- A + )(A + + if) + |a 
(A_- A+)(A + + a) + |fi_| : 



We now calculate the index of refraction for the f2 + beam with the aid of Eq. (I8.9I ). 
We know from the derivation of the two-level refraction index in Sec. 18 .21 that Pi 

" + ~ i+ l^ Tr < aH - e » + <4 + <8 - 44) 

For the Q + beam we have • £q + oc |+) (e|, so that 

„ + « i-l^l A -~ A + (845) 

+ 2fe (A_-A + )(A + + ^) + |ft_| 2 1 ^ 

The prefactor gV 2 / (2he ) is typical for optically thin gases and also appears in 
the two-level result. The term after the prefactor is responsible for EIT. 

For A_ = A + we have n = 1, i.e., the same refraction index as in free space; 
in particular, there is no absorption and the phase velocity is c. However, the 
system does not completely behave as a vacuum since the right-hand side is trivial 
only exactly at resonance. The group velocity ^ is different from its vacuum 
value. 

Let us consider the case A_ = (pump laser in resonance): 

n + = 1 ~ \7 , —\ I in 12 (8 - 46) 

2ne {ujq - oj + ){u) + - u + i\) + 



We have for the wave vector 



k + (cu+) = —n + {u + ) (8.47) 

c 



3 The derivation of this equation is completely analogous to that of Eq. 08.31b . 
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It is then easy to calculate the group velocity, 



du 



dk 



+ 



+ 



n + u + dn 



c c duj 



+ 



+ 



1 | u q qV 2 1 
c c 2he \Q- 1 2 



(8.48) 



c 



(8.49) 



1 + 



gp 2 up 

2he |0_| 2 



In principle this would allow us to achieve y gr = (for f2„ — > 0). However, in this 
limit the expansion (f2 + Q_) we have made is not valid anymore. Nevertheless, 
this limit has been achieved experimentally. 



Example uj = 10 15 s~\ V = ea , q w 10 21 m" 3 , w 10 8 s' 1 . This results in 
a group velocity of 



The first experiment by the group of Lene Hau lETTI achieved v gr = 17m/s in 
1999. A year later the group of Philipps sogar was able to stop light completely 
P4l . Stopping of light can be understood as a transfer of atoms from |+) to |— ) 
followed by switching off the pump laser. The information about the probe pulse 
is then coherently and reversibly stored in the medium. The fi + pulse can be 
restored after up to a millisecond ll35l . 



v, 




(8.50) 
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Chapter 9 
Photonic band gaps 



Photonic band gaps (PBG) are completely analogous to the energy bands of elec- 
trons in a periodic potential which are well-known in solid state physics. Their 
physical meaning is that in a periodic dielectric photons cannot propagate at cer- 
tain frequencies (inside a band gap). One can create such band gaps by a periodic 
modulation of the index of refraction. To simplify the notation we introduce the 
relative dielectric constant e r , 

si X) 

e r (x) = = n 2 (x) n: index of refraction (9.1) 

which we assume to be time-independent (e r oc e(x) = 0). We then have 

D = e E + P = e(x)E = e e r (x)E(x) (9.2) 
and Maxwell's equations become 

B = -curl E D = 5 c 2 curl B (9.3) 

which leads to a "symmetric" wave equation for B: 

(a 2 + c 2 curl ^rcurl ] B(x, t) = (9.4) 

V e r{X) J 

Consider now a periodic index of refraction (or dielectric factor) e r (x) = e r (x + 
L), which for simplicity only depends on xQ. To emphasize the analogy between 
photonic and elektronic band gaps we will treat the wave equation with concepts 
borrowed from quantum mechanics. 



'in the following we will nevertheless use the notation e r (x) for a dependence on all compo- 
nents if the equation is generally valid. 



180 



Photonic band gaps 



We first define a conserved scalar product (see also Sec. 12.1 .61) . which will 
allow us to introduce Hermitean operators: 

(B\B') := -i§ / d 3 x (B*B' - B*b) (9.5) 



d t (B\B') oc / d 3 x (B B' - B*B ) (9.6) 



inserting this into the wave equation (19.41 ) 

d 3 x 



d 3 x 



(-c 2 curl -^curl B*)B' - B*(-c 2 curl -^curl B') 

e r (x) s r (x) 

(9.7) 

B*(-c 2 curl — ^curl B') - £*(-c 2 curl — ^curl B') 

(9.8) 



In the last step we have performed two partial integrations to move the curl 
operator to the primed part. We also used that e r (x) does not depend on 
time and appears between the two curl operators. 

• The "Hamiltonian" operator H := —c 2 cm\e r (x)~ 1 cm\ is Hermitean, i.e., 
(B'\HB) = (HB'\B). In addition, H is periodic, H(x + L) = H(x). 

To derive the eigenmodes of the wave equation (19.41 ) we prooced analogously 
to solid state theory and first consider the Translation operator translation operator 

T xo = e'^ ^ = e XlyVx (9.9) 



Effect of operator: (T X J) (x) = f(x + x ) (9. 10) 

Proof: T x J(x) = e x °^°>f(x) 

oo 



_ sr^ — o — x — j., n _ jayiQj-.Rei^e um x 
^— ' n\ 

n=0 

= f(x + x ) (9.11) 

H commutes with T L , 

(HT L - T L H)f(x) = Hf(x + L) - H(x + L)f(x + L) = , 
because H is periodic. Hence, H and T L share a common system of eigenstates. 
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The eigenstates of Tl are the momentum states e %kx with eigenvalues T^e^ = 
e tk(x+L) _ e tkL e ikx_ jfols, means the eigenvalues are degenerate for all momenta 
of the form k = q + n E Z with eigenvalue exp(iqL). 
The eigenstates of H can be written as a superposition of degenerate eigenstates 
of Tl, 



OO 

; 2/] 7T 



Tl — — OO 



The quantum number q can be chosen to be in the interval [— f , f ] y. ^ is called 
quasi momentum (= momentum modulo 2ix/L). 

Eigenvalue problem: (j^ — H J tp = 0, Ansatz ij;(x, t) = e ±luJt ip(x) =^ 

— u) 2 ip(x) = Hip(x) 

Here we consider the most simple case: 

e r = 1 + oecos — x 



with Se -C 1 



2 1 

H = —c curl — r curl 



1 +5ecos (^x) 
I 



-c 2 curl ( 1 — 5e cos ( — x \ \ curl 



H + SeH, 



with 



'2tt 

i^i = c z curl cos I —x I curl 



i/ — — c curl curl 

i 

~L 

Consider a light wave along the x axis, 

B(sc) = B(x)e y 

This ansatz fulfills the Maxwell equation divB = 0. The action of the Hamilto- 
nian becomes 



1 \ / d 1 d 

-curl —-curl B = r^TT 3 ' e 

e r (x) / \aie r (x) ax 



2 Another choice can be made by changing the sumation index. For instance, for q' = q + 2w/L 
we get the same expansion as above if ip q>n is replaced by tp q ' t n-i- 
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We now employ perturbation theory, B(x) = (x\B), \B) = \B ) + Se \Bi), where 
\B ) is a momentum eigenstate p . 

{H + 8eH x ) (\B ) + 5e {B,}) = -J 1 (\B Q ) + 6e {B^) with oo 2 = uj* + deooj 



Equation for zeroth order: 



Ho \p ) = -Uq \p ) 



cpo 
% 



Equation for 1st order: 

H \B X ) + H x \B ) = -uj\ \B Q ) - ul \B X ) 

To solve these equations is a bit cumbersome. The structure of the solution is as 
follows. The change of the eigenvalue can be derived from —u\ = (B \Hi\B ) . 
The state has the form 

\B X ) ~ H x | B ) 

/2tt \ . 
~ rot cos I —x J rot \p ) 



+ e"^) bo) 

I 27T 27T.. 

bo + y" "> + bo - -jj h ' 



(?p 2 

H b) = — ^- b) 
= w%(p) b) 

For p « 4p the component c|po — ^ I is nearly resonant to c|po | while po + ftfc 
has a much higher energy =>- effectively only couples p to p ~ ^ f° r Po ~ if 
[§ We therefore obtain an effective two-level system. If |p 1 is very different from 
y, coupling is negligible. The energy bands are shown in Fig. l9.2[ 

Consider the two-level system for p ~ -y. Denoting the matrix element 
(po| -f^i bo — ^) a s ft-i one obtains the equations 

\B) ^a\p )+(3\po-hk) 



3 or p and + Kk f° r Po ~ — 
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Figure 9.1: Energy-(quasi-)momentum relation for electromagnetic waves in free 
space. The dashed line corresponds to the energy E(p) = c\p\ as a function of 
momentum p, solid lines the energy as a function of quasi momentum q. In free 
space, this is merely a difference in the use of quantum numbers: Instead of one 
continuous variable p one uses the continuous quasi momentum q (which in free 
space corresponds to a finite momentum interval) and the discrete band index n. 




-k/2 



Figure 9.2: Energy quasi momentum for electromagnetic waves in a periodic di- 
electric. At the position in momentum space where the free states are resonant 
one obtains band gaps. 
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5eHi \p ) w Sehi \p - hk) 
5eHi \p - hk) ps 8eh\ \p ) 
and the eigenvalue problem 

H\B) = -lu 2 \B) 

H\P>\= ( ~^o(Po) Seh \ ( a \ 
1 ; V 8eh\ -co 2 (p - hk) J \ (3 J 

with eigenvalues 

ul(p ) +u 2 (po - hk) i r - 2 
2 ± 2 V (^olPo) - ^o(Po - ^)) + 45e 2 

where the "detuning" is given by uj 2 (po) — u 2 (po — hk) ~ in which 

the momentum p = 4j- + Ap is expressed relative to ^A: /2. Then Ap plays a 
similar role as the real detuning appearing in the two-level system of Sec. 16.31 
By varying Ap one obtains, in complete analogy to Landau-Zener transitions, an 
avoided crossing. This is the reason for the appearance of band gaps. 
Significance of photonic band gaps: 

• Light with frequency inside a band gap cannot propagate through the medium. 
This can be used for mirrors, glass fibres, filters, . . . 

• Theoretical predictions: spontaneous emission of atoms placed in the pe- 
riodic medium will be strongly modified ['23- However, photonic crystals 
with band gaps stretching over the entire angular range of photon momenta 
have only been produced recently ll23ll . 
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10.1 Green's functions 

Consider a general linear inhomogeneous differential equation (DEQ) of the form 



It is not important that t and x represent space and time, they can be replaced 
by any number of independent variables. A Green's function is a special solution 
with the inhomogeneity 5(x — x')5(t — t'). 



The Green's function serves two puposes: 

1 . If one knows G (t, x; t'x'), the solution of the general inhomogeneous DEQ 
is given by 



Proof: 

{L t>x F) (t, x) = L tjX J dx'dt'G(t, x; tf, x') ■ h(t', x') 
= J dx'dt'L t!X G(t, x; t', x') ■ h(t\ x') 

2. The Green's function can also be used to solve the initial value problem 



(L t , x F)(t,x) =h(t,x) . 



(L t ,jG)(t, x- 1', x') =5(x- x') 5(t - t') . 




(10.1) 



(L t , x F)(t,x) = 



F(t ,x)=F (x) . 



(10.2) 
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Suppose the Green's function fulfills at equal times 

G(t,x;t,x') = 6{x-x') . (10.3) 

Then 

F(t, x) = J G(t, x; t , x') ■ F (x')dx' (10.4) 



Proof: 



{L t>x F){t,x) = J dx'L t>x G{t,x;t ,x')-F {x') 

= J dx'5(x- x) 6(t - to)Fo(x') = OVt > t 



In addition, 



F{t ,x) = J dx'G(t ,x;t ,x') ■ F (x') 



dx 5 (x — x') F()(x') = Fnix) 



This concludes the proof. 

It is often not too hard to find some Green's function for a given problem, for 
t ^ t' and x ^ x' it simply fulfills the homogeneous equation. One therefore 
can often compose G by using different homogeneous solutions for t < t' and 
t > t'. As a rule of thumb, if the highest derivative that appears in L is of order 
n, the n — 1th derivative of G must be discontinuous in order to cancel the delta 
distribution inhomogeneity. 

There are many different Green's functions for a given problem which differ 
by their boundary conditions. Basically, the choice of homogeneous solutions 
away from the delta inhomogeneity is fixed by the boundary conditions. One often 
is especially interested in finding the retarded (or advanced) Green's function. In 
this case, the n — 1th time derivative is usually proportional to the step function 
9(t — t') (or 9(t' — t)). Bear in mind that the derivative of the step function is the 
delta distribution. 



10.1.1 Green's operators 

In a more formel way, Green's functions are actually the integral kernel of a linear 
operator that acts on a Hilbert space. This is not often appreciated, but it is implicit 
in its DEQ due to the appearance of the delta distribution. This implies that G 
itself is also a distribution, and this actually fits nicely to its usual applications. 



10.1 Green's functions 



187 



Considering G as a distribution means that we can only consider integrals like 

(6{t,t')f){x) := J dx'G{t,x;t',x')f{t',x') . (10.5) 

This equation defines a linear operator G which maps f(t',x') to f(t,x). This 
operator has some nice formal features. If Eq. (110.31) is valid then 



(G(t,t)f)(x) = I dx'G(t,x;t,x')f(t',x" 



dx'5(x-x')f(t,x') 
lf(t,x). (10.6) 



Hence, G(t,t) = l. 

Likewise, one can introduce the operator 



(m,x)--=J^m,^mtw). <>o.7> 

whichs leads to 

{L tiX Gh) (t, x) = L tjX J dx'dt'G(t, x; t' , x)h{t' , x) 

dx'dt' L t)X G{t, x\ t' ', x')h(t', x') 



h(t,x) (10.8) 



which means LG = 1. The Green's operator is therefore formally the inverse of 
the DEQ operator. However, one has to be careful with this interpretation. Every 
solution of the DEQ corresponds to an eigenfunction of L with eigenvalue zero. 
Hence, this operator is not invertible and the above interpretation is only formal. 

An example for Green's operators can be constructed for the Schrodinger 
equation. The equation 

(id t - H)G(t, t') = 5(t - t') (10.9) 

has the retarded solution 

G(t, t') = -i9(t - t')e-^-^ H (10. 10) 



188 



Appendices 



Proof: 

(idt-H)(-i)9(t-lf)e^-^ H = 5(t - t')e- i{t - t,)H 

-iH6(t - t')e- i{t - t,)H - H{-i)B{t - t')e- i{t - t,)H 

= 5(t-t') (10.11) 

The initial condition is G(t,t) = —il if we set 0(0) = 1. The corresponding 
Green's function can be derived from 

= (x\ J dx'dx"\x')(x'\G(t,t )\x")(x"\ij(t )) 
= J dx"(x\G{t,t )\x")(x"\4>(t )) 
= J dx"G(t,x;to,x")iJj(to,x') (10.12) 
where we have set G(t, x; t , x") = (x\G(t,t )\x"). 



10.2 Angular momentum algebra 

At the beginning a word of caution: in angular momentum algebra different con- 
ventions are used in different textbooks, which are distinguished by their choice 
of phases of the angular momentum eigenstates. We follow the convention of 
Ref. \4\, which also seems to be used in the Mathematica software package. Other 
good references on angular momentum in quantum physics are Refs. [2] und 01 . 

Angular momentum is defined as a vector L whose three components are op- 
erators fulfilling the commutation relations 

[L a , L b ] = ihe abc L c . (10.13) 

It can be shown that [L 2 , Li] = ist. Hence, one can find a common set of 
eigenstates of L 2 and one of the components, which usually is chosen to be L z : 

L%m) = H 2 l{l + l)\l,m) , L s \l,m) = hm\l,m) (10.14) 

Here, m runs from —I to I and I — 0, 1/2, 1, 

It is advantageous to introduce the operators L± :— L x ± iL y , which obey 
[L±, L z ] = ^hL± and [L + , L_] = 2hL z . It follows that 

^hL±\l,m) = [L±,L z ]\l,m) 

= mhL±\l,m) — L z L±\l,m) (10.15) 
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Hence, L±\l,m) is proportional to \l, m ± 1). Since L\_ = X_ and L 2 
L-L + + hL z we find 

(l,m\L 2 \l,m) = h 2 l(l + l) 

= (l,m\L 2 z + L_L + + hL z \l,m) 



h l m l + h 2 m+\(l,m + l\L + \l,m)\ z (10.16) 



and thus 



(/', m'\L + \l, m) 

The phase factors of the states \l,m) can be chosen in such a way that the matrix 
elements are real. We have introduced the notation 

km ■= y/l(l + l)-m(m + l) (10.18) 

since this factor appears frequently. A very useful consequence of Eq. (110.171) 
is that for a given I all the states \l, m) can be obtained from \l, —I) by repeated 
application of the operator J + . 



= (l,m\L^\l',m') = hfi jTn 5 hV 5 m+ i,m' (10.17) 



10.2.1 Addition of two angular momenta 

When adding two (commuting) angular momentum operators L and S to J = 
L + S, one has to address the question of choosing a suitable basis for the product 
of the respective Hilbert spaces. A possible choice would obviously be the tensor 
product of the individual bases, 

\l, mf, s, m s ) := |/, mi) ® \s, m s ) . (10.19) 

However, it is usually more convenient to work with eigenstates of J 2 and J z 
instead. Because [L 2 , J,] = [S 2 , Jj] = these states are characterized by four 
quantum numbers 

b\r%,M) • (10.20) 

The matrix elements (/, mf, s, m s \j, nij, (I, s)) =: (I, mi; s, m s \j, rrij) are called 
Clebsch-Gordan coefficients and allow to switch from one basis to the other: 

\j, mj, (I, a)) = ^ {I, my, s, m s \j, m h ) \l,mi;s,m s ) (10.21) 

mi,m s 

By applying J z one can easily see that m s + mi = mj is valid. The numerical 
values for Clebsch-Gordan coefficients can be derived from the linear system of 
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equations 

(l,mi;s,m s \J + \j,mj,(l,s)) = hfj, mj (l,mi; s,m s \j,mj + 1) (10.22) 

= (I, mi', s,m s \(L + + S + )\j,mj, (l,s)) 

= (j,mj, (/, s)|(£_ + S-)\l,mi; s,m s )* 

= hfi tmi -i(l,mi - l;s,m s \j,mj) 
+ftf s ,m s -i{l, mi; s, m s - l\j, rrij) 

and the analogous equation for J_, 

f j>m .-i(l,mi]s,m s \j,mj - 1,) = fi, mi {l,mi + l;s,m s \j,mj) (10.23) 

as well as the normalization condition 

1 = (7,"ii,(i,«)b',mj,(Z,s)> 

Generally, Clebsch-Gordan coefficients are only non-zero for I + s > j > \l — s\. 
In the literature one also finds Wigner's 3j-symbol instead of Clebsch-Gordan 
coefficients, to which they are related by 

(Is 7 \ (-l) l ~ s + m j 

) ■= 7= ^(f,m t ;s,m fl |j,m 7 -) . (10.25) 

Wigner's 3j-symbol is more symmetric than Clebsch-Gordan coefficients. It is 
invariant under cyclic permutations of the three columns. Under an anti-cyclic 
permutation, or under simultaneous change of the sign of mi, m s and rrij, it ob- 
tains a factor of (— iy+ s +i. 

10.2.2 Tensor operators and Wigner-Eckart theorem 

Tensor operators T^ mk are defined by their transformation properties under ro- 
tations. They correspond to an irreducible representation of the rotation group. 
Equivalently, they can be defined as operators which fulfill the following commu- 
tation relations with a given angular momentum operator J, 

hm k T k>mk (10.26) 

ftfk,m k T k ,m k +l (10.27) 
hfk,m k -lTk,m k -l ■ (10.28) 



[J+, Tk,m k ] = 
[^-)^,mj = 
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Tensor operators are of great importance for atomic selection rules and allow an 
elegant derivation of the latter by means of the Wigner-Eckart theorem. 
To prove the Wigner-Eckart theorem we first consider the states 

\f^( h k)): = J2{j"^';\j, mj ;k,m k ) T k , mk \j, mj ) (10.29) 

m k ,mj 

The action of J z on these states is given by 

Jz\j",m'j(j,k)) = (f, m j\j, m fi k i m k) x 

m k ,m,j 

([Jz,T k ,m k ] +T k>mk J z )\j,mj) 

m k ,mj 

(hm k T k>mk + T k>mk hrrij)\j, rrij) 
= hm"\f,m/j(j,k)) (10.30) 

where the last step follows from the fact that Clebsch-Gordan coefficients are non- 
zero only for m k + rrij = m". The action of J + reads 

J + \j",m'j(j,k)) = Y {f i m "\h m ii k i m k) ([J+,T k ,m k ] +T kjmk J+)\j,m j ) 

m k ,m,j 

= Y k,m k ) 



x 



m k ,rrij 



( 



\ 



frfk,m k T k)7nk + rrij) + T k:mk hf jjmj \j,riy+J) 
\ =*»'„ = m 'j J 

,m k — l T k rrik 



m k ,m,j 



+h ^2 (j"> m j\j' m 3 ~~ UK^k) fj, mj -i T k , mk \j,mj) 

m k ,mj 



m k ,m,j 



+ (]", rn"\j, rrij - 1; k, m k ) fj ,m 3 -i)T k)rnk \j, rrij) 

m k ,mj 

Kfj», m » \j'',mf+l(j,k)) 



(10.31) 
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where we have used Eq. (110.231) in the second-to-last step. Analogously, we find 

for J_ 

J-\j",m'<(j,k)) = hf),,^ \j»,mj^l(j,k)) (10.32) 

The action of the angular momentum operators on the states m"(j, k)) is the 
same as on the angular momentum eigenstates m'j). Therefore, these states 
need to be proportional to each other. Since for given j" the eigenstates m") 
can be constructed by applying J + to —j"), this can also be done for the 

states m"(j, k)). For this reason, the proportionality factor between the states 
is independent from m'j and we find 

\f^k)) = (-l) 2fc ^™^ |/,m;) . (10.33) 

The symbol \T k \\j) is called reduced matrix element and contains the depen- 
dence of the proportionality factor on j", k and j. The factor of (— l) 2k / \ / 2j" + 1 
is convention. It is easy to see that the reduced matrix element can be calculated 
by 

(j"\\Tk\\j) = (-l) 2k y/2f+l E (j",rn , ;\j,m J ;k,m k )(j",m';\T k)mk \j,m J ) 

m k ,rrij 

(10.34) 

The matrix elements of tensor operators can be derived from Eq. (I10.33I ). Be- 
cause of the completeness of the angular momentum eigenstates we have 



^ (j, rrij] k, m k \j", m") m'j(j, k)) = ^ (j, mf, k, m k \j", m") x 

j",m" 

Yl 0'"' m 'j b'» m P fc > m 'k) T k,m' k \j, m 'j) 
= ^' mfe ^' m 'k) T k,m' k \j, 

= T k>mk \j,mj) (10.35) 

and thus 

0", ™'j\Tk,m k \j, mj) = ^ (j, mf, k, m k \j", m") (f, m'j\j" \ m'!(j, k)) 

j",m'! 



This is the Wigner-Eckart theorem. 



(-l) 2fc (j, my, k, m k \f, m',) (10.36) 

\ Ly + 1 
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10.2.3 Addition of three angular momenta 



When adding three or more angular momenta one has to be careful with the choice 
of basis. Consider the case F = I + L + S. If we first add J = S + L and then 
F = I + J, the procedure of the previous section would result in the basis states 
{|/, mf, (i,j(s, I)))}, i.e., eigenstates of F 2 , F z , I 2 , and J 2 . However, if we 
first add T, = I + S and then F = S + L, we implicitly would use the basis 
{|/, nif, (cr(i, s), I))}. The two sets of basis vectors are different. In a sense the 
addition of angular momenta is therefore not associative (although the addition of 
the operators is of course associative). Which of the two bases is more suitable 
depends entirely on the physical situation at hand. 

Detailed construction of the basis states: 
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Multiplication of the vectors results in 

{f,mf,(i,j(s,l))\f',m' f ,(a(i,s),l)) = (i,mi;j,mj\f,m f ) {cr,m a ;l,mi\f \ 

mi,mi,m a 

(s, m s ; I, mi\j, rrij) (i, mf, s, m s \cr, m a ) 
= (2/ + lV(2j + l)(2t7 + l) « 

m j , TTLfj 

m,i,mi,m s 

13 f \( a i r 



■it- 



mi rrij —nif J \ mi —m^ 
si j \ / i s a 



m s mi —mj I \ m^ m s —m a 
v/(2j + l)(2o-+ 1) x 

'_iy+s+i+f \ i s a 
I f j 



The last line defines Wigner's 6j-Symbol, which agrees with the so-called Racah 
symbols W(- ■ ■ ) up to a sign: 

{))]} = {-l) l+s+l+f W{i, 8,1, f; a, 3) . (10.40) 

The expression for the 6j-Symbol as a sum over Clebsch-Gordan coefficients is 
equivalent to Eq. (C.32) of Ref. 0. 



10.2.4 Matrix elements for composed systems 

The content of this section is of great importance for the actual calculation of 
multipole matrix elements for atoms and molecules. One often has to resolve the 
problem that the atomic energy eigenstates are (to a good approximation) eigen- 
states of the hyperfine-spin operator F, but the dipole matrix elements have to be 
calculated with respect to the electronic states alone (without nuclear spin /). Gen- 
erally one needs to express the reduced matrix element (/'(«', j') \ \Tjf' \\f(i, j)) of 
a tensor operator , that only acts on the total electronic angular momentum 

J, through (j'\\T®\\j). 

To achieve this one starts from Eq. (110.341) and expands, with the aid of 
Eq. dlO.211) . the hyperfine states \f,rrif) in terms of the nuclear spin eigenstates 
\i, m^ and the eigenstates \j, mj) of te total electronic angular momentum. One 
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then arrives at 

(f{i\j'W^\\f{hj)) = (-l) 2fe V2f + l Y, (f\rn)\f,m f ,k,m k )(f\m)\T^ k \f,rri f ) 

m k ,mf 

m k ,mf mi,m,j m'^ra'- 

(i, m f ; j, mj-l/, m/> (z', m-, j", m^T^I*, m*, j, m,-) 
= (-1)»^ + T ^ E X] m / ;fc J m fc ></' ) rn' / |t , X;j' ) m;> 

m,k,mf rrn,rrij m'^m'. 

Si,i'6 mi>m ^ (i, mi, j, rrijlf, m f ) (?, m'j |T^ fc \j, rrij) 
= V 2 f + E (/'' m '/l/' m /; fc ' m fc)(/'' m /l i '' m i;/' m j) x 

mi,mj,m'j 

(j'\\T^\\j) 

(i, rrii, j, rrijlf, m f) — /0 | , 5 (j> m i! fc > m fcl/' m j) 
^{^(/||Tf||i)x 

E (fi m fi k i m k\f\ m 'f)(h m i',fi m/ j\f'i m 'f) X 

m k ,mf 

mi,mj ,m'j 

(i,mi]j,mj\f,mf)(j,mj; k,m k \f ,m'j) (10.4 

The sum in the last two lines can be written as a 6j symbol by using Eq. (110.391) 
with /' = /, m'f — mf. In the sum over the Clebsch-Gordan coefficients only 
m/ is not a summation index. Therefore / in Eq. (110.391 ) corresponds to /' in 
Eq. (110.4 II) . All other relations can be derived from the position of the indices in 
the Clebsch-Gordan coefficients. One thus deduces 



i j f \ (-l)^W 



Y («>«;/. m il/'' m /) (f,nif,k,m k \f',m' f ) 



j' 

(j,mj;k,m k \j',rn'j) (i, rrn; j, rrijlf, m f ) (10.42) 



rrii,frik ,rrij 



and therefore 

(fii'J'W^WfiiJ)) = v / (2/ + l)(2f + l)^(-ir +J+fe+/ 'x 

(/ll^f ||i>{ * |, } (10.43) 
This is equivalent to Eq. (C.90) of Ref. 0. 
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10.3 Distributions 

Distributions are mathematical objects which, strictly speaking, are only defined 
if they appear in an integral. For instance, S(x — xq) is only defined through its 
effect inside the integral 

dx f(x)5(x - xo) := f(x ) 

Here, f(x) is a so-called Test function. These are functions which, speaking in 
a sloppy way, "always behave well". More formally, they are C°° functions with 
compact support, or at least they converge sufficiently fast to zero for |x| — > oo). 
Distributions are no functions. For example, S(x — x ) is not well defined as a 
function at the position x . 

The integral D(u) := J °° dr e~ luJT is also only defined in the sense of distribu- 
tions: 



/oo /*oo /*oo 

dcu D(uj)f(uj) = / dto f(to) / dr e 
-oo J~oo JO 



oo JO 

00 Pc 
du) f(uj)(ir5(uj) — i — ) 

= nm-i hm ( r + r d -m) 

The transition from the first to the second line relies on a mathematical theorem. 
The principal value in this expression means that one has to use the same small 
parameter e on both sides of the singularity of the integrand. If one calculates the 
limiting value of an integral without principal value, one has to admit different 
small parameters on both sides of the singularity. 
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